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!•  INTRODUCTION 


The  AEDC  Dust  Erosion  Tunnel  (DET)  is  currently  using  a Huels-type, 
high  voltage  arc  air  heater  to  accelerate  the  test  section  oust  cloud  for 
determination  of  the  surface-recession  characteristics  of  reentry  vehicle 
(RV)  materials  in  an  abrasive  enviroument.  The  maximum  attainable  particle 
velc'tty  is  presently  less  than  7000  ft/s.  The  use  of  hydrogen  or  helium 
instead  of  air  as  the  driving  gas  has  the  potential  of  higher  gas  enthalpies 
thus  increasing  this  upper  velocity  limit.  However,  no  data  were  available 
in  the  open  literature  on  the  performance  characteristics  of  the  Huels-type 
arc  heater  operating  on  hydrogen  or  helium  to  confirm  the  higher  gas 
enthalpies. 

The  objective  of  this  work  was  to  map  the  performance  of  a Huels-type 
arc  heater  using  gaseous  hydrogen  at  pressures  of  10,  15,  25  and  50  atm  and 
helium  at  pressures  of  25,  50,  75  and  100  atm  with  arc  currents  of  250,  320 
and  400  A.  A supersonic  nozzle  having  a 0. 162  in«  diam  throat  was  used  to 
constrict  the  flow,  and  measurements  of  the  stagnation  enthalpy,  stagnation 
pressure,  and  arc  power  were  made.  These  performance  data  were  compared 
with  data  using  air  and  with  computer  coded  correlations  derived  to  permit 
scaling  of  the  arc  heater.  Electrode  erosion  was  documented  during  these 
performance  tests.  Minor  modifications  to  the  basic  arc  heater  were  made 
to  achieve  stable  operation  over  the  full  test  matrix  with  optimum  performance 
as  a goal. 

Use  of  a Linde  N-250  arc  heater  satisfied  all  test  requirements.  This 
heater  was  less  than  one-third  the  size  of  the  AEDC  DET  heater.  MDRL  had 
previously  conducted  more  than  700  tests  on  air  using  a Linde  N-250  arc  heater 
with  nozzle  throat  diameters  ranging  from  0.142  to  0.250  in.  Unrestricted 
flow  tests  with  a 0.375  in.  diam  front  electrode  had  also  been  run.  One 
scale  factor  of  significance  with  high  voltage  arc  heaters  is  the  ratio  of 
the  nozzle  throat  area  to  the  front  electrode  gas  flow  area.  In  order  to 
maintain  the  ratio  present  in  the  DET  arc  heater,  the  N-250  required  a 
0.162  in.  diam  throat.  Since  the  heater  had  been  operated  on  air  with  throats 
near  that  size,  inspection  of  those  data  and  the  heater  specifications 
indicated  that  operating  on  air  at  pressures  from  5 to  50  atm  and  arc  currents 
from  200  to  400  A was  feasible. 

The  AEDC  DET  scale  N = 1.73  high  voltage  arc  heater  has  a nozzle  throat 
diameter  of  0.5625  in.  Operation  of  this"  scale  heater  was  satisfactory  for 
air  as  the  working  gas,  but  there  was  serious  doubt  al.out  this  configuration 
being  optimal  for  use  with  hydrogen.  The  scaling  relations  for  air  used  by 
Linde  had  proven  to  be  erroneous, ^ so  it  was  doubtful  that  whey  could  accurately 
be  applied  to  hydrogen.  Thus,  selection  of  the  optimum  scale  heater  for 
hydrogen  required  scaling  studies  as  well  as  performance  measurements. 

A minimum-energy-addition,  arc-model- scaling  computer  code  had  been 
used  successfully  by  MDRL  on  air  arc  heaters.  The  addition  of  the  proper 
thermodynamic  and  transport  properties  for  helium  and  hydrogen  plus  an 
improved  model  for  the  Huels-type  arc  heater  was  a major  part  of  this  program. 
Scaled  performance  calculations  were  then  made  with  the  improved  codes. 
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2.  TEST  APPARATUS 


The  arc  heater  used  in  these  tests  was  a nominal  250  kW  Huels-type,  high- 
voltage  unit.  It  was  mounted  and  tested  in  the  McDonnell  Douglas  Research 
Laboratories  (MDRL)  High  Impact  Pressure  (HIP)  facility.^  ' The  HIP  subsystems, 
controls,  and  data  acquisition  system  were  used  along  with  special  gas  controls. 
The  instrumentation  used  was  a combination  of  standard  and  special  equipment 
selected  for  specific  test  requirements-  Data  reduction  routines  were  written 
to  convert  the  raw  signals  into  printed  engineering  units  using  the  HIP 
computerized  data  acquisition  system. 

2.1  THE  N-250  ARC  HEATER 

Figure  1 shows  the  N-250  arc  heater  used  in  these  tests.  Two  tandem 
cylindrical  electrodes  were  separated  by  a central  gas- injection  chamber  and 
insulated  from  each  other  by  a threaded  Delrin  cylindrical  Tnsulator.  The 
5.6  in.  long  rear  electrode  (anode)  had  an  inside  diameter  of  0.625  in.  The 
front  electrode  (cathode)  had  an  inside  diameter  of  0.375  in.  The  cathode 
length  was  5 in.  in  most  tests  and  8 in.  in  a few  tests,  designa  ed  as  N-250-5 
or  N- 250-8.  Both  electrodes  and  the  center  chamber  were  water  cooled.  Four 
tubular  gas  injectors  were  used.  The  gas  was  injected  tangentially  to  the 
heater  walls  in  a counter-clockwise  direction  looking  downstream.  The  resul- 
tant radial  pressure  gradient  stabilized  the  arc  on  the  heater  centerline. 

The  injection  pressure  ratio  (inector-pressure/arc-chamber  pressure)  was 
varied  from  1.1  to  2.2  in  these  tests. 

A magnetic  field  coil  was  used  on  the  rear  electrode  to  rotate  the  arc 
termination  and  prevent  arc  transfer  to  the  rear  plug. 

A nozzle  module  was  bolted  to  the  downstream  end  of  the  heater  to 
constrict  the  flow.  The  nozzle  throat  diameter  was  0.162  in.  for  all  tests 
def/ignated  by  N-250- 5-0. 162  or  N-250-8-0. 162.  The  nozzle  design  is  discussed 
in  Section  2. 2- 

Figure  2 shows  the  N-250  arc  heater  mounted  in  the  MDRL  HIP  test  cell- 
Each  component  was  individually  water  cooled  and  monitored  for  energy  losses. 
The  differential  temperature  sensors  and  turbine  flowmeters  can  be  seen 
beneath  the  arc  heater  test  stand-  The  test  set-up  allowed  rapid  disassembly 
for  weighing  and  cleaning  the  electrodes  post-test.  Much  of  the  gas  flow 
instrumentation  was  fastened  to  the  rear  panel  on  the  test  stand,  so  that  the 
entire  basic  system  was  portable. 
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Figure  1 N'250  arc  heater 


AEDC.TR-76-25 


Figure  ^ N>250  arc  heater  mounted  in  the  HIP  facility 
2.2  ARC  HEATER  NOZZLE 


The  nozzles  used  in  these  tests  (Figure  1)  were  s^'^led  to  duplicate  the 
AEDC  DET  arc  heater  throat-to-electrodc  area  ratio. 


d*  = 0.162  in.  (3) 

The  nozzle  inlet  diameter  matched  the  electrode  diameter  (0.375  in.),  and  the 
exit  diameter  was  sized  for  proper  expansion  from  25  atm,  30,000  Btu/lb 
reservoir  conditions  to  atmospheric  pressure.  This  yielded  ^ 0.296  in.  diam, 
Mach  2.75  test  stream.  The  convergence  and  divergence  of  the  nozzle  were 
conical  for  ease  of  machining;  a short  cylindrical  throat  was  'ncorporated 
for  erosion  protection  and  improved  source  flow. 

2 

The  nozzle  throat  cooling  was  designed  to  withstand  6000  Bcu/ ft  s heat 
flux  without  nucleate  boiling  with  a coolant  flow  rate  of  30  g;jm.  The  coolant 
flow  velocity  at  the  throat  exceeded  170  ft/s,  which  was  sufficient  to  sweep 
vapor  bubbles  away  at  the  higher  heating  rates.  The  first  and  second  nozzles 
experienced  local  heating  rates  to  approximately  9800  Btu/ft^s  before  failure, 
and  the  third  nozzle  liner  experienced  heating  rates  to  6700  Btu/ft^s  without 
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failure  . The  calculated  combined  thermal  and  mechanical  stress  at  the 
throat  was  20,000  psi  at  the  design  heat  flux.  The  nozzle  throat  wall 
thickness  was  initially  0.0315  in.  All  of  the  hydrogen  testing  was  completed 
with  two  such  nozzles.  A third  nozzle  was  built  with  a 0.047  in.  throat  wall 
thickness  for  the  high-pressure  air  and  helium  tests. 

2.3  INSTRUMENTATION  AND  DATA  ACQUISITION 

The  arc  heater  was  instrumented  to  determine  the  following  parameters. 


0 

Arc  current 

0 

Anode  power  loss 

o 

Arc  voltage 

o 

Cathode-body  power  loss 

0 

Arc  chamber  pressure 

0 

Nozzle  power  loss 

o 

Gas  flow  rate 

0 

Injector  pressure 

These  parameters  were  scanned  at  100  scans/s  on  hydrogen  and  10  scans/s 
on  helium  and  air.  These  data  w^re  averaged  and  the  arc  power,  energy  balance 
enthalpy,  thermal  efficiency,  and  sonic  flow  enthalpy  were  calculated  and 
printed  out  approximately  once  each  second  of  the  test. 

Table  1 lists  the  primary  instruments  and  their  accuracy.  The  consider- 
able electrical  noise  encountered  when  operating  on  hydrogen  required  thorough 
shielding  of  all  instruments,  increased  sampling  rates  and  averaging  larger 
sampled  batches  in  order  to  acquire  acceptable  data.  No  noise  problems  were 
encountered  on  helium  or  air.  The  sampling  rates  and  batches  were  reduced 
without  increasing  the  data  scatter  significantly. 

The  arc  voltage  on  hydrogen  oscillated  significantly  and  caused  strong 
random  fluctuations  in  the  printed  data.  Several  attempts  to  reduce  these 
oscillations  were  made  by  changing  the  sampling  rate  and  the  averaging  batch 
size;  however,  voltage  spil^.es  persisted.  A satisfactory  technique  was  derived 
by  filtering  the  signal  with  a digital  voltmeter  isolation  amplification 
network.  This  was  not  necessary  for  the  helium  or  air  te^sts. 

The  arc  chamber  pressure  was  measured  on  the  side  wall  of  the  center  body. 
Previc't'v  tests  have  shoTO  that  this  value  may  be  higher  than  the  stagnation 
pressure  at  the  nozzle  entrance.  This  may  partially  explain  some  of  the 
deviations  found  between  the  sonic  flow  and  energy  balance  enthalpies. 

The  anode  and  noz*^le  had  individual  coolant  differential  temperature 
sensors  and  water  flowmeters  for  measurdng  the  power  losses  to  these  components. 
The  cathode  and  body  coolant  were  supplied  through  the  same  port,  and  since 
the  body  power  losses  ar^*  minimal  it  was  decided  to  combine  them  with  those 
of  the  cathode. 

The  hydrogen  flow  rates  were  measured  using  a sonic  nc  j.e  flowmeter  with 
a 0.055  in.  diam  throat.  The  flowmeter  was  calibrated  by  tne  MCAIR  Bureau 
of  Standards  at  air  flow  rates  from  0.004  to  0.123  Ib/s.  The  orifice  co- 
efficient on  air  wai  derived  from  28  data  oolnts  and  was  accurate  within 
+ 1.5%  with  a confidence  level  of  99.5%.  Conversion  of  the  air  calibration 
to  hydrogen  was  accomplished  using  the  techniques  given  in  References  6 and  7. 
Godt'  found  for  critical  nozzle  flowmeters  having  throat  diameters  Irom  0.0308 
to  0.0707  in.  that  the  relative  values  of  the  nozzle  discharge  coefficients 
demonstrate  that  a flowmeter  calibrated  on  one  gas  can  be  used  to  mcjasure  flow 
rates  of  the  other  gases  within  + D. 1%. 
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The  helium  flow  ratet^  were  measured  using  a 0,070  in*  diam  sonic  nozzle 
flowmeter  for  arc  pressures  to  75  atm.  Because  of  gas  supply  limitations 
with  conventional  storage  bottles,  the  100  atm  test  flow  rates  were  based  on 
a calibration  of  the  four  0.0625  in.  diam  injectors.  The  injectors  were 
calibrated  on  helium  in  series  with  the  0.070  in.  diam  sonxc  nozzle.  The 
subsonic  mass  flow  rate  followed  the  general  equation  for  subcritical  real 
gas  flow  through  a thin-plate  orifice  within  + 0.56%. 


The  air  flow  rates  were  measured  using  the  0.055  in.  diam  sonic  nozzle 
for  arc  pressures  from  50  to  100  atm  and  0.043  in.  diam  sonic  nozzle  for  arc 
pressures  from  10  to  25  atm. 

Tests  were  made  on  all  flow  systems  to  determine  if  the  flowmeters  were 
choked.  Corrections  were  necessary  on  some  of  the  high  current  hydrogen 
tests  where  the  injector  pressure  exceeded  73%  of  the  flowmeter  pressure. 

All  primary  instrumentation  signals  were  recorded  on  discs  in  the 
Research  Data  Acquisition  System  (RDAS).  Data  for  the  hydrogen  tests  were 
scanned  at  100  scans/s;  data  for  the  air  and  helium  tests  were  scanned  at 
?0  scans/s.  The  disc- recorded  information  was  transferred  to  tape  post- test, 
processed  into  engineering  units,  and  printed  in  the  format  discussed  in 
Appendix  A. 

2,4  ARC  HEATER  SUBSYSTEMS 


The  arc  heater  v/as  tested  in  the  MDRL  HIP  facility.  The  water,  power, 
and  air  subsystems  were  more  than  sufficient  for  this  scale  heater.  The 
standard  HIP  arc  heater  CMOC-300)  is  four  times  the  size  of  the  N-250,^  ' 

A 13  0 ballast  resistor  was  necessary  for  this  program.  The  resistor 
consisted  of  60  water-cooled  stainless  steel  0.25  in.  diam  tubes  with  0.020  in. 
wall  taickness;  each  tube  was  10  ft  long.  The  tubes  were  manifolded  in  two 
banks  of  30,  and  each  bank  was  separately  cooled.  Connections  of  0.9,  1,8, 

3.3,  7 and  13  were  possible.  The  maximum  current  capability  was  in  excess 
of  400  A.  This  balxast  resistor  provided  the  immediate  voltage  reserve 
required  by  an  arc  excursion  without  the  inherent  time  delay  of  the  power 
supply.  The  7 Q ballast  connection  was  used  on  all  reported  hydrogen, 
helium  and  air  tests. 

Cooling  water  was  supplied  by  a threr^stage  centrifugal  pump  capable  of 
1200  gpm  at  600  psi.  The  normal  flow  rates  were  23  gpm  to  the  nozzle,  34  gpm 
to  the  cathode-body  and  20  gpm  to  the  anode.  The  ballast  resistor  required 
120  gpm  and  the  spin  ceil,  2 gpm. 

Power  was  supplied  by  two,  three-phase,  full-wave  rectifier  modules 
with  saturable  reactor  current  controls.  The  units  were  connected  in  series 
for  an  open  circuit  voltage  of  10.2  kV.  The  maximum  power  capability  of 
these  two  units  was  5 MW  for  6 rain.  Additional  identical  power  units  were 
connected  in  series  without  activation  to  utilize  the  inductance  of  th^ir 
large  smoothing  reactors  (8  mH).  The  number  of  additional  units  varied  from 
1 to  3.  Six  power  units  were  available. 
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The  hydrogen  and  helium  gas  supplies  consisted  of  multiple  type-lA 
cylinders  manifolded  as  shown  in  Figure  3.  The  gas  flowed  through  an  on-off 
solenoid  valve,  a check  valve,  and  a dome-loaded  regulator  to  the  primary 
flowmeter,  then  through  a second  check  valve  to  the  arc  heater  injectors. 
Argon  or  nitrogen  was  used  to  start  the  heater  during  the  hydrogen  tests, 
and  the  dome  regulator  load  gas  was  either  argon  or  nitrogen.  The  normal 
number  of  hydrogen  cylinders  was  eight.  As  many  as  20  helium  cylinders  were 
manifolded  for  the  high  pressure  tests. 


Argon-nitrogen 


i 


Figure  3 Gas  flow  system  for  arc  heater  operation  on  hydrogen 
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Air  was  supplied  from  the  HIP  system  which  is  capable  of  6000  psi  supply 
pressure  and  Tlow  rates  in  excess  of  4 Ib/s.  The  air  was  routed  to  the  same 
regulator  and  flowmeter  system  used  for  hydrogen  and  helium. 

The  RIP  Research  Data  Acquisition  System  (RDAS),  used  for  all  tests, 
consists  of  a Varian  model  620/f-100  computer  with  32,000  (16-bit)  words  of 
core  memory,  a disc  memory  with  a capacity  of  2.34  million  (16-bit)  words, 
a 7-track  magnetic  tape  unit  with  556  bpi,  an  electrostatic  printer/plotter, 
a card  reader  and  a teletype.  The  system  currently  has  48  data  channels  and 
is  expandable  to  256.  The  RDAS  has  a maximum  sequential  scanning  capability 
of  50,000  channels/s  with  a 100  ns  software  overhead  between  records.  It  can 
also  be  programmed  to  scan  selected  channels  on  an  individual  basis  at  a 
maximum  rate  of  33,000  channels/s  with  a 60  us  software  overhead  between 
records.  This  mode  of  scanning  can  be  used  for  the  acquisition  of  temperature- 
time data  to  calculate  heat  flux  Incident  on  a model  subjected  to  the  high- 
temperature  plasma  environment. 

Binary  information  including  sign  and  parity  checks  is  the  primary 
output.  The  format  is  suitable  for  digital  editing,  sorting,  computing, 
and  plotting  by  electronic  computers.  The  data  are  packed  in  IBM  format,  on 
7-channel,  magnetic  tape  at  a packing  density  of  556  bits/in. 

The  overall  system  accuracy,  from  the  signal  input  to  the  digital  tape 
output,  is  approximately  + 0.10%.  This  accuracy  takes  into  account  all  items 
Including  the  multiplexing  unit. 

The  digital  readout  is  available  for  the  display  of  any  given  input 
channel  for  signal  conditioning  set-up  and  check-out,  as  well  as  viewing 
any  channel  during  an  arc  heater  test.  This  readout  was  used  to  monitor 
the  arc  pressure  for  these  tests. 
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3.  TEST  TECHNIQUE  AND  MATRIX 


The  techniques  used  to  operate  the  arc  heater  in  this  program  were 
basically  those  developed  for  the  HIP  facility  over  the  past  decade,  plus 
special  safety  procedures  for  hydrogen  operation.  The  test  matrix  included 
operation  at  three  arc  currents  at  pressures  from  10  to  50  atm  on  hydrogen, 
25  to  100  atm  on  helium,  and  10  to  100  atm  on  air. 

3.1  ARC  HEATER  OPERATION 


Preparation  of  the  arc  heater  for  testing  included  careful  assembly  and 
mounting  followed  by  a full-flow  water  check,  static  gas-pressure  check,  and 
an  insulator  resistance  check.  The  procedures  outlined  in  Reference  8 were 
followed  where  applicable. 

The  normal  arc  heater  starting  procedure  for  a hydrogen  test  was  as 
follows: 

1.  Observe  the  hydrogen  detector  alarm  to  ensure  a safe  test  cell 
hydrogen  level. 

2.  Turn  on  low-speed  exhaust  fan  in  the  HIP  test  cell. 

3.  Visually  check  the  gas  flow  system  and  arc  heater. 

4.  Open  the  argon  or  nitrogen  gas  supply  bottles  (hydrogen  bottles 
closed) . 

5.  Check  the  functioning  of  all  flow  system  components  to  ensure 
proper  control  and  operation  using  argon. 

6.  Follow  applicable  HIP  facility  operating  procedures  in  preparation 
for  the  test. 

7.  When  all  subsystems  including  the  data  acquisition  system  are  ready 
for  the  test,  turn  on  the  high-speed  test-cell  exhaust  fan.  Check 
to  see  that  all  interlocks  are  satisfied. 

8.  Close  the  argon,  nitrogen  and  jiydrogen  on-off  valves. 

9.  Manually  open  the  hydrogen  bottle  valves  and  check  the  supply 
pressure  in  the  bottles. 

10.  Observe  the  flowmeter  pressure  with  the  main  solenoid  vnive  closed 
to  ensure  no  leakage  through  the  system. 

11.  Open  the  argon-nitrogen  on-off  valve  and  adjust  the  regulator  to 
the  desired  starting  flow  rate, 

12.  Turn  on  the  power  supplies.  (For  nitrogen  starts,  engage  striker.) 
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13.  Stabilize  the  arc  heater  on  argon  or  nitrogen,  then  open  the  hydrogen 
on-off  valve  and  close  the  argon-nitrogen  on-off  valve. 

14.  Stabilize  the  arc  heater  on  hydrogen  and  acquire  the  desired  data 
by  regulating  the  hydrogen  gas  flow  rate  and  arc  current.  Hold 
each  point  for  30  s minimum. 

15.  Turn  off  the  power  supplies.  (Loss  of  current  automatically  switches 
back  to  argon  or  nitrogen. ) 

16.  Aij.ow  the  high-speed  fan  to  continue  blowing  while  the  hydrogen 
bottles  are  closed  manually  and  the  system  is  purged  \d.th  argon  or 
nitrogen. 

17.  Observe  the  hydrogen-detector  alarm  to  ensure  a safe  test  cell 
hydrogen  level  prior  to  entering. 

18.  Secure  all  systems  and  reduce  data. 

Rigid  safety  precautions  and  procedures  were  followed  for  hydrogen  arc 
heater  operation. 

All  gaseous  hydrogen  was  stored  in  an  eight-bottle  rack  along  the  outside 
of  the  building.  Bottled  hydrogen  was  acquired  from  the  pan  stock  crib  on 
an  as-needed  basis.  No  more  than  eight  full  bottles  were  allowed  on  the 
premises  at  one  time.  No  smoking  was  allowed  in  the  vicinity  of  the  HIP 
test  cell  or  hydrogen  storage. 

All  personnel  were  familiar  with  the  industrial  safety  procedures  for 
handling  hydrogen  and  adhered  to  them  at  all  times. 

The  hydrogen  supply  system  was  constructed  using  valves,  fittings,  and 
pipe  made  of  materials  that  were  compatible  with  hydrogen.  Connections  were 
minimized  to  avoid  leaks. 

The  complete  hydrogen  supply  system  was  pressure  checked  at  least  weekly 
using  2500  psi  helium.  The  system  was  pressurized  for  at  least  10  min,  and 
the  pressure  decay  had  to  be  less  than  1 psi/min.  Leaks  greater  than  this 
rate  were  located,  corrected,  and  rechecked  prior  to  operation. 

All  hydrogen  bottles  were  closed  manually,  and  the  system  was  purged 
with  inert  gas  when  the  next  test  was  to  be  more  than  one  hour  later  and/or 
when  the  testing  was  complete  for  the  day. 

The  hydrogen  bottles  were  returned  to  the  pan  stock  crib  when  further 
testing  was  not  scheduled  within  a one  week  period. 

Arc  heater  operation  on  helium  did  not  require  any  special  precautions. 
The  arc  was  started  on  helium  by  turning  on  the  power  supply  with  a low  flow 
rate  (0.005  Ib/s).  The  desired  arc  pressure  and  arc  current  were  then  set 
and  held  for  at  least  30  s.  Usually  three  arc  current  levels  were  tested  at 
a fixed  pressure  during  a single  test. 
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Operation  on  air  required  an  argon  start.  As  the  arc  current  was 
initiated  the  main  air  valve  opened.  When  the  air  supply  pressure  exceeded 
the  argon  supply  pressure,  a check  valve  closed  off  the  argon.  The  test 
procedure  after  starting  was  the  same  as  for  helium. 

Following  most  of  the  tests,  the  electrodes  v^ere  removed,  weighed, 
inspected,  cleaned,  reweighed,  and  re-installed  or  replaced*  Section  4 
discusses  the  resultant  contamination  measurements.  The  nozzle  throat 
diameter  was  measured,  and  the  nozzle  was  cleaned  for  the  next  test. 

3.2  DATA  ACQUISITION 


The  experimental  data  were  acquired  using  the  Instrumentation  and 
subsystems  described  in  Sections  2.3  and  2,4.  The  procedure  consisted  oi.  the 
following  steps: 

1.  Load  acquisition  system  into  the  computer  via  the  boot-strap  loader. 

2.  Enter  by  teletype  the  date,  run  number,  number  of  channels  (21) 
scanned,  and  scan  interval  (0,01  s for  hydrogen  and  0,1  s for 
helium  and  air). 

3.  Start  the  system  acquisition  and  verify  that  it  is  functioning 
during  the  test, 

4.  Stop  the  system  acquisition  post-test  and  load  the  vortex  system 
into  the  computer. 

5.  Transfer  disc  data  to  magnetic  tape. 

6.  Call  the  engineering  data  process  program  and  print  out  data. 

7.  Pack  as  many  tests  as  possible  on  a single  magnetic  tape  when 
time  permits. 

In  addition  to  the  primary  data  acquisiton  system,  it  was  also  the 
practice  to  record  the  arc  voltage,  arc  current,  system  pressures,  system 
water  flow  rates,  and  arc  heater  description  on  a separate  data  sheet. 
Electrode  weights,  arc  track  locations,  and  other  pertinent  data  were  also 
recorded. 

A complete  list  of  Instruments  and  their  calibrations  was  maintained 
and  updated  as  changes  were  made  (Table  1) . 
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TABLET  INSTRUMENTATION 


Parameter 

Instruments 

Calibrated 

range 

Overall 

system  accuracy 
(%) 

Test 

no. 

Mfg 

Model 

Arc  voltage 

G.E. 

JTG-1 

0 - 3000  V 

0.6 

35- 

109 

Arc  current 

G.F. 

JDC-1 

0 - 900  A 

1.5 

35- 

109 

Arc  pressure 

CEC 

326-1 

0 - 1500  psia 

0.6 

35- 

64 

CEC 

326-1 

0 • 1500  psia 

0.6 

65- 

86 

CEC 

326-12 

0 - 1500  psia 

0.3 

87- 

109 

Flowmeter  pressure 

Statham 

PG285TC 

0 - 5000  psig 

1.0 

35- 

86 

CEC 

326-8 

0 - 3500  psia 

0.2 

87- 

109 

Injector  pressure 

CEC 

313-1 

0 • 5000  psig 

0.6 

35- 

86 

CEC 

326-12 

0 • 2500  psia 

0.3 

87- 

109 

Gas  temperature 

MDC 

I/C 

0 - 500°F 

0.8 

35- 

109 

HY-CAL 

RTS4135 

0-  125°F 

0.5 

35- 

86 

HY-CAL 

RTS  4135 

0-  l25°F 

0.5 

87- 

109 

Water  temperature’ 

MDC 

C/A 

0 - 600°F 

0.7 

35- 

109 

Anode  AT 

DELTA-T 

100 

0-  100°F 

1,0 

35- 

109 

Cathode-body  AT 

DELTA-T 

100 

0-  100°F 

1.0 

35- 

109 

Nozzle  at 

DELTA-T 

100 

0-  100°F 

1.0 

35- 

109 

Anode  w 

Brooks 

HP-16M 

0 - 50  gpm 

1.5 

35- 

109 

Cathode-body  w 

Brooks 

HP-16M 

0 - 50  gpm 

1.5 

35- 

109 

Nozzle  w 

Brooks 

HP-16M 

0 • 50  gpm 

1.5 

35- 

109 

Body  w 

Brooks 

HP-8M 

0 - 6 gpm 

1.5 

35- 

109 

3.3  TEST  MATRIX 


Table  2 shows  the  test  matrix  for  this  program  and  lists  the  test 
numbers  for  each  condition.  Each  test  condition  was  repeated  except  the 
240  and  400  A points  at  75  and  100  atm  on  air.  The  arc  current  levels  were 
established  by  stability  limitations  and  the  specification  of  a minimum 
increment  of  25%.  These  currents  correspond  to  scaled  (I/N)  DET  arc  heater 
currents  of  800,  1100  and  1400  A which  are  moderate  to  high  for  that  heater. 
The  enthalpies  are  thus  representative  of  the  peak  values  attainable  with  a 
Huels-t)rpe  heater  such  as  the  one  used  in  the  DET. 

Numerous  additional  tests  were  made  on  all  three  gases,  and  some  of 
these  data  are  included  in  the  analysis  section  of  this  report.  The  tests 
shown  in  Table  2 represent  the  best  data  acquired  at  the  primary  conditions 
of  interest. 

Several  tests  were  made  in  the  early  part  of  the  program  on  nitrogen. 
Specifically,  tests  3 to  54  included  stabilization  on  nitrogen  and  then  a 
switch  to  hydrogen.  The  nitrogen  data  were  not  of  primary  interest  and  thus 
were  not  fully  reduced  nor  analyzed. 

The  striker  start  on  nitrogen  was  initially  thought  to  be  necessary  to 
minimize  the  transient  from  the  starting  gas  to  hydrogen.  This  procedure 
was  later  proven  unnecessary  when  argon  starts  were  successful. 
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TABLE  2 ARC  HEATER  OPERATION  ON  HVDROGEN.  HELIUM  OR  AIR: 
MATRIX  TEST  NUMBER/s 


Arc 

current 

(A) 

Arc  pressure  (atm) 

) 

10 

15 

25 

50 

75 

100 

240 

37,38 

59,60 

56,57 

_ 

— 

1 

Hydrogen 

310 

37,38 

59,60 

56,57 

64,74 

- 

- 

i 

400 

37,38 

59,60 

56,57 

- 

- 

\ 

240 

77,79 

75,76 

78,84 

07 

93, 

108, 

Air 

320 

77,80 

75,76 

78,83 

9'>  J7 

93,94 

108,109 

400 

77,80 

75,76 

78,84 

107 

93, 

108, 

\ ^ 

240 

— 

— 

95,106 

0-39 

100,104 

103,105 

i 5 

Helium 

320 

— 

— 

95,106 

98,99 

100,104 

103,105 

1 

400 

- 

- 

95,106 

98,99 

100,104 

103,105 

i 1 

! 

\ 

i 


t 
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4.  ARC  HEATER  PERFORMANCE  CHARACTERISTICS 


The  absolute  and  relative  performance  of  a Huels-t^pe  arc  heater 
operating  on  hydrogen,  helium  or  air  is  discussed  in  this  section.  The  same 
basic  arc  heater  with  only  minor  changes  was  operated  on  all  three  gases  at 
the  same  arc  currents  and  overlapping  arc  pressures.  The  primary  performance 
parameters  discussed  herein  include  the  gas  stagnation  enthalpy  and  pressure, 
the  arc  voltage,  the  energy  losses,  the  electrode  erosion,  the  resultant 
test  stream  contamination,  and  the  inferred  arc  length  required  for  scaling 
purposes. 

The  performance  on  each  gas  is  discussed  separately  in  the  same  format; 
comparisons  are  made  at  the  end  of  this  section, 

4. L HYDROGEN  ARC  CHARACTERISTICS 


The  arc  heater  was  operated  on  hydrogen  at  arc  currencs  of  240,  310  and 
400  A at  pressures  of  10,  15  and  25  atm  and  at  310  A and  50  atm.  The  arc 
parameter  values  were  significantly  different  from  those  operating  on  air  or 
helium;  yet  the  general  trends  were  similar.  The  arc  fluctuations  were  quite 
strong  on  hydrogen  requiring  additional  ballast  and  Inductance  in  the  power 
circuit.  The  hydrogen  enthalpies  were  as  much  as  seven  times  those  of  air 
and  four  times  those  on  helium,  while  the  hydrogen  arc  voltages  were  a factor 
of  two  higher  than  those  for  the  other  two  gases. 

Figure  4 shows  the  arc  heater  firing  on  hydrogen.  At  240  A current  and 
25  atm  pressure  the  exhaust  jet  was  a deep  blue  at  the  nozzle  exit  changing 
to  white  with  a light  green  border  in  the  subsonii  portion  of  the  plume  when 
the  hydrogen  was  burning.  The  supersonic  portion  of  the  exhaust  jet  appeared 
more  stable  than  the  arc  current-voltage  oscillations  indicated. 

A prolonged  test  condition  could  not  be  achieved  without  the  use  of 
approximately  16  mH  additional  Inductance  and  7 U ballast  resistance  in 
series  with  the  arc.  The  arc  current  oscillations  without  these  stabilizing 
circuit  components  were  significantly  stronger,  and  in  some  cases  the  arc 
was  totally  unstable.  At  10  atm  and  240  A,  the  arc  current  oscillation  (as 
recorded  on  an  oscillograph)  was  reduced  from  + 23%  to  + 13%  with  the  addition 
of  the  ballast  resistor.  No  major  instabilities  were  encountered  using  the 
added  inductance  and  ballast  resistance. 

Table  A-1  in  Appendix  A presents  a summary  of  the  hydrogen  data  acquired. 
The  tabulated  values  were  extracted  from  the  printed  data  sheets.  Fach  point 
represents  a stable  test  condition  where  sufficient  t^rae  was  allowed  for 
equilibration  of  the  major  parameters,  especially  the  differential  temperature 
sensors. 

The  nozzle  throat  became  with  a thin  layer  of  copper  melt  flow 

during  a test;  the  thickness  of  this  layer  was  dependent  on  run  time,  arc 
pressure  and  the  arc  current  level.  As  the  throat  size  decreased,  the  arc 
voltage  decreased  and  reductions  were  necessary  in  the  hydrogen  flow  rate  to 
maintain  a fixed  arc  pressure.  In  one  test  (64)  the  nozzle  throat  diameter 
decreased  from  0.162  in.  to  0. 14A  in.  Normally  the  diameter  reduction  was 
)ess  than  3%. 
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Fioure  4 Arc  heale:  operating  on  hydrogen  (25  atm;  240  A) 

4* 1*1  Hydrogen  Enthalpy 

The  enthalpy  levels  achieved  on  hydrogen  are  shown  In  Figure  5 as  a 
function  of  arc  current  and  pressure.  Also  shown  are  the  enthalpy  levels 
predicted  by  the  hydrogen  scaling  algorithm  (HYARC)  described  in  Section  5. 
The  experimental  enthalpy  was  determined  from  an  energy  balance  on  the  arc 
heater  combined  with  the  measured  hydrogen  flow  rate. 

The  enthalpy  varied  from  27,700  to  41,600  Btu/lb  over  the  range  tested. 
In  general,  the  enthalpy  characteristics  were  as  predicted  except  at  the  low 
pressure  and  highest  arc  current.  For  these  conditions,  a problem  was 
experienced  where  the  erosion  of  the  entrance  to  the  front  electrode  was 
quite  severe  and  the  gas  vortex  was  disrupted,  resulting  in  a shorter  arc 
length  and  lower  arc  voltage  than  predicted.  The  gradual  decay  of  enthalpy 
as  the  pressure  is  increased  at  a fixed  arc  current  is  normal  for  a Huels- 
type  arc  heater.  The  increase  in  arc  voltage  as  pressure  Increases  does  not 
compensate  for  the  higher  gas  flow  rate  and  energy  losses,  and  thus  the 
enthalpy  dro^-j  gradually. 

The  repeatability  of  the  energy  balance  enthalpy  improved  with  increased 
presr.ure.  The  enthalpy  at  10  atm  was  repeated  within  approximately  + 11%  of 
the  average  value,  and  at  15  atm  it  was  repeated  within  approximately  + 7%. 

At  25  atm,  the  repeatability  imprc’ed  with  increased  arc  current.  At  240  A 
and  25  atm,  the  enthalpy  was  repeated  within  + 6%  of  the  average  value, 
whereas  at  400  A it  was  the  same  within  measurement  error  in  both  tests.  At 
50  atm  the  enthalpy  was  repeated  within  + 4%  of  the  average  value. 
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Figure  5 Hydrogen  enthalpy  performance  of  the  N-?50  arc  heater 

The  sonic  flow  enthalpy  was  computed  for  each  data  point*  The  enthalpy 
is  useful  provided  the  radial  variation  of  gas  properties  is  not  significant 
through  the  nozzle  throat  and  the  throat  size  remains  fixed.  The  relation 
used  was  derived  from  the  basic  sonic  flow  equation  and  equilibrium  thermo- 
dynamic properties  of  hydrogen.  These  properties  were  incorporated  in  a 
subroitine  to  the  data  reduction  program.  The  sonic  flow  ratio  (m/A*  Pq)  was 
calculated  for  equilibrium  hydrogen  at  pressures  from  5 to  50  atm  and 
temperatures  from  5000  to  10,000®R.  These  values  were  then  curve-fitted  to 
within  + 5%  resulting  in  the  following  equation. 

A*  P 
m 


This  equation  used  the  experimental  flow  rate,  pressure,  and  original  throat 
size  to  compute  the  stagnation  temperature  v;hich  in  turn  yielded  the  sonic 
flow  enthalpy  from  a subroutine. 


Comparisons  between  the  energy  balance  enthalpies  and  sonic  flow 
enthalpies  are  shown  in  Figures  6 through  8.  At  the  low-current  level,  the 
agreement  is  fair,  but  as  the  arc  current  is  increased  and  test  time  elapses, 
the  two  values  diverge  significantly.  This  discrepancy  was  a strong  indica- 
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tion  that  the  nozzle  throat  area  was  being  reduced  by  resolidified  electrode 
melt.  Post-- test  inspection  of  the  nozzles  verified  the  throat  area  reduction* 
Since  the  calculated  sonic  flow  temperature  varied  as  the  nozzle  throat 
diameter  to  the  2.46  power,  a diameter  reduction  from  0.162  in.  to  0,144  in. 
results  in  a 34%  increase  in  calculated  sonic  gas  temperature  and  a typical 
Increase  from  30,000  Btu/lb  to  70,000  Btu/lb  in  the  calculated  sonic  flow 
enthalpy.  Thus  a factor  of  2.33  increase  in  calculated  sonic  flow  enthalpy 
could  result  from  a 0.009  in.  thick  layer  of  melt  flow  in  the  nozzle  throat. 
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In  test  61,  n steady  arc  current  of  240  A and  a pressure  of  25  atm  was 
held  for  nearly  3 min.  In  this  test,  the  energy  balance  enthalpy  was 
32,000  Btu/lb  at  the  beginning  and  30,000  Btu/lb  at  the  end.  The  sonic  flow 
entiialpy  rose  from  22,000  Btu/lb  to  85,000  Btu/lb.  This  cliange  was  indicative 
of  a final  throat  diameter  of  0.132  in.  which  is  somewhat  less  than  the  0,144  in. 
minimum  measured. 

^ Hydrogen  Arc  Voltage 

The  arc  voltage  operating  on  hydrogen  was  nearly  a linear  function  of  the 
hydrogen  flow  rate  (Figure  9)  and  nearly  an  inverse  square  root  function  of 
the  arc  current  (Figure  10).  The  arc  voltage  increased  from  550  to  3600  V 
as  the  hydrogen  flow  rate  increased  from  0.004  to  0.02  Ib/s.  The  increased 
flow  rate  caused  an  increased  arc  pressure.  The  arc  voltage  increased 
linearly  with  pressure  as  shown  in  Figure  11  where  the  correlation  equation 
was  found  to  be: 


V = 67-2  + 140  (I  = 300  A),  (6) 


This  equation  along  with  the  arc  current  dependence  in  Figure  10  was  used  in 
the  basic  HYARC  program  to  establish  an  arc  length  correlation.  The  results 
of  that  correlation  are  shown  in  Figure  12  where  the  arc  voltages  from  the 
complete  test  matrix,  including  peak  values  during  each  test,  are  compared 
with  the  improved  HYARC  output.  This  comparison  is  discussed  in  Section  5. 


Figure  9 Effect  of  hydu.^en  fiow  rate  on  arc  voltage 
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Figure  10  Effect  of  arc  current  on  hydrogen  arc  voltage 


The  arc  voltage  decreased  during  a test  for  several  reasons.  The  nozzle 
throat  area  reduction  was  accompanied  by  a hydrogen  flow  rate  reduction  (to 
maintain  a fixed  pressure)  which  in  turn  reduced  the  arc  extension  forces 
and  the  arc  column  cooling,  resulting  in  lower  arc  voltage.  In  addition  the 
roughened  electrodes  in  the  restrike  region  created  high  electric  field  points 
that  retained  the  arc  terminations  longer,  thus  lowering  the  average  arc 
length  and  voltage.  Finally,  the  roughened  entrance  region  disturbed  the  gas 
vortex  resulting  in  more  rapid  column  spreading  which  caused  early  restrike 
and  short,  low-voltage  arcs. 

Oscillations  in  the  arc  voltage  were  analyzed  early  in  the  testing  to 
eliminate  noise  on  the  data  record.  A storage  osc-'lloscope  was  used,  and 
noise  frequencies  of  150  kHz  with  a modulation  of  2.5  kHz  were  Identified. 

The  signal-to-noise  ratio  was  as  low  as  0.5.  An  isolation  amplifier  was  sub- 
sequently used  to  eliminate  the  noise  from  the  data  record.  The  oscillograph 
traces  of  the  arc  voltage  indicate  oscillations  in  excess  of  + 12%  of  the 
average  value  at  50  atm  and  300  A. 

Inspection  of  the  oscillograph  traces  at  high  '5peed  revealed  a sawtooth 
arc  voltage  resulting  from  the  arc  rotation  and  downstream  extension  followed 
by  abrupt  shunting  upstream.  This  is  typical  of  Hi;els-type  arc  heater  voltage 
behavior. ^ 
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Arc  pressure  (atm) 
Figure  12  Hydrogen  arc  voltages 


4. 1. 3 Energy  Losses  on  Hydrogen 

The  energy  losses  to  the  anode,  cathode-body,  ana  nozzle  operating  on 
hydrogen  are  given  in  Table  A-1.  Anode  (rear  electrode)  losses  varied  from 
11, 7%  to  23.67  of  the  total  loss  over  the  full  matrix  of  test  conditions. 

The  magnitude  of  the  anode  losses  Increased  with  arc  current  and  pressure  as 
shown  in  Figure  13.  The  anode  percent  of  the  total  losses  Increased  with 
arc  current  at  10  atm  pressure  but  decreased  at  15  and  25  atm.  The  anode 
percent  of  the  total  losses  Increased  with  pressure  from  12.7%  to  13.9%  at 
the  310  A level.  As  the  gas  flow  rate  increased  (higher  pressure),  the  rear 
arc  termination  moved  further  back  exposing  a larger  area  to  arc  radiation. 
Since  there  was  cold  reverse  flow  along  the  rear  electrode  wall,  the  losses 
x^ere  predominantly  by  radiation.  The  location  of  the  rear  arc  tracks  combined 
with  the  anode  losses  indicate  radiative  heating  rates  from  845  to  1945 
Btu/ft^s  at  400  A - 10  atm  and  310  A - 50  atm  respectively.  These  values 
were  correlated  xd.th  the  HYARC  radiation  flux  calculation  as  discussed  in 
Section  5. 

The  energy  losses  to  the  cathode-body  varied  from  53.5%  to  71.5%  of  the 
total  losses#  The  absolute  cathode-body  losses  increased  with  arc  current 
and  pressure  as  shown  in  Figure  13.  Above  15  atm  pressure  the  percent  loss 
attributable  to  the  cahtode  increased  x^ith  arc  current  and  pressure.  The 
predictions  shox^  by  HYARC  are  discussed  in  Section  5.  In  general,  the  data 
correlate  well. 
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Arc  pressure  (atm) 


Figure  13  Electrode  energy  losses  on  hydrogen 

Separate  body  loss  measurements  were  not  made  but  a reasonable  estimate 
is  that  they  comprise  6%  to  7%  of  the  total  losses  when  considering  that 
the  body  coolant  flow  includes  the  entrance  corner  and  upstream  cathode 
flange  losses  plus  the  body  losses. 

An  estimate  of  the  experimental  convective  heating  was  obtained  by 
comparing  cathode-body  losses  using  an  8 in.  long  cathode  and  a 5 in.  long 
cathode  for  the  same  arc  current  and  pressure.  Thus,  the  losses  in  the  last 
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3 in,  of  the  longer  cathode  are  primarily  convective  since  the  arc  attachment 
is  contained  within  the  first  5 in.  The  experimental  heat  flux  to 
the  last  3 in.  of  the  cathode  was  940  Btu/ft'^s  at  25  atm  and  240  A.  The  HYARC 
predicted  convective  heat  flux  was  1150  Btu/ft^s.  The  comparative  fluxes  at 
15  atm  and  240  A were  920  Btu/ft^s  measured  and  837  Btu/ft^s  predicted.  A 
discussion  of  the  convective  correlation  used  is  given  in  Reference  7. 

The  nozzle  (Section  2.2)  energy  losses  varied  from  13.8%  to  25%  of  the 
total  energy  loss.  In  contrast  to  the  anode  and  cathode-body,  the  nozzle 
losses  remained  relatively  constant  with  increased  arc  current  at  pressures 
^ 15  atm.  The  nozzle  losses  increased  with  hydrogen  flow  rate  to  the  0.8 
power.  They  correlated  within  + 20%  using  the  equation  shown  in  Figure  14 
as  a function  of  flow  rate  and  enthalpy  difference  between  the  bulk  gas  and 
the  nozzle  wall.  This  equation  was  used  in  the  HYARC  computer  code.  An 


Hydrogen  flow  rate  (Ib/s) 

Figure  14  Correlation  of  nozzle  heat  transfer  on  hydrogen 
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estimate  was  made  of  the  peak  throat  heat  flux  encountered.  A Bartz^^  heat 
flux  distribution  was  assumed  throughout  the  nozzle  and  integrated  over  the 
surface  area.  The  experimental  nozzle  loss  was  then  divided  by  the  integrated 
heat-flux-ratio/area  product  to  determine  the  nozzle  throat  flux.  The 
estimated  peak  flux  of  9800  Btu/ft^s  occurred  at  50  atm  and  310  A.  One  30  s 
test  was  completed  (64)  with  this  heat  flux,  but  the  nozzle  throat  failed  in 
the  second  (71)  and  third  (74)  tests.  The  duration  of  the  third  test  was 
sufficient  to  obtain  repeatability  data. 

Figure  15  shows  the  nozzle  after  test  57  at  25  atm  pressure.  There  was 
some  indication  of  nucleate  boiling  despite  the  fact  the  estimated  throat 
flux  was  less  than  6000  Btu/ft^s  for  this  test.  The  white  area  on  the  hot 
backwall  of  the  entrance  region  to  the  nozzle  was  caused  by  precipitation  of 
water  impurities  in  a poorly  cooled  region.  The  two  nozzle  failures  in 
later  tests  occurred  at  the  juncture  of  the  straight  throat  section  and  the 
divergent  section  where  the  peak  heat  flux  and  a stress  concentration  were 
combined.  A dark  ring  can  be  seen  in  Figure  15  at  this  point. 


Figure  15  Arc  heater  nozzle  after  25  atm  hydrogen  tests 


4. 1. 4 Electrode  Erosion  and  Test  Stream  Contamination 

The  electrode  weight  losses  were  divided  by  the  total  hydrogen  mass  flow 
between  weighings  to  determine  the  test  stream  contamination.  Several  tests 
were  made  at  the  240  A level  at  pressures  from  10  to  25  atm,  and  contamina- 
tion wis  determined  for  a single  current  level.  However,  many  of  the  tests 
were  made  at  three  arc  current  levels,  and  thus  the  contamination  was  represen- 
tative of  that  occurring  at  a median  arc  current.  Figure  16  showb  the  test 
stream  contamination  on  hydrogen.  At  240  A the  contamination  level  is  much 
lower  at  low  pressure  than  the  higher  current.  The  trend  shows  considerably 
less  difference  at  the  higher  pressure  levels. 
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Figure  17  shows  the  post-test  entrance  region  of  two  typical  front 
electrodes  (cathodes)  for  pressures  of  15  and  25  atm  and  arc  currents  of 
250  - 400  A.  At  15  atm,  evidence  of  shore  arcing  and  entrance  radius  erosion 
is  visible,  whereas  at  25  atm-  the  primary  erosion  is  further  inside  the 
cathode  barrel.  The  electrode  used  in  the  50  atm  test  did  not  experience 
the  entrance  damage.  Arc  tracks  were  visible  from  the  entrance  to  5.5  in. 
inside  the  barrel  with  only  minor  erosion. 


Pressure:  15  atm  OFHC  copper  Pressure:  25  atm 
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Tests  were  made  at  25  atm  and  240  - 400  A using  an  80%  silver,  20% 
copper  cathode.  Fig., re  17  compares  the  er.trance  reigion  of  that  electrode 
to  the  previously  tested  OFHC  copper.  The  erosion  was  more  severe,  but  the 
test  stream  contamination  (defined  as  the  electrode  mass  loss  as  a percent 
of  the  total  gas  flow  rate)  was  reduced  from  0*28%  to  0.22%.  The  reason  for 
the  reduced  contamination  was  increased  melt  solidification.  The  nozzle  and 
downstream  cathode  regions  were  heavily  coated  with  silver-copper  melt. 

Figure  18  shows  the  rear  electrode  Canode)  erosion  characteristics.  This 
anode  was  used  for  tests  1 through  34  at  10  to  15  atm  pressure  and  240  to  400  A. 
The  anode  used  for  tests  35  through  74  looked  very  similar  with  less  entrance 
corner  damage.  It  was  subjected  to  pressures  from  10  to  50  atm  and  arc 
currents  from  240  to  400  A.  The  anode  average  test  stream  contamination  for 
the  complete  test  matrix  (tests  36  to  74)  was  0.12%.  The  anode  exhibited 
erosion  and  melt  flow  over  the  first  1.38  in.  inside  the  barrel;  beyond  that 
Doint  there  were  circular  surface  patterns  to  3.00  in.  inside  the  barrel 
with  essentially  no  erosion. 


Post  tests:  1-34  Current:  240-400  A 

Pressure:  10-1 5 atm  OFHC  copper 

Figure  18  N-250-5.0.162  rear  electrode  (anode)  after  operation  on  hydrogen 


4.1.5  Hydrogen  Arc  Length 

Determination  of  the  hydrogen  arc  length  received  considerable  attention 
in  order  to  derive  a correlation  equation  for  the  HYARC  scaling  program. 
Careful  measurements  of  the  arc  track  locations  were  made  throughout  the 
testing.  However,  for  economical  reasons  it  was  not  feasible  to  change 
electrodes  for  each  new  condition  nor  to  test  each  arc  current  separately. 
Thus,  arc  tracks  accumulated  over  a wide  band  in  each  electrode  making  it 
difficult  to  distinguish  the  new  tracks  from  the  old,  especially  midway  down 
the  0.375  in.  diam  cathode. 
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Figure  19  presents  the  arc  length  data  as  a function  of  arc  current  and 
pressure.  The  lengths  shown  are  the  combined  depths  in  the  anode  and  cathode 
plus  the  0.25  in.  gap  between  electrodes.  The  240  A points  shown  are  primarily 
from  single  current  tests  as  are  the  310  A points.  The  400  A points 
represent  estimated  lengths  based  on  the  brightest  tracks  from  multi-current 
tests.  The  HYARC  lengths  shown  are  discussed  in  Section  5. 


Figure  19  Hydrogen  arc  lengths 


4.  2 HELIUM  ARC  CHARACTERISTICS 


The  arc  heater  was  operated  on  helium  at  nominal  arc  currents  of  240, 

320  and  400  A at  nominal  pressures  of  25,  50,  75  and  100  atm  with  a few  inter- 
mediate levels.  The  operating  characteristics  were  similar  to  those  on  air 
with  better  stability,  significantly  less  electrode  erosion  and  higher  gas 
enthalpies. 

Table  A-2  in  the  appendix  summarizes  the  complete  helium  test  matrix 
Including  duplicate  points.  The  tabulated  values  were  extracted  from  the 
printed  data  sheets.  Each  point  represents  a stable  test  condition  where 
sufficient  time  was  allowed  for  equilibration.  No  nozzle  throat- area 
reduction  problems  were  encountered  on  helium.  The  throat  diameter  remained 
constant  at  0.162  in. 

All  helium  tests  were  made  using  one  pair  of  electrodes  and  one  nozzle. 
The  higher  pressure  tests  utilized  four  0.062  in.  diam  gas  injectors  (tests 
103  through  106).  The  latter  test  (106)  was  a duplicate  test  at  25  atm.  All 
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Other  tests  utilized  four  0.052  injectors.  The  larger  injectors  reduced  the 
supply  pressure  requirements  to  a readily  available  level  (2600  psi).  The 
injection  pressure  ratio  for  the  small  injectors  (0.052  in.)  ranged  from  1.18 
to  1.37;  the  lower  value  was  25  atm  and  400  A and  the  higher  value  was  75  atm 
and  240  A.  In  general,  the  higher  the  gas  enthalpy  at  a fixed  arc  pressure, 
the  lower  the  injection  nressure  ratio.  The  injection-pressure  ratio  with  the 
larger  injectors  (0.062)  ranged  from  1.15  to  1.28.  For  a given  test  condition, 
the  ratio  was  lowered  anproximately  6%  using  these  injectors. 

In  order  to  complete  test  105  at  100  atm  and  three  arc  current  levels 
each  for  30  s,  it  was  necessary  to  manifold  20  helium  Bottles.  The  Initial 
attempt  at  100  atm  Ctest  101)  fell  short  C90  atm)  because  of  supply  depletion. 
The  second  attempt  Ctest  103)  using  the  larger  Injectors  reached  and  held 
the  95  atm  level  with  12  supply  bottles.  The  third  attempt  (test  105)  using 
20  bottles  was  successful  ClOO  atm). 

4. 2. 1 Helium  Enthalpy 

The  energy  balance  enthalpies  achieved  on  helium  are  shown  in  Figure  20 
as  a function  of  arc  pressure  and  current.  The  helium  enthalpies  were  a 
factor  of  four  lower  than  the  hydrogen  enthalpies.  The  helium  flow  rates  were 
a factor  of  two  higher  than  those  on  hydrogen,  and  the  helium  arc  voltages 
were  more  than  a factor  of  two  lower  than  the  hydrogen  arc  voltages.  Thus, 
for  similar  thermal  efficient  tes,  the  factor  of  four  difference  in  the 
enthalpies  was  reasonable. 

The  peak  helium  enthalpy  was  9700  Btu/lb  at  25  atm  and  the  minimum  was 
4700  Btu/lb  at  100  atm.  At  the  400  A level,  the  factor-of-four  increase  in 
arc  pressure  caused  only  a 24%  drop  in  the  energy  balance  enthalpy;  at  the 
200  A level,  the  drop  was  more  severe,  38%. 
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The  predicted  enthalpies  using  the  HEARC  computer  code  are  also  shown  in 
Figure  20.  This  code  is  discussed  in  Section  5. 


Agreement  between  energy  balance  enthalpy  and  sonic  flow  enthalpy  was 
better  on  helium  than  on  hydrogen  primarily  because  of  a constant  throat 
area.  General  agreement  of  + 5%  was  found  using  the  sonic  nozzle  mass 
flowmeter  and  small  injectors.  When  using  the  larger  injectors  and  a sub- 
sonic mass  flow  calibration,  the  enthalpies  differed  by  8%  to  21%  with  the 
sonic  flow  values  consistently  higher.  The  sonic  flow  enthalpy  on  helium  is: 


h 


sf 


(7) 


A comparison  of  tests  95  to  106  and  100  to  104  from  Table  A-2,  indicated 
that  the  suspicious  parameter  was  the  mass  flow  rate.  Both  comparisons  were 
at  fixed  pressure  and  throat  area,  and  both  indicated  better  agreement  in  the 
enthalpies  when  the  sonic  mass  flowmeter  was  used.  However,  checks  of  the 
subsonic  Injector  flow  calibrations  did  not  reveal  any  significant  errors. 


4.2.2  Helium  Arc  Voltage 

The  helium  arc  voltages  were  lower  by  a factor  of  two  and  less  sensitive 
to  gas  flow  rate  than  the  hydrogen  arc  voltages.  Figure  21  shows  the  helium 
arc  voltages  as  a function  of  flow  rate  and  arc  current  for  the  full  test 
matrix  from  25  to  100  atm.  Above  50  atm  (0.04  Ib/s) , the  voltage  was  very 
insensitive  to  arc  current  and  depended  primarily  on  the  mass  flow  rate.  The 
helium  data  shown  were  correlated  using  the  following  equation: 


V = 


, - -.4  . 0.72 

1.5  X 10 


(8) 


The  worst  discrepancies  were  less  than  8%  over  the  full  range  of  the  data. 

Most  of  the  data  were  correlated  much  closer  using  Equation  (8). 

Figure  22  compares  the  correlated  arc  voltage  from  HEARC  with  the 
experimental  values.  The  arc  length  used  in  HEARC  was  based  on  post- test  arc 
track  locations  and  experimental  arc  voltages.  Considering  this  fact,  agree- 
ment is  quite  good.  These  comparisons  are  discussed  further  in  Section  5. 

The  pressure  dependence  was  approximately  the  0.736  power,  whereas  the  current 
dependence  at  constant  pressure  varied  from  the  0.20  power  at  100  atm  to  the 
0.36  pox7er  at  25  atm.  Note  that  this  arc  current  effect  was  largely  absorbed 
by  the  mass  flow  correlation,  which  is  typical  for  a Huels-type  arc  heater. 2 

The  arc  voltage  did  not  vary  with  time  at  a given  pressure  and  current. 
The  low  electrode  erosion  and  constant  throat  area  maintained  a very  steady 
arc  voltage.  The  voltage  oscillations  were  mild.  The  variation  in  the  point 
average  over  a 0.85  s period  varied  less  than  6%  for  a steady  current  and 
pressure. 
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Helium  flow  rate  (Ib/s) 


Figure  21  Helium  ore  voltage  characteristics 
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Figure  22  Comparison  of  measured  and  predicted  helium  arc  voltages 


4.2.3  Energy  Losses  on  Helium 

The  energy  losses  to  the  anode,  cathode-body,  and  nozzle  operating  on 
helium  are  given  in  Table  A-2.  Anode  losses  varied  from  18%  to  32%  of  the 
total  losses  over  the  full  matrix  of  the  test  conditions.  The  variation  of 
these  losses  with  current  and  pressure  is  shown  in  Figure  23.  The  magnitude 
of  the  anode  losses  consistently  increased  with  arc  current  and  arc  pressure, 
as  did  the  percent  of  the  total  losses.  The  anode  loss  dependence  on  current 
was  approximately  the  1.2  power,  and  the  dependence  on  pressure  was  nearly 
linear.  As  the  anode  termination  moved  rearward  with  increased  mass  flow 
(and  pressure),  it  encountered  a stronger  radial  magnetic  field  from  the 
spin  coil.  The  anode  area  exposed  to  radiation  thus  was  not  linear  with 
flow  rate;  hence,  the  total  anode  loss  became  less  pressure  dependent  as 
shown  in  Figure  23. 

The  front  electrode  (cathode)  losses  varied  from  56%  to  70%  of  the 
tbtairrosses  over  the  full  test  matrix.  The  magnitude  of  the  cathode  losses 
increased  with  arc  current  and  pressure  as  shown  in  Figure  24.  The  cathode 
percent  of  the  total  losses  decreased  with  increased  pressure  and  was  rela- 
tively insensitive  to  increased  pressure  and  to  arc  current,  ihe  predicted 
cathode  losses  from  HEARC  also  are  shown  in  Figure  24.  Good  agreement  was 
achieved  at  the  middle  pressure  range  as  discussed  in  Section  5.  In  contrast 
to  the  hydrogen  data,  these  data  do  not  indicate  a significant  body  loss. 

The  nozzle  energy  losses  on  helium  varied  from  12%  to  16%  of  the  total 
losses.  The  magnitude  of  the  nozzle  losses  increased  with  arc  current  and 
pressure,  whereas  their  percentage  of  the  total  loss  was  relatively  insensitive 
to  both  parameters. 
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Figure  23  Comparison  of  measured  and  predicted  anode  energy  losses  on  helium 


The  nozzle  energy  losses  were  correlated  as  a function  of  the  gas  flow 
rate  and  energy  balance  enthalpy  as  shown  in  Figure  25.  A wall  temperature 
of  600®F  was  » ‘ ' in  the  correlation.  No  attempt  was  made  to  vary  the  wall 

temperature  . heat  flux.  All  of  the  data  except  test  95  correlated  within 
+ 5%  using  t following  equation: 


n n/ 

= 0.04 


h ) 

w 


(9) 


$ 


Equation  (9)  closely  follows  typical  turbulent  convective  heat  transfer 
correlations.  The  weaker  dependence  on  mass  flow  rate  (0.7  power  rather 
than  the  usual  0.8  power)  was  somewhat  arbitrary.  Inclusion  of  the  data  from 
test  95  results  in  a closer  correlation  with  m^*®,  but  the  accuracy  is 
reduced  to  + 15%. 
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The  peak  nozzle  loss  occurred  at  100  atm  and  400  A.  That  less  was 
converted  tc  a peak  throat  heat  flux  of  6700  Btu/ft^s  using  a Bartz  ^ heat 
flux  distribution  integrated  over  the  nozzle  surface  area.  No  damage  to 
the  nozzle  was  encountered  in  the  helium  tests. 


25  50  75  100 

Arc  pressure  (atm) 


Figure  24  Comparison  of  measured  and  predicted  cathode 
energy  losses  on  helium 
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Helium  flow  rate,  m|^g  (Ib/sl 

Figure  25  Correlation  of  arc  heater  nozzle  heat  transfer  on  helium 


4.2.4  Electrode  Erosion  and  Test  Stream  Contamination 


At  the  beginning  of  the  helium  testing,  the  cathode  weight  was  510.4  g 
and  the  anode  weight  was  510.3  g.  After  15  tests  on  helium  at  pressures  to 
100  atm  and  currents  to  400  A,  the  cathode  weight  was  503.8  g,  a loss  of 
6.6  g,  and  the  anode  weight  was  510.2  g,  a loss  of  0.1  g.  The  total  helium 
used  in  these  tests  was  approximately  67.2  lb.  Thus  the  average  test  stream 
contamination  (electrode  mass  loss)  for  all  tests  was  0.02%  of  the  total  gas 
flow  rate.  Figure  26  shows  the  variation  of  the  helium  test  stream  contamina- 
tion with  arc  pressure.  The  higher  mass  flow  rate  at  higher  pressures  moves 
the  arc  termination  more  rapidly,  reducing  electrode  material  loss  and  melt 
flow. 


The  anode  surface  virtually  undamaged  in  the  helium  test  matrix. 
The  surface  condition  varied  slightly  with  the  arc  pressure  and  current. 
The  cathode  surface  was  roughened  but  did  not  lose  significant  material. 
The  entrance  corner  and  downstream  barrel  were  undamaged. 
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Figure  26  Helium  test  stream  contamination 


4*  2. 5 Helium  Arc  Length 

Figure  27  shows  the  arc  track  locations  in  the  anode  as  a function  of 
helium  arc  pressure.  Since  the  arc  current  was  varied  during  a given  test, 
its  positioning  effect  was  not  documented.  However,  the  surface  roughness  is 
normally  more  severe  at  high  arc  current,  which  was  found  to  be  nearest  the 
downstream  end  of  the  anode.  It  is  characteristic  of  Huels-type  arc  heaters 
that  arc  shunting  is  enhanced  by  both  higher  arc  current  and  higher  arc 
potential  gradients.^  Extension  phenomena  include  gas  forces  from  the  reverse 
flow  and  Lorentz  forces  from  the  centering  field  coil  rear  of  the  arc  termina- 
tion. The  data  indicate  that  the  extension  forces  predominate  to  60  atm  at 
which  point  the  shunting  and  restrike  phenomena  begin  to  reduce  the  arc  length 
in  the  anode  barrel. 

Typical  cathode  arc  tracks  are  shown  in  Figure  28.  The  maximum  arc 
extension  was  linear  with  pressure  (flow  rate)  ranging  from  1.50  to  4.25  in. 
do\^nstream  of  the  cathode  barrel  entrance.  The  cathode  was  5 in.  long  in 
all  these  tests.  The  minimum  restrike  position  is  also  shown  in  Figure  28. 

Arc  column  shunting  (restrike)  occurs  near  the  entrance  (0.75  to  1.5  in.). 
Except  for  the  100  atm  tests,  the  restrike  tracks  were  in  essentially  the 
same  position.  The  arc  attachment  region  was  divided  into  two  zones:  the 

rough  zone  probably  characterizes  the  400  A arc  erosion,  and  the  spiral  track 
zone  probably  represents  the  240  A arc  extension  tracks. 
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Figure  27  Helium  arc  anode  positions 

An  accurate  definition  of  the  total  arc  length  was  not  feasible,  but 
approximations  could  be  inferred  from  the  arc  track  locations.  Figure  29 
shows  the  range  of  possible  arc  lengths  from  the  arc  tracks.  The  maximum 
possible  length  is  the  sum  of  the  maximum  anode  plus  cathode  positions  plus 
the  electrode  gap  (0.25  in.).  The  minimum  possible  length  is  the  sum  of 
the  shortest  positions  and  the  gap.  The  possible  arc  length  range  is  quite 
wide  (+  2 in.).  The  most  probable  arc  length  combined  the  minimum  anode 
track  depth  with  the  maximum  cathode  track  depth  and  vice  versa.  That  narrowed 
the  range  as  shown  in  Figure  29-  The  average  arc  length  was  relatively 
constant  above  50  atm  with  the  possible  arc  length  span  increasing  with 
pressure. 

The  arc  length  required  for  an  arc  voltage  match  using  HEARC  is  also 
shown.  There  was  reasonable  agreement  between  the  most  probable  arc  lengths 
and  required  arc  lengths  at  the  higher  pressures.  The  inferred  arc  lengths 
do  not  include  any  arc  twisting  which  has  beer  reported  in  the  past.  Future 
studies  should  include  arc  visualization  for  a better  definition  of  the  arc 
characteristics. 
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Figure  29  Helium  arc  lengths 


4,3  AIR  ARC  CHARACTERISTICS 


The  arc  heater  was  operated  on  air  at  nominal  arc  currents  of  240,  320 
and  400  A at  nominal  arc  pressures  of  10,  15,  25,  50,  75  and  100  atm. 

Table  A-3  summarizes  the  complete  air  test  matrix,  including  duplicate 
polnts>  The  tabulated  values  were  extracted  from  the  prl*ited  data  sheets, 
and  each  point  represents  a stable  test  condition  where  sufficient  time  war, 
allowed  for  equilibration* 
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The  7 Q ballast  resistance  and  the  additional  line  inductance  were  used 
for  all  air  tests.  Early  tests  on  nitrogen  indicated  that  some  of  the  lower 
pressure  conditions  could  have  been  run  without  these  additional  circuit 
components,  but  they  did  add  to  the  stability  and  allowed  direct  stability 
comparisons  to  be  made  over  the  full  test  range  of  the  other  gases. 

Tests  75  through  78  were  made  using  four  0.032  in.  diam  gas  injectors. 
The  injection  pressure  ratio  varied  from  1.10  to  1.34.  Tests  79  through  84 
were  made  using  four  0.025  in.  diam  gas  injectors.  The  Injection  pressure 
ratio  for  these  tests  varied  from  1.79  to  2.13.  Some  performance  improvement 
was  seen  at  the  10  atm  level  using  the  smaller  injectors,  but  there  was  no 
significant  difference  at  the  25  atm  level.  Tests  92  through  94  and  107 
through  109  were  made  using  four  larger  injectors  (0.036  in.  diam)  to  reduce 
the  gas  system  pressure  level.  The  injection  pressure  ratio  varied  from  1.24 
to  1.  40. 

All  of  the  air  tests  in  Table  A-3  were  made  using  the  5 in.  long  cathode 
and  the  0.162  in.  diam  nozzle  throat.  The  eroded  electrode  material  coated 
the  nozzle  throat  resulting  in  reduced  gas  flow  rates  for  a fixed  pressure 
and  discrepancies  between  energy  balance  and  sonic  flow  enthalpies.  Normally 
the  post-test  throat  diameter  was  approximately  0.158  in.,  but  after /tests 
76  and  78,  it  was  less  than  0.140  in.  The  throat  was  cleaned  after  each  test. 

4.3.1  Air  Enthalpy 

The  energy  balance  enthalpies  achieved  on  air  are  shown  in  Figure  30 
as  a function  of  arc  pressure  and  current.  The  air  enthalpies  were  a factor 
of  5 to  7 lower  than  those  on  hydrogen  and  41%  to  69%  lower  than  those  on 
helium  for  the  same  arc  current  and  pressure^ 
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The  peak  air  enthalpy  was  7200  Btu/lb  at  10  atm,  and  the  minimum  was 
3400  Btu/lb  at  an  order  of  magnitude  higher  pressure.  At  a given  arc  current, 
a factor  of  ten  increase  in  pressure  decreased  the  enthalpy  by  as  little  as 
33%  and  at  most  44%. 

TTie  predicted  enthalpy  values  from  AIRARC  I are  also  shown  in  Figure  30. 
These  predictions  were  based  on  arc  lengths  derived  from  the  experimental 
arc  tracks.  The  actual  air  flow  rates  were  used  instead  of  the  sonic  flow 
relation  to  avoid  throat-area  change  problems.  The  nozzle  energy  loss  corre- 
lation developed  in  Section  4.3.3  was  added  to  the  calculated  arc  region 
losses  to  determine  convective  losses  downstream  of  the  arc.  Agreement  is 
good  at  the  higher  pressures,  with  the  worst  agreement  occurring  at  the 
high  arc  current  and  low  pressure.  The  factors  involved  in  this  comparison 
are  further  discussed  in  this  section  and  Section  5. 

Agreement  between  energy  balance  enthalpy  and  sonic  flow  enthalpy  in 
certain  tests  below  25  atm  was  reasonably  good.  These  tests  did  not  result 
in  significant  throat-area  reduction.  In  a majority  of  the  tests  on  air, 
particularly  above  25  atm,  the  sonic  flow  enthalpy  was  consistently  higher  by 
25%  or  more.  A coated  nozzle  throat  diameter  of  0. 155  in.  would  have  caused 
this  enthalpy  difference.  Normally  the  post-test  throat  diameter  was  0.158  in. 
after  the  high-pressure  tests. 


4.3.2  Air  Arc  Voltages 


The  air  arc  voltages  varied  from  400  to  3090  V as  the  pressure  was 
increased  from  10  to  100  atm.  Figure  31  shows  the  experimental  values  as  a 
function  of  arc  current  and  pressure.  As  expected,  the  arc  voltages  were 
higher  at  the  lower  nrc  currents.  The  sensitivity  to  arc  current  decreased 
from  10  to  25  atm,  then  increased  again  as  the  pressure  was  increased  to 
100  atm. 

The  arc  voltage  and  current  were  very  stable  over  the  full  range  of 
test  conditions.  The  10-scan  average  (0.85  s elapsed  time)  voltage  at  240  A 
and  100  atm  varied  less  than  + 1%.  Above  15  atm,  the  arc  voltage  repeatability 
was  remarkable  compared  with  other  Huels-type  heaters.  At  50,  75  and  100  atm, 
the  repeat  test  voltages  were  identical  within  the  accuracy  of  the  measurements. 

The  predicted  arc  voltages  shown  in* Figure  31  were  based  on  arc  lengths 
inferred  from  the  arc  tracks  post-test.  These  lengths  are  discussed  in 
Section  4.3.5.  The  predicted  arc  voltages  could  possibly  have  been  matched 
closer  by  using  a technique  similar  to  that  used  in  HYARC,  but  the  agreement 
shown  in  Figure  31  was  judged  sufficiently  accurate.  This  technique  was 
also  a better  test  of  the  scaling  accuracy  of  AIRARC. 
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Arc  pressure  (atm) 


Figure  31  Air  am  voltages 

4.3.3  Energy  Losses  on  Air 

The  energy  losses  to  the  anode,  cathode-body,  and  nozzle  are  given  in 
Table  A- 3 for  the  complete  air  test  matrix.  The  anode  losses  varied  from 
14%  to  26%  of  the  total  losses  and  consistently  increased  with  arc  current 
and  pressure.  The  cathode-body  losses  varied  from  58%  to  74%,  and  the  nozzle 
losses  varied  from  11%  to  21%  of  the  total  energy  losses. 
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The  predictions  shown  in  Figure  32  were  made  using  AIRARC  I which 
combines  volumetric  radiation  (optically  thin)  and  a standard  convective 
correlation.  Although  the  absolute  magnitudes  of  the  predictions  are  not 
significantly  in  error,  particularly  in  the  middle  pressure  range,  the 
functional  dependence  on  pressure  appears  high  (0.92  power).  The  data  were 
quite  repeatable  and  consistent.  A much  closer  correlation  of  the  data  was 
achieved  using  AIRARC  II  as  discussed  in  Section  5. 

The  air  nozzle  energy  losses  (Figure  33)  were  correlated  using  the  same 
equation  as  for  the  hydrogen  nozzle  losses: 

Q„  = 0.0469  a (h^  - h ) . (11) 

a b w 

Equation  (11)  has  the  same  general  turbulent  convection  form  as  the  helium 
nozzle  loss  equation  and  standard  pipe  correlations.  The  coefficient  and 
mass  flow  rate  exponent  are  slightly  different. 

The  peak  nozzle  loss  occurred  at  100  atm  and  400  A.  That  loss  was 
converted  to  a peak  throat  heat  flux  of  5900  Btu/ft^s  using  a Bartz^  heat 
flux  distribution  integrated  over  the  nozzle  surface  area.  No  damage  was 
encountered  to  the  nozzle  during  the  air  tests. 


Air  flow  rate,  m3  (Ib/s) 

Figure  33  Correlation  of  the  air  nozzle  energy  toss 
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4. 3*4  Electrode  Erosion  and  Test-Stream  Contamination 

The  test-stream  contamination  on  air  was  determined  for  a wide  range  of 
conditions  using  pre-  and  post- test  electrode  weights.  Figure  34  depicts 
the  results  as  a function  of  arc  pressure.  The  level  was  generally  between 
0.16%  and  0.27%.  This  was  of  the  same  order  as  that  on  hydrogen  at  25  atm 
but  considerably  lower  than  that  on  hydrogen  at  10  atm  and  much  higher  at 
50  atm.  The  air  contamination  was  an  order  of  magnitude  higher  than  that 
on  helium. 

All  reported  tests  at  pressures  of  25  atm  and  below  were  made  utdng  an 
80%  silver,  20%  copper  anode.  From  Figure  34  there  appears  to  be  a significant 
reduction  in  test-stream  contamination  at  the  lower  pressures  using  the  silver 
alloy.  This  agrees  with  past  experience  with  the  Huels-type  arc  heater.  9 


Arc  pressure  (atm) 

Figure  34  Test  stream  contamination  on  air 
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4, 3. 5 Air  Arc  Length 


The  arc  length  in  the  N-250  operating  on  air  was  required  as  an  input 
to  the  AIRARC  scaling  program.  Measurements  of  the  arc  track  locations  were 
made  when  discernable.  The  heavy  oxide  buildup  and  melt  flow  caused  consider- 
able ambiguity  in  these  measurements.  Figure  35  shows  the  best  arc  length 
estimates  made  from  the  arc  track  locations.  There  was  a general  trend  for 
the  arc  to  lengthen  as  the  gas  flow  rate  (pressure)  was  Increased.  The  most 
probable  arc  length  range  shown  pairs  the  anode  res trike  position  with  the 
cathode  extension  position  and  vice  versa. 

The  arc  lengths  required  by  AIRARC  for  a good  match  of  the  experimental 
voltages  are  also  shown  in  Figure  35.  For  pressures  from  15  to  100  atm,  these 
lengths  fall  within  the  span  of  the  most  probable  arc  lengths. 

The  6 in.  arc  length  selected  prior  to  testing  was  based  on  previous 
scaling  studies  discussed  in  Section  5.  As  can  be  seen  in  Figure  35,  the 
6 in.  length  was  a fair  approximation  over  the  higher  pressure  range  (>  25  atm). 

It  should  be  reiterated  that  the  use  of  the  arc  tracks  for  an  inferred 
arc  length  does  not  Include  consideration  of  any  twisting  which  has  been 
reported  in  the  past. ^ ^ 


Figure  35  Arc  length  operating  on  air 
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4.4  COMPARISON  OF  HYDROGEN,  HELIUM  AND  AIR  ARC  CHARACTERISTICS 

The  performance  characteristics  of  the  N-250  arc  heater  on  hydrogen, 
helium  or  air  have  been  described  in  the  preceding  sections  and  some  comparisons 
have  been  made.  This  section  provides  a convenient,  compact  comparison  of  the 
major  performance  parameters. 

4. 4. i Enthalpy  Comparisons 

Table  3 summarizes  the  energy  balance  enthalpies  measured  on  the  three 
gases  tested  over  the  complete  test  matrix.  Direct  comparisons  can  be  made 
at  25  and  50  atm  among  all  three  gases.  The  hydrogen  enthalpies  were 
approximately  four  times  those  on  helium  and  five  to  seven  times  those  on 
air  for  the  same  arc  current  and  pressure.  At  the  high  pressures,  the  helium 
enthalpies  were  41%  to  69%  higher  than  air. 


TABLE  3 COMPARATIVE  N-250-5  0.162  ARC  HEATER  PERFORMANCE  ON  Hj,  He  OR  AIR 


Arc 

current 

(A) 

Arc  pressure  (atm) 

Test 

gas 

10 

15 

25 

50 

75 

100 

Average  energy  balance  enthalpy  (Btu/lb) 

H2 

30,200 

33,300 

30,900 

- 

- 

- 

240 

He 

- 

- 

7,300 

5,800 

5,100 

4,800 

Air 

5,600 

4,900 

4,300 

3,700 

3,500 

3,400 

1^2 

3;'  800 

35,700 

35,000 

31,200 

- 

~ 

320 

He 

- 

- 

8,300 

7,200 

6,200 

6,100 

Air 

5,900 

5,700 

5,100 

4,400 

4,100 

3,900 

H2 

32,900 

38,900 

37,800 

- 

- 

- 

400 

He 

- 

- 

9,300 

8,300 

7,600 

7,400 

Air 

7,000 

6,200 

5,700 

4,900 

4,500 

4,400 

The  performance  advantage  of  hydrogen  is  shown  in  Figure  36  where  the 
significantly  higher  enthalpies  from  Table  3 are  compared  graphically.  Also 
shown  are  the  performance  predictions  of  the  three  versions  of  the  minimum 
energy  addition  scaling  program. 
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Figure  36  Comparison  of  enthalpy-pressure  performance  on  H2,  He  or  air 
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4. 4. 2 Arc  Voltage  Comparisons 

Figure  37  compares  the  arc  voltages  on  hydrogen,  helium  and  air  for  the 
same  and  similar  arc  currents  and  pressures.  The  hydrogen  arc  voltages  were 
more  than  a factor  of  two  higher  than  those  on  air,  and  the  air  voltages 
were  approximately  10%  higher  than  those  on  helium. 

The  voltage  and  current  oscillations  on  hydrogen  were  more  severe  than 
those  on  air,  with  the  helium  arc  characteristics  being  the  most  stable  of 
the  three.  Noise  frequencies  of  150  kHz  with  a modulation  of  2.5  kHz  were 
identified  on  hydrogen.  An  isolation  amplifier  was  necessary  to  obtain  an 
accurate  hydrogen  arc  voltage. 

4. 4. 3 Test-Stream  Contamination  Comparison 

The  electrode  erosion  and  resultant  test-stream  contamination  varied 
significantly  with  the  test  gas  and,  on  hydrogen,  with  the  test  pressure. 
Figure  38  compares  the  total  test  stream  contamination  (anode  plus  cathode 
mass  loss)  as  a percent  of  the  gas  mass  flow  rate  for  all  three  gases. 

At  low  pressure,  the  arc  was  relatively  short  on  hydrogen  and  was  not 
evenly  distributed  on  the  electrode  surface.  The  result  was  severe  erosion 
and  high  contamination.  The  low-pre. .are  air  tests  were  made  using  a silver- 
copper  anode  which  appears  to  have  helped  reduce  the  total  contamination. 
Contamination  on  air  was  more  severe  than  on  hydrogen  or  helium  at  the  higher 
pressures.  Anode  erosion  at  75  and  ICO  atm  on  air  was  quite  severe,  whereas 
the  erosion  was  negligible  on  helium  for  the  same  conditions. 
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5.  ARC  HEATER  SCALING 


Scaling  of  Huels-type  arc  heaters  received  considerable  attention  a 
decade  ago^^lO,  and  simple  rules  were  derived  from  data  on  small  arc  heaters. 
The  use  of  these  simple  scaling  rules  led  to  sizable  design  errors  for  larger 
scale  heaters. ^ 1 The  power  capability  of  the  heaters  was  under-estimated  by 
a factor  of  four,  and  the  scaled  arc  voltages  were  high,  leading  to  predicted 
enthalpies  considerably  above  those  actually  achieved.^ 

1 7 

Later  work  eveloped  an  analytical  model  for  high  pressure,  high  flow 
rate  arc  heaters  which  included  radiation  from  the  arc  column  and  turbulent 
convective  heating  of  the  gas.  This  model,  for  fixed-length  arc  heaters,  was 
computer  coded  and  shown  to  yield  reasonable  agreement  with  limited  experi- 
mental data.  A study  by  ^^)RL2  resulted  in  the  addition  of  the  sonic  flow 
relation  and  wall  convection  losses,  and  in  improvements  of  the  transport 
and  thermodynamic  properties  used  in  the  above  model. 

Comparisons  with  data  from  fixed- length  arc  heaters  were  favorable  in 
the  moderate  pressure  regime  (25  to  125  atm).^^  Comparisons  with  Huels-cype 
arc  heater  (MDC-200)  data  were  surprisingly  favorable  in  the  higher  pressure 
regime^^  (50  to  250  atm),  but  perhaps  fortuitous  since  the  model  did  not 
consider  the  variable  arc  position. 

In  order  to  assess  the  accuracy  of  the  modified  code  for  scaling  Huels- 
type  arc  heaters,  comparisons  with  experimental  arc  voltages  from  four 
different  scale  heaters  were  made. ^ ^ These  heaters  included  the  N-250 
(N  = 0.5)^^  used  in  this  study,  the  Linde  N-1000  (N  = 1)]^,  the  MDC-200 
(N  = 2)15  the  Linde  N-4000  (N  = 2)^^,  and  the  USAF-AFFDLRENT  heater  (N  = 4) . 
The  nozzle  throat  diameters,  throat  air  flow  rates  and  aic  currents  were 
normalized  using  the  arc  heater  scale  to  allow  direct  comparisons.  The  arc 
lengths  were  scaled  linearly  as  L/N  = 12  in.  Figures  39  and  40  illustrate 
the  results  of  these  comparisons  for  45  and  11  atm.  At  45  atm  pressure,  the 
scaled  voltages  range  from  76%  to  103%  of  the  measured  values.  At  11  atm, 
they  range  ^rom  84%  to  104%.  These  data  cover  a wide  range  of  variables  and 
arc  heater  scale,  thus  indicating  that  the  arc  length  can  be  scald  linearly. 

The  effect  of  the  individual  operating  variables  is  shown  in  Figure  41 
for  an  N = 1 heater 1^  where  the  arc  voltage  comparison  is  presented  as  a 
function  of  scaled  arc  current,  mass  flow  rate,  and  nozzle  throat  diameter. 
Again  the  arc  length  was  scaled  linearly.  The  strongest  effect  is  due  to 
mass  flow  rate  where  a factor  of  eight  increase  causes  only  a 23%  deviation. 

Having  established  that  the  improved  computer  code  could  be  used  with 
confidence  for  scaling,  attention  was  next  devoted  to  further  improvements 
that  better  modeled  the  lluels-type  heater  and  accurately  matched  the  air 
data  accumulated  in  this  experimental  program.  This  version  of  the  code  was 
labeled  AIRARC. 

The  gas  transport  properties,  thermodynamic  properties  and  other  features 
of  the  code  were  then  modified  for  hydrogen  (HYARC)  and  helium  (HEARC) . 
Correlation  of  the  experimental  data  on  these  gases  using  the  revised  codes 
plus  the  confidence  in  scaling 'Uir  data  implies  that  reliable  estimates  can 
be  made  for  larger  scale  hydrogen  and  helium  arc  heaters.  Such  estimates 
have  been  made  and  are  included  in.  this  section. 


64 


Al  RARC  scaled/actual  arc  voltage  Al R ARC  scaled/actual  arc  voltagt; 


Figure  40  Comparison  of  scaled  and  actual  voltages  in  Huels-type  arc  air  heaters 
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Figure  41  Effect  of  scaling  parameters  on  AIR  ARC  predicted  voltages 


5.1  SCALING  CODE 

12 

The  arc  heater  scaling  code  Is  based  on  a model  for  high  pressure,  high 
flow  rate  operation.  Energy  input  is  by  ohmic  heating  of  the  arc  column  and 
losses  Include  radiation  and  convection.  The  bulk  gas  is  heated  by  turbulent 
convection  from  the  arc.  The  solution  of  the  model  is  based  on  the  principle 
of  energy  conservation,  minimum  energy  addition,  total  gas  energy  balance, 
and  optional  nozzle  throat  sonic  flow.  The  algorithm  is  essentially  the  same 
for  each  operating  gas,  the  difference  being  the  use  the  appropriate 
thermodynamic  and  transport  properties. 

The  input  parameters  required  for  the  algorithm  are: 

o Arc  current 
o Chamber  pressure 
o Arc  heater  geometry 

Constrictor  diameter 
Front  constrictor  length 
Nozzle  throat  diameter 
Arc  length. 

Based  on  the  above  information,  the  algorithm  predicts  the  following: 


0 

Enthalpy 

o 

Convective  heat  flux 

0 

Voltage 

n 

Constrictor  heat  load 

0 

Mass  flow  rate 

o 

Rear  electrode  heat  load 

0 

Efficiency 

0 

Electrical  gradient 

0 

Electrical  power  input 

0 

Power  in  the  gas 

0 

Stagnation  temperature 

0 

Arc  temperature 

0 

Radiative  heat  flux 

0 

Arc  diameter 
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5.1.1  Physical  Model 


The  physical  model  for  high-pressure  arc  heater  operation  is  based  on 
the  analysis  of  Richter.  ^ ^ The  justification  and  rationalization  of  this 
approach  is  presented  in  Reference  1 2 . The  incorporation  of  total  gas 
energy  balance  for  determining  bulk  gas  temperature  and  convective  losses  to 
the  constrictor  wall,  and  of  nozzle  sonic  flow  relation  for  determining  mass 
flow  rate  are  from  Painter. 1 ^ 

The  physical  model  characterizing  a high  pressure,  high  flow  rate  arc 
heater  is  illustrated  in  Figure  42.  The  arc  column  is  considered  as  a central 
conductive  body  (wire  like)  at  a uniformly  high  temperature  and  is  the  only 
source  of  ohmic  heat.  The  bulk  gas  in  the  annulus  surrounding  the  arc  is 
heated  by  convective  turbulent  transfer  from  the  arc.  In  addition  to  heating 
the  gas,  the  arc  column  also  loses  energy  directly  to  the  constrictor  wall  by 
volumetric  radiation.  The  arc  and  bulk  gases  are  considered  to  be  optically 
thin.  The  bulk  gas  loses  energy  tc  the  constrictor  wall  by  convection. 


Figure  42  High  pressure  arc  heater  mode! 
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The  algorithm  has  been  adapte*  to  characterize  Huels-type  configurations. 
Figure  43,  by  allowing  for  the  foil  ^wing  types  of  energy  losses  which  are 
incorporated  into  the  detailed  energy  balance  equatio.;s. 

1.  Kear  Electrode  - This  loss  is  by  radiation  frr-j  the  arc  only.  Since 
a portion  of  the  room- temperature  inlet  gas  has  reverse  flow  along 
the  electrode  wall,  hence  shielding  it  from  the  hot  core  gas,  there 
is  no  convection.  The  bulk  gas  is  still  heated  convectively  from 
the  arc,  but  i ^re  is  no  convective  loss  to  the  wall  of  the  rear 
electrode. 

2.  Front  Electrode  - These  losses  are  by  radiation  from  the  arc  only, 
over  the  distance  that  the  arc  extends  into  the  constrictor. 
Convection  from  the  hot  bulk  gas  is  allowed  over  the  entire 
constructor  length. 

3.  Nozzle  “ Losses  to  the  nozzle  are  includeci  by  using  an  empirical 
convection  relation  based  on  experimental  t:ata. 


Rear  electrode  (anode)  Front  electrode  (cathode)  Nozzle 


This  arrangement  recognizes  some  of  the  realities  of  the  Huels-type 
configuration,  namely  that  convective  losses  to  the  wall  of  the  front  electrode 
do  occur  from  the  hot  gas  downstream  of  the  arc  termination,  that  the  rear 
electrode  has  no  convective  losses  resulting  from  the  reverse  flow  of  the 
cool  inlet  gas,  and  that  there  are  losses  to  the  nozzle. 

Tlie  front  and  rear  arc  lengths  and  the  front  constrictor  lengths  were 
input  pa* imeters  the  algorithm,  so  that  other  arc  heater  configurations 
could  also  be  simulated. 

The  following  equations  show  the  development  of  the  model.  The  first 
part  of  the  analysis  is  presented  for  completeness  only;  it  may  also  be  found 
in  Reference  12. 
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The  power  balance  on  the  arc  column  per  unit  length  is: 

P = P + P , 
ire 

The  electrical  power  input  to  the  arc  per  unit  length  is: 


P.  = . (13) 

71  ^ 2 D ^ 
o -r  D a 

4 a 

The  power  radiated  per  unit  length  of  the  arc  column  (to  the  constrictor 
walls)  is: 

P =U7D^  = YD^.  (14) 

r 4 a a 

The  power  convected  from  the  arc  to  the  bulk  gas  per  unit  length  is: 

P = h(T  - T ) TT  D . (15) 

c a g a 

The  heat  transfer  coefficient  from  the  arc  to  the  gas  is: 


ag  D-D 
c a 


The  Nusselt  number  is: 


Nu  = 0.0265  (Re)^’^  (Pr)°‘^  . 
The  Reynolds  number  is: 


7T  U TD  + D ) • 

c a 

Summarizing,  the  resultant  expression  for  P is; 

D (D  + D )^’^ 

p - 7 c a .. 

r - ^ 2 9 ’ V. 

^ D ^ - D “ 

c a 

where  Z is  a parameter  independent  of  arc  diameter  and  is  defined  as: 


Z = 0.0265 


, . \0.8 

4 m \ 


(Pr)®'^  X 71  (T  - T ) . 

a g' 


The  arc  diameter  is  obtained  by  solving  the  power  balance  equation. 
Equation  (12),  as  rewritten  in  terms  of  the  arc  diameter. 
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» D (D  + D ) 

; % 

D-D 
c a 


0.2 


= 0 


(21) 


Equation  (21)  Is  solved  for  given  values  of  arc  temperature,  bulk  gas 
temperature,  mass  flow  rate,  arc  current,  pressure,  and  constrictor  diameter. 
The  algorithm  solves  this  equation  for  arc  diameter  using  a Newton-Raphson 
iteration  technique  starting  with  an  initial  approximation  for  D^.  This 
solution  is  performed  in  subroutine  ARCDIA. 

The  arc  temperature  and  hencA  electrical  power  input.  Equation  (13) , 
is  determined  by  finding  the  arc  temperature  for  which  the  power  input  is  a 
minimum.  This  corresponds  to  the  physical  principle  that  processes  in  nature 
will  always  occur  for  the  least  possible  expenditure  of  energy.  It  is  from 
this  principle  that  the  algorithm  derives  the  title  of  minimum  energy  addition. 
Figure  44  illustrates  the  variation  of  power  input  with  arc  temperature.  A 
minimum  power  is  clearly  evident.  The  algorithm  determines  this  minimum 
power  by  stepping  the  arc  temperature  in  a direction  that  reduces  the  power 
input.  When  the  minimum  is  passed,  the  temperature  step  size  is  reduced  and 
proceeds  in  the  opposite  direction.  The  searching  process  continues  until 
a small  pre-set  temperature  step  is  reached.  The  determination  of  the  arc 
temperature  for  which  the  power  input  is  a minimum  is  performed  in  subroutine 
ARCTEM. 


The  bulk  gas  temperature  is  determined  from  a total  gas  power  balance 
for  the  arc  column  and  the  bulk  gas  annulus.  The  gas  temperature  is 
iterated  until  the  total  gas  power  based  on  arc  and  gas  annulus  mass  flow 
rates  and  enthalpies  agree  with  the  power  based  on  power  input  less  losses. 

The  total  gas  power  based  on  input  less  losses  is  obtained  from: 


P = + m H.  -(q  TT  D L + q tt  D L + Q J , 

i i ^^r  c a ^c  c c * 


(22) 


where  P = total  power  in  gas  based  on  input  less  losses. 


g 


Bal 


The  radiative  heat  flux  is  obtained  from  the  arc  column  volumetric 
radiation  loss: 


q 


r 


(23) 
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The  convective  heat  flux  from  the  bulk  gas  to  constrictor  wall  is 
obtained  as  follows: 


q ^ h (T  - T ) , 
gw  g w ' 


(24) 


where  T = 590  K and  h is  the  turbulent  film  heat  transfer  coefficient: 
w gw 


h 


gw 


Nu  A 

D-D 
c a 


(25) 


The  total  power  in  the  gas,  based  on  arc  column  and  gas  annulus  flow 
rates  and  enthalpies,  is  obtained  by  assuming  a uniform  velocity  for  the  arc 
and  the  gas.  This  assumption  is  in  contrast  to  that  made  previously  for 
computing  the  turbulent  heating  of  the  gas  by  the  arc  where  the  velocity  was 
implicitly  assumed  to  be  zero.  The  assumption  of  non-zero  arc  column  velocity 
is  required  to  obtain  an  arc  column  mass  flow  rate  which,  however  small,  is 
required  to  obtain  the  power  in  the  arc  column.  This  is  not  insignificant 
because  of  the  high  arc  column  enthalpy.  The  velocity  is  obtained  from: 


p A + p A 
g g a a 

The  flow  rates  for  the  arc  column  and  gas  annulus  are  hence: 

m = p A V , and 
a a a 


m = p A V . 
g g g 


(26) 


(27) 

(28) 


The  total  power  in 


the  gas  based  on  enthalpy  and  mass  flow  is: 


P 

g 


Ien 


m + H m 
g a a 


(29) 


The  computation  of  these  two  expressions,  Equations  (22)  and  (29),  for 
total  gas  power  is  performed  in  subroutine  GASPOW,  The  bulk  gas  temperature 
is  iterated  by  the  method  of  false  position  until  total  gas  power  based  on 
Equations  (22)  and  (29)  agree  within  a specified  tolerance.  The  iteration 
process  is  performed  in  the  main  program. 
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The  mass  flow  rate  is  determined  by  assuming  the  flow  to  be  sonic  at  the 
nozzle  throat: 


The  mass  flow  rate  is  compared  with  that  predicted  by  the  above,  and  iterated 
until  the  two  values  agree  within  a specified  tolerance.  These  computations 
are  performea  in  the  main  program. 

The  stagnation  enthalpy  is  the  total  poTOr  in  the  gas  divided  by  the 
mass  flow  rate: 


H 

o 


(31) 


The  stagnation  temperature,  T , is  the  temperature  that  corresponds  to  the 
stagnation  enthalpy,  H , at  pressure  P.  This  determination  is  done  in  sub- 
routine TEMP.  ° 


The  voltage  gradient  is  obtained  from: 


G 


a A 
a a 


The  arc  voltage  is  thus: 
V = G La  . 


(32) 


(33) 


where  the  arc  length  L is  an  input  parameter  to  the  algorithm. 

a 

Efficiency  is  defined  as  the  power  acdec  to  the  gas  divided  by  the  input 
electrical  power: 


P " m H 

n = -^-5 * 100  . (34) 

i 

The  rear  electrode  heat  loss  is  assumed  to  be  due  only  to  radiation, 
since  for  a Huels-type  heater  the  reverse  flow  of  the  inlet  gas  at  room 
.temperature  would  preclude  convective  transfer.  Hence, 


q 7T  D L 
^r  c r 


where  is  the  length  of  the  arc  in  the  rear  electrode. 


(35) 
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The  front  electrode  Cconstrictor)  heat  loss  is  due  to  radiation  over  that 
portion  of  the  arc  extending  into  this  electrode,  and  convection  from  the  bulk 
gas  over  the  entire  electrode  length.  Hence: 


Q^=q  tiD  L-  + q ttD  L . (36) 

^r  c f ^c  c c 

The  total  arc  length  is: 

L^.lp  + Lj  • ™ 

Equations  (32)  through  (36)  are  evaluated  in  subroutine  OUTPUT. 

5.1.2  Operating  Gas  Characterization 

The  rainimuin-energy-addition,  arc  heater  characterization  code  requires 
information  on  the  gas  thermodynamic  and  transport  properties  as  a function 
of  temperature  and  pressure,  nozzles  losses,  and  front  and  rear  electrode 
arc  lengths.  This  section  outlines  the  source  of  such  data  and  describes 
the  correlations  utilized. 

The  helium,  hydrogen,  and  air  thermodynamic  and  transport  data  were 
obtained  from  the  sources  shown  in  Table  4 over  the  pressure  and  temperature 
ranges  indicated.  These  data  were  incorporated  into  the  algorithm  by  means 
of  table  interpolation  on  both  temperature  and  pressure  for  helium  and 
hydrogen  with  the  exception  of  the  electrical  conductivity  of  hydrogen  which 
was  curve  fit.  AIRaRC  used  curve  fit  routines  and  closed  form  equations. 
Transport  properties  not  available  above  30  atm  were  obtained  by  extrapolation 
using  the  pressure  dependent  relationship  of  the  data  below  30  atm. 

The  Prandtl  number  used  for  each  gas  is  shown  in  Table  5.  A single 
constant  value  was  chosen  because  its  variation  over  the  temperature  and 
pressure  ranges  of  interest  was  not  great,  plus  the  fact  that  it  was  raised 
to  the  0.3  poorer,  (Pr)^*^,  which  tended  to  further  minimize  any  variations. 

The  ratio  of  specific  heats,  used  in  the  sonic  flow  relation  is  shown 
in  Table  5.  A single  constant  value  was  used  for  helium,  and  a pressure- 
temperature  curve  fit  was  used  for  hydrogen.  An  empirical  sonic  flow  relation 
was  used  for  air  that  precluded  the  need  for  y. 

The  room  temperature  inlet  enthalpy  used  for  each  gas  is  also  shown  in 
Table  5. 

Volumetric  radiation  for  air  and  hydrogen  as  functions  of  temperature 
and  pressure  were  obtained  from  Yos.  por  helium,  radiation  was  determined 
functionally  as  that  resulting  from  free-bound  transitions.  Kogan points 
out  that  the  main  mechanism  of  radiative  heat  transfer  for  a plasma  with  low 
electron  temperature  (of  the  order  of  electron  volts)  is  the  process  of 
recombination  of  electrons  with  ions  as  opposed  to  the  case  of  high  temperature 
plasma  (T  ^ 100  eV)  where  electron  bremsstrahlung  is  thamain  contribution. 

The  functional  relationship  for  free-bound  radiation,^ ^ ^ ® is: 

u = N ^ . C38) 

/T  ® 
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TABLE  4 THERMODYNAMIC  AND  TRANSPORT  PROPERTY  SOURCES  AND  RANGES 
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Gas  Prandtl  no.  Inlet  enthalpy  Specific  heat  ratio 

Helium  0.672  670  Btu/lb  1.6667 

Hydrogen  0.650  1800  Btu/lb  Curve  fit 

Air  0.680  130  Btu/lb  Not  used 
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The  minimum-energy- add it ion  arc  heater  scaling  algorithm  has  been 
correlated  to  the  experimental  data  of  this  program.  The  parameters 
(observables)  were  arc  voltage,  power  input,  efficiency,  front  and  rear 
electrode  losses,  enthalpy,  and  front  and  rear  electrode  arc  tracks.  The 
correlation  goal  was  to  obtain  as  good  an  agreement  as  possible  between 
algorithm  predictions  and  all  the  above  cbservables  over  the  entire  pressure- 
current  test  matrix.  The  algorithm  parameters  "adjusted**  to  obtain  this 
correlation  were:  relative  radiation  level;  the  temperature  used  to  evaluate 

the  thermal  conductivity  in  the  expression  for  energy  transfer  between  arc 
column  and  bulk  gas,  and  likewise  in  the  expression  for  energy  transfer 
between  the  bulk  gas  and  constrictor  wall;  and  front  and  rear  arc  lengths. 

The  rationale  for  this  approach  was  that  by  fitting  the  algorithm  model 
to  experiment,  a degree  of  confidence  for  scaling  to  other  heater  designs  was 
added.  This  approach  also  tended  to  remove  possible  discrepancies  of  the 
theoretically  calculated  thermodynamic  and  transport  properties.  The  results 
of  this  correlation  process  for  hydrogen  and  helium  are  described  below. 

Table  6 shows  the  observables  used  in  adjusting  each  of  the  algorithm 
parameters.  The  scatter  in  the  observed  arc  track  data  was  generally  too 
great  to  assign  well  defined  values.  Hence,  arc  length  was  made  an  adjustable 
parameter.  If  arc  lengths  had  been  well  defined  (for  each  pressure  and 
current) , then  the  correlation  processes  would  have  been  easier  and  probably 
more  accurate.  For  example,  if  the  rear  arc  length  were  known  with  certainty, 
then  the  radiation  level  could  have  been  directly  related  to  rear  power  losses. 
As  it  was,  a trade-off  was  required  to  interrelate  rear  length  and  radiation 
with  the  losses. 


TABLE  6 CORRELATION  PARAMETERS 


Algorithm  parameters 

Observable  correlated 

Total  arc  length,  Lj 

Arc  voltage  (power  input),  arc  tracks 

Rear  electrode  arc 
Length,  Lr 

Rear  electrode  losses,  arc  tracks, 
Lr  = Ly  -Lp 

Front  electrode  arc 
length,  Lp 

Front  electrode  losses,  arc  tracks 
Lp  = L-p  - Lr 

Radiation  level 

Rear  electrode  losses,  arc  voltage, 
efficiency 

Temperature  for 
evaluating  X^g 

Efficiency  (power  in  gas  after  arc 
voltage  is  correlated),  enthalpy 

Temperature  for 
evaluating  Xg^y 

Front  electrode  loss 

The  predicted  value  for  each  observable  is  dependent  upon  more  than  one 
adjustable  algoritlim  parameter;  hence  several  **adjustment"  iterations  were 
required  to  determine  the  final  values.  These  final  values  represent  a 
subjective  best-fit  of  the  computed  results  to  the  experimental  observables. 
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The  arc  length  results  are  shown  in  Table  7 and  in  Figures  19,  29,  and 
35  for  comparison  with  observed  arc  track  data.  For  heaters  other  than  scale 
N = 0.5,  the  arc  lengths  were  scaled  from  L/N  values  equal  to  those  obtained 
in  this  work.  For  example,  a scale  N = 1 heater  front  and  rear  arc  lengths 
was  scaled  as  txd.ce  those  observed  for  this  work.  For  fixed  arc  length 
heaters,  the  physical  electrode  gap  should  be  input  to  the  algorithm. 


TABLE?  ARC  LENGTH  CORRELATIONS:  N-250-5/8-0.1 62 


Helium 

Lp  = 4.45  + 0.0106  (P(j  - 25)  • Lr  (in.) 
Lr=Po/(0.288Po  + 4.22)  (in.) 

Hydrogen 

/320\ 

Lj  = 1.3  + (4.51  + 0.086)  (Pq  - 15)  (—  1 (in.) 

/40o\°-22 

Lp  = (3.2 +0.063  (Pq -10)  ( — ) (i"-) 

Lr  = Lt  - Lf  (in.) 

Air 

La  =5.0 + 0.02  Pq  (in.) 

The  temperature  at  which  the  thermal  conductivity  was  evaluated  for 
computing  the  power  transfer  from  the  arc  column  to  the  bulk  gas  was  determined 
to  be: 


T = 0.8  (T 

a 

- T ) + T 
g g 

Hydrogen 

(39) 

T = 0.7  (T 

- T ) + T 

rr  rr 

Helium 

(40) 

a g'  g 

The  temperature  at  which  the  thermal  conductivity  was  evaluated  for 
computing  the  power  transfer  from  the  bulk  gas  to  the  constrictor  wall  was 
determined  to  be: 

T=0.5(T  -T)+T  Hydrogen  (41) 

gw  w 

T = T Helium  (42) 

g 

The  nozzle  energy  loss  expressions  used  in  the  algorithm  for  use  in 
calculating  an  overall  gas  power  balance  are  shoxm  in  Table  8.  These  were 
correlated  from  the  experimental  data  of  this  program,  Section  4.  Generally, 
the  nozzle  losses  were  5-6%  of  the  input  power  and  were  not  critical  for 
reasonable  algorithm  predictions.  When  scaling  to  different  size  heaters, 
the  applicability  of  this  expression  should  be  examined. 
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TABLE  8 NOZZLE  LOSS  CORRELATIONS  FOR  HYARC,  HEARC  AND  AIRARC:  N-250-5/8-0.162 


Helium 

Qn  = 0.04  (hj,  - 1325)  3tu/s 

Hydrogen 

Qn  = 0.0469  (mH)®  ®(h5  - 3600)  BUi.'s 

Air 

Qn  = 0.0469  (rha)®-®  (»»b  ~ 250)Biu/s 

Radiation  levels  were  correlated  to  the  power  loss  of  the  rear  electrode. 
Rear  electrode  losses  were  assumed  to  be  entirely  radiative  because  reverse 
flow  of  the  inlet  gas  at  room  temperature  precludes  any  convective  transfer 
to  the  electrode  walls.  The  resultant  constant  for  the  helium  radiation. 
Equation  (38),  is  C - 4.34  x 10”^^  for  radiation  in  W/cm^  and  electron  density 


r3 


For  hydrogen,  Yos^^  radiation  levels  were  increased  by  a factor  of 


in  cm' 

2.4.  Yos  radiation  was  computed  according  to: 


U„  = 17.1  X 10~^^  /T  , (43) 

Yos  ^ 

which  was  an  expression  for  free-free  radiation  as  corrected  for  free-bound 
and  modified  by  an  empirical  correction  factor  based  on  radiation  data  for 
shock  heated  argon.  As  discussed  by  Kogan,  ^ 7 f^^^e-bound  radiation  predominates 
for  arc  heater  conditions.  He  has  computed  for  hydrogen: 

= 5.32  X 10“^^  N^/ZT  . (44) 

22  2 7 

Comparing  Yos’s  value  with  Kogan's  for  hydrogen  arc  temperatures  of 
11,500  K,  we  obtain: 


”fb  3.11x10''  . 

_ — = =2.71,  (45) 

xos 

which  nearly  corresponds  to  the  factor  of  2.4  found  in  this  work. 

The  helium  radiation  is  8.1  times  that  for  hydrogen.  This  can  be  partly 
explained  by  the  higher  ionization  potential  of  helium,  and  possibly  by  non- 
unity  Gaunt  factors.  In  addition,  the  radiation  constant  (Equation  44)  varies 
as  the  square  of  the  effective  charge  on  the  nucleus,  Z^.  A value  of  unity 
was  used  for  hydrogen.  For  singly  ionized  helium,  the  effective  nuclear 
charge  is  1 ^ Z £ 2 since  its  remaining  electron  cannot  effectively  shield 
the  inner  core.  Thus,  up  to  a factor  of  4 could  be  accounted  for  in  this 
manner. 
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5. 1. 3 Discussion 


As  with  any  physical  modeling  process,  trade-offs  must  be  made  between 
solving  simpler  models  which  include  the  principal  physical  processes  (and 
hence  more  cost-effective)  and  solving  more  detailed  models  v/hich  may  include 
parameters  only  vaguely  known  and/or  so  involved  that  significant  phenomena 
such  as  radial  pressure  gradients  and  swirl  must  be  ignored  to  obtain  a 
solution. Thus  simpler  approaches,  particularly  when  correlated  with  data, 
can  often  be  adequate  for  engineering  estimation  purposes. 

The  physical  model  upon  which  this  algorithm  is  based  includes  the 
significant  principal  mechanisms  for  energy  input,  arc  column  radiation  loss, 
turbulent  convective  heating  of  the  gas,  bulk  gas  power  level  determination, 
bulk  gas  wall  losses,  and  mass  flow  rate  determination. 

This  approach,  however*  has  certain  inherent  limitations.  The  analysis 
is  a single-point  solution. ^ ^ It  evaluates  a stable  operating  condition  at 
? single  location  and  projects  this  solution  for  the  entire  arc  length.  In 
practice,  however,  high  pressure  constricted  arcs  are  not  fully  developed. 

The  energy  exchange  between  the  arc  and  bulk  gas  annulus  is  based  on  turbulent 
heat  transfer  correlations  that  inadequately  consider  many  of  the  details  of 
the  heating  process.  For  that  portion  of  the  arc  in  the  rear  electrode  where 
reverse  flow  occurs,  the  program  assumes  the  entire  flow  rate  for  calculating 
the  convective  gas  heating  rate;  in  practice,  only  a fraction  of  the  inlet 
gas  experiences  reverse  flox^  in  the  rear  electrode.  Another  undetermined 
parameter  for  the  rear  electrode  is  the  effective  annular  cross  sectional 
area  for  the  flow.  S^.nce  the  flow  is  reversed,  an  effective  diameter 
corresponding  to  the  front  constrictor  diameter  has  been  utilized.  If  in 
practice  this  diameter  was  smaller,  it  would  tend  to  offset  the  reduced  flow 
rate  in  the  Reynolds  number  determination. 

5.2  AIRARC  COMPARISONS  AND  SCALING 


Computations  using  AIRARC  were  made  for  the  full  range  of  the  experimental 
N-250  data  from  this  program.  The  initial  comparisons  indicated  that  the  nozzle 
throat  area  change  during  test  caused  the  sonic  flow  rates  to  be  low,  thus 
yielding  high  predicted  gas  enthalpies.  The  sonic  flow  section  of  the  code 
was  bypassed,  and  the  experimental  gas  flow  rates  were  input.  The  resultant 
performance  predictions  are  shown  in  Figures  30  through  32  and  Table  10  (see 
Table  9 for  nomenclature). 

The  key  predictions  were  the  total  enetgy  losses  since  the  experimental 
flow  rates  were  used  and  the  arc  length  was  selected  to  best  correlate  the 
experimental  arc  voltages  and  arc  tracks.  Although  the  magnitudes  of  the 
energy  loss  predictions  (Figure  32)  were  not  significantly  in  error,  particu- 
larly in  the  middle  pressure  range,  the  functional  dependence  on  pressure 
appears  high.  The  experimental  losses  Increased  as  the  0.57  power  of  the 
arc  pressure  while  AIRARC  predicted  an  exponential  dependence  of  approximately 
0.92.  This  is  reflected  in  the  predicted  enthalpies  which  exceed  the  measured 
values  by  as  much  as  30%  at  the  10  atm  pressure  level.  At  higher  pressures, 
the  enthalpy  predictions  were  much  better,  partially  because  of  small  over- 
prediction  of  the  arc  voltages 
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TABLE  9 AIRARC  OUTPUT  NOMENCLATURE 


Parameter 

Definition 

P 

ATM 

Arc  pressure  (atm) 

I 

AMP 

Arc  current  (A) 

DIAC 

IN 

Constrictor  diam  (in.) 

ARCL 

IN 

Arc  length  (in.) 

OlAT 

IN 

Nozzle  throat  diam  (in.) 

BGT 

DEGR 

Bulk  gas  temperature  (°R) 

V 

VOLT 

Arc  voltage  (V) 

PIN 

KW 

Air  power  input  (kW) 

GASP 

KW 

Power  to  gas  (kW) 

E 

V/IN 

Arc  potential  gradient  (V/m.) 

RAO 

B/F2S 

9 

Radiant  heat  flux  to  wall  in  arc  region  (Btu/ft  s) 

CNV 

B/F2S 

Convective  heat  flux  to  wall  (Btu/ft  s) 

QE 

B/S 

Energy  loss  to  nozzle  (Btu/s) 

M 

LBM/S 

Air  fiow  rate  (Ib/s) 

H 

B/LB 

Stagnation  enthalpy  (Btu/lb) 

EFF 

PERC 

Thermal  efficiency  (%) 

ARCT 

DEGR 

Arc  temperature  (°R) 

ARCD 

IN 

Arc  diam  (in.) 

I 

V: 

'I 

j 

I 

f 

i 

t 

1 

1 

s 

3 

] 

j 

<5 


I 
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p 

ATM 

V 

1 

AMP 

PI^ 

TABLE  10  PREDICTED  N-250  PERFORMANCE  BY  AIRARC  I 

DIAC  AriCL  DIA  T BGT 

I'i  I < DEG  H 

GASP  E RAD  C^V  QE  M ^ EFF  ARC  T 

ARC  D 

VOLTS 

KW 

Kw  ' 

V/I  B/F2S 

3/S 

LB/S  B/LB  PErtC 

DEG  R 

IM 

100 

3137 

240 

753 

0.375 

460 

7.0  0. 162 
448  2847 

8500 
1292  29 

0.  146 

3150 

61 

23093 

0.073 

100 

2783 

320 

391 

0.375 

531 

7.0  0.  162 
398  3787 

9606 
1371  34 

0.  137 

38  38 

60 

228  01 

0.092 

100 

2565 

400 

1026 

0.375 

602 

7.0  0.  162 
366  4673 

10537 
1446  39 

0.  132 

4470 

59 

22570 

0.  no 

75 
238  0 

240 

571 

0.375 

350 

6*5  0.  1 62 
366  2276 

8639 
1015  22 

0.  107 

3259 

61 

22987 

0.030 

75 

2126 

320 
68  0 

0.375 

409 

6.5  0.  162 
327  2976 

9727 
1081  28 

0.  102 

3975 

60 

22672 

0.  101 

75 

1954 

400 

782 

0.375 

461 

6.5  0*1 62 
301  3704 

10682 
1130  32 

0.09  7 

4654 

59 

22458 

0.121 

50 

1685 

240 

404 

0.375 

250 

6.0  0.  162 
281  1614 

8851 
736  13 

0.  071 

3516 

62 

22790 

0.091 

50 

1508 

320 

483 

0.375 

294 

6.0  0* 162 
251  2140 

999  1 
787  22 

0.  068 

4286 

61 

22520 

0.  114 

50 

1388 

400 

555 

0.375 

332 

6.0  0.  162 
231  2700 

109  46 
816  25 

0.  065 

5045 

60 

22345 

0.  136 

25 
99  1 

240 

238 

0.375 

152 

5.5  0.  162 
180  857 

9672 
451  13 

0.  035 

4252 

64 

22334 

0.116 

25 

875 

320 

280 

0.375 

176 

5.5  0.  162 
159  1169 

109  69 
463  15 

0.  032 

5347 

63 

22103 

0.  146 

25 

807 

400 

323 

0.375 

2G0 

5.5  0.1  62 
147  1508 

12023 
470  18 

0.  030 

6441 

62 

2199  1 

0.  172 

15 

697 

240 

167 

0.375 

112 

5.3  0.  162 
132  507 

10549 
319  10 

0.  OPl 

5161 

67 

21878 

0.  140 

15 

619 

320 

198 

0.375 

129 

5.3  0.162 
1 17  740 

11849 
321  12 

0.019 

6549 

65 

21811 

0.  173 

15 

582 

400 

233 

0.375 

148 

5.3  0.  162 
110  1029 

12866 
326  14 

0.  018 

7909 

64 

21912 

0.197 

Seo  Table  9 for  nomenclature 
See  Table  A 3 for  experimental  data 
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Substitution  of  a radiation  heat. transfer  correlation  (derived  in  an 
earlier  electrode  heat  transfer  study for  the  volumetric  radiation  of  Yos, 
yielded  much  closer  correlation  of  the  electrode  energy  losses  in  the  N-250 
operating  on  air.  Figure  45  and  Table  11  compare  the  empirical  version  of 
AIRARC  CAIRARC  II)  to  the  N-250  electrode  energy  losses.  The  functional 
dependence  on  pressure  is  closely  correlated,  and  the  magnitudes  are  quite 
accurate  over  the  full  range  of  the  data.  The  arc  voltage  and  resultant 
enthalpy  are  also  correlated  by  the  revised  code  with  improved  accuracy,  as 
shovm  by  comparing  Tables  11  and  23. 

TABLE  1 1 PREDICTED  N-250  PERFORMANCE  BY  AJRARC  H 


p 

I 

01  AC 

AHCL 

DIA  T 

3GT 

ATrt 

AMP 

IM 

I ^ 

' DEG  R 

V 

PI  V 

GASP 

£ 

HAD  C9V  QE 

6 

EFF 

ARC  T 

ARC  D 

VOLT 

KW 

KV 

v/i  g 

B/F2S  B/S 

LB>i/S  B/LB  PERC 

DEG  R 

IV  ' 

100 

240 

0*37 

7.0 

0*  162 

6572* 

2662 

666*  V 

467*  1 

409 

1764 

1 166 

29 

0*  146 

3194 

63*0 

27048 

0*05 

100 

320 

0*37 

7.0 

0*  162 

9676* 

2527 

606*5 

535*7 

361 

2453 

1271 

34 

0*  137 

3875 

66*3 

268  51 

0*07 

100 

400 

0.37 

7*0 

0*  162 

10564* 

232V 

V3l*5 

604.9 

333 

3170 

1346 

39 

0*  132 

4493 

64*9 

26633 

0*03 

75 

240 

0*37 

6*5 

0*  162 

6774* 

2237 

536*6 

361* 

344 

1528 

940 

23 

0*107 

3362 

67*3 

26648 

0*  06 

75 

320 

0*37 

6*5 

0*  162 

9669* 

IV92 

637.6 

416*  0 

307 

2126 

1012 

29 

0*  102 

4061 

65*6 

26440 

0*  07 

75 

400 

0*37 

6*5 

0*162 

10612* 

1632 

732.V 

470*5 

262 

2746 

1064 

33 

0*097 

4743 

64*2 

26283 

0*  09 

50 

240 

0*37 

6*0 

0*  162 

9000* 

16A7 

3V5*4 

256*6 

275 

1246 

664 

25 

0*071 

3605 

64*9 

26170 

0.  07 

50 

320 

0*37 

6*0 

0*  162 

10199* 

U75 

471. V 

306*3 

246 

1736 

744 

23 

0.  066 

4460 

64*9 

25951 

0*08 

bU 

400 

0*37 

6*  0 

0*  162 

1 1172* 

135V 

543*b 

345*7 

226 

2247 

777 

26 

0.065 

5245 

63*6 

25793 

0*  10 

25 

240 

0*37 

5*5 

0*  162 

993U* 

1040 

24V  *6 

163*3 

169 

664 

428 

14 

0*035 

4550 

65*4 

25192 

0*  09 

25 

320 

0*37 

5*5 

0*  162 

11276* 

V22 

295*0 

168*5 

166 

1232 

443 

16 

0*  032 

5709 

63*9 

25017 

0*  1 1 

25 

400 

0*37 

5*5 

0*  >2 

12376* 

651 

340*6 

213*3 

155 

159, 

454 

19 

0.030 

6666 

62*6 

24699 

0*  13 

15 

240 

0*37 

5*3 

0*  162 

10962* 

772 

165*2 

120*5 

146 

666 

307 

1 1 

0*  021 

5564 

65*  1 

24303 

0*  1 1 

15 

320 

0*37 

5*3 

0*  162 

12325* 

66  V 

220*6 

140*6 

130 

956 

312 

14 

0*019 

7146 

63*8 

24196 

0*  13 

15 

400 

0*37 

5*3 

0*  162 

13379* 

647 

259*0 

163*2 

122 

1243 

318 

16 

0*016 

6716 

63*0 

24145 

0.  15 

Seo  , eblo  9 for  nomenclature 
See  Table  A-3  for  expenmontat  data 
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Figure  45  Comparison  of  AtRARCHand  N-250  electrode  energy  losses 


Further  verification  of  fhe  scaling  accuracy  of  AIRARC  was  gained  through 
comparisons  with  actual  DET  arc  heater  data.  The  DET  heater  is  a scale  1. 73 
Huels-type  unit  with  a 0.5625  in.  diam  throat.  Two  sets  of  data  were  available 
for  pressures  from  29  to  93  atm  and  arc  currents  from  570  to  750  A. 


AEOC-TR-76-25 


The  first  set  of  data  was  acquired  using  a 34  in.  long  cathode.  The 
second  set  used  24,  27  and  30  in.  long  cathodes.  There  were  significant 
differences  in  the  performance  of  the  heater.  The  new  shorter  electrodes 
yielded  higher  arc  voltages  and  higher  enthalpies  from  improved  thermal 
efficiency.  Table  12  compares  AIRARC  I to  the  higher  performance  data 
(shorter  cathodes).  The  arc  length  selected  as  based  on  Figures  39  and  40 
where,  for  an  N = 1.73  arc  heater  and  a scaled  nozzle  throat  diameter  of 
d*/N  = 0.324  in.,  the  scaled  arc  length  (L/N  = 12  in.)  must  be  increased  by 
18%  for  a good  voltage  match.  This  yielded  a 24.5  in.  arc  length.  The 
resultant  arc  voltage,  enthalpy  and  thermal  efficiency  are  in  good  agreement 
over  the  complete  test  matrix. 

TABLE  12  COMPARISON  OF  AIRARC  I AND  DET  ARC  HEATER  DATA  (SHORT  CATHODES) 


p 

AM 

I 

A4P 

D1  AC 
I 4 

AriCL 

14 

DIA  I 
I 4 

BGl 
DEG  H 

V 

VOL  IS 

PI  4 
KW 

GASP 

r<W 

E 

V/I 

HAD 

B/F2S 

C4V  GE 

3/S 

A '{  EFF 

LB/S  B/L3  PERC 

ARC  T 
DEG  R 

ARC  D 
14 

30 

3AJ1 

(3900) 

740 

2539 

U300 

1506 

24*  5 
140 

0.  562 
667 

7617 
336  73 

0*504 

2966 

(3330) 

59 

(58) 

21569 

0*252 

2<5 

32t>3 

(3650) 

750 

2440 

1*300 

1476 

24-5 

133 

0*562 

642 

6 026 
329  73 

0*474 

3109 

(3360) 

61 

(59) 

21496 

0*262 

29 

32d5 

(3470) 

765 

2513 

1*300 

1542 

24*  5 
134 

0*562 

661 

6120 
336  76 

0*464 

3i72 

(3190) 

61 

(59) 

21473 

0*264 

5J 

5290 

(5610) 

700 

3703 

:*300 

2167 

24*5 

216 

0*562 

1344 

6777 
476  93 

0.936 

2356 

(2530) 

59 

(60) 

2206  7 

0*  192 

51 

5312 

(5540) 

715 

3796 

1*300 

2264 

24*5 

217 

0*562 

1377 

6634 
495  97 

0*969 

2363 

(2480) 

60 

(50) 

22*03 

0*  193 

52 

5239 

(5210) 

730 

3625 

1*300 

2326 

24-5 

214 

0*  562 
1406 

7006 
501  96 

0*946 

2473 

(2510) 

61 

(62) 

22013 

0*  193 

6A 

6225 

(6220) 

710 

4420 

1*300 
266  0 

24*5 

254 

0.  562 
1650 

6622 
573  106 

1*211 

2246 

(2360) 

61 

(64) 

22222 

0*  173 

76 

7359 

(7340) 

650 

4764 

1*300 

2925 

24*5 

300 

0.  562 
1606 

6154 
621  112 

1*463 

2014 

(2000) 

61 

(61) 

22227 

0.  153 

76 

7332 

(7490) 

675 

4949 

1*300 

3060 

24*5 

299 

0*562 

1661 

6265 
646  116 

1*523 

2062 

(2130) 

62 

(63) 

22205 

0-  162 

92 

6704 

625 

5440 

1*300 

3234 

24*5 

355 

0*562 

2010 

566  0 
672  116 

1*636 

16  13 

59 

22655 

0*  133 

(8480)  (1870)  (63) 

Seo  Table  9 for  nomor^claturo 
( ) DET  arc  heater  d&ta^ 
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Table  13  compares  AIRARC  II,  the  empirical  version  of  AIRARC,  with  the 
DET  arc  heater  data  using  the  long  (34  in.)  cathode.  The  scaled  arc  length 
was  the  same  as  that  used  with  AIRARC  I,  24.5  in.  The  agreement  is  good 
over  the  full  range  of  the  data  for  arc  voltage,  enthalpy  and  thermal 
efficiency.  Thus,  both  versions  of  AIRARC  can  be  used  with  reasonable 
confidence  for  scaling  air  arc  heaters  to  N = 1.73. 


TABLE  13  COMPARISON  OF  AIRARC  H AND  DET  ARC  HEATER  DATA  (LONG  CATHODE) 


p 

ATrt 

I 

AMP 

D 1 Hv 
19 

AHCL 

X9 

DIA  T 
19 

BGT 
DEG  R 

V 

VOLT 

Piy 

KW 

GASP 

KW 

E 

9/19 

PAD  C9V  QE  M K EFF 

B/F2S  3/S  LBM/S  B/LB  PErtC 

ARC  T 
DEG  R 

ARC  D 
19 

50 

(4770) 

680 

3399 

U30 

2029 

24«5 

204 

0.562 

U80 

6314. 

406  84 

0.967 

2136  59.7 
(1910)  (55) 

26677 

0.  13 

52 

^<$95 

(4870) 

750 

3672 

1*30 

2157 

24.5 

200 

0.562 

1345 

6498. 

422  91 

0.975 

2245  53.7 
(2120)  (55) 

26688 

0.  14 

77 
659  0 
(6660) 

655 

4317 

1*30 

2657 

24.5 

269 

0.562 

1404 

5635. 

512  100 

1.523 

1797  61.6 
(1730)  (58) 

27329 

0.  10 

77 

6394 

(6500) 

705 

4508 

1*30 

2740 

24.5 

261 

0.562 

1527 

5813. 

525  104 

1.503 

1372  60.3 
(1810)  (57) 

27284 

0.  11 

93 

773tJ 

(8150) 

590 

4565 

1*30 

2894 

24.  5 
316 

0.562 

1371 

5147. 

544  104 

1.887 

1596  63.4 
(1680)  (63) 

27673 

0.09 

92 

7351 

670 

4925 

U30 

3066 

24.5 

300 

0.562 

1576 

5420. 

572  . 1 12 

1.856 

1708  62.3 

27593 

0.  10 

(7720)  (1750)  (60) 

Seo  Table  9 for  nomenclatur') 

( ) DET  arc  heater  data^ 


5.3  HEARC  COMPARISONS  AND  SCALING 

The  results  of  the  algorithm  for  helium  are  presented  in  Figures  20,  22, 
32,  24,  and  46  as  a function  of  pressure  with  current  as  a parameter.  For 
comparison,  experimental  results  are  also  shown.  The  predicted  performance 
is  also  showh  in  tabular  form  in  Table  15  (see  Table  14  for  nomenclature). 

The  experimental  results  are  tabulated  in  Table  22. 

The  resulting  **matched"  voltages  are  shown  in  Figure  22.  The  predicted 
enthalpy  and  mass  flow  rates  compared  to  experiment  are  shown  in  Figures  20 
and  36.  Note  that  where  the  calculated  enthalpy  deviates  slightly  (3%)  from 
experiment,  such  as  at  the  100  atm  400  A point,  the  calculated  flow  rate  also 
deviates  slightly  (7%).  If  the  calculated  gas  power  is  divided  by  the 
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experimental  f.low  rate,  the  resulting  enthalpy  is  closer  to  the  experimental 
value.  This  possibly  indicates  that  the  effective  nozzle  throat  diameter  in 
the  experiment  was  slightly  less  than  0.162  in.,  the  value  used  by  the 
algorithm  sonic  flow  equation  for  determining  mass  flow  ra^e*^. 


Figure  46  Comparison  of  measured  and  HEARC  predicted  mass  flow  rates 


The  rear  electrode  power  loss  is  shown  in  Figure  23.  These  were 
essentially  "matched**  to  experiment  by  adjusting  the  constant  in  Equation  (38) 
and  by  varying  the  rear  arc  length. 

The  front  electrode  losses  are  shown  in  Figure  24.  Note  that  the  total 
power  loss  of  the  front  plus  rear  electrode  is  in  better  agreement  with 
experiment  than  either  one  separately.  This  indicates  compensating  errors  in 
predictions  for  front  and  rear  electrodes  and  perhaps  also  the  difficulty  in 
determining  the  fraction  of  the  total  arc  length  within  each  electrode. 

The  calculated  efficiencies,  Table  15,  are  in  good  agreement  with 
experimental  values. 

The  semi-empirical  scaling  code  HEARC,  having  been  correlated  with 
the  N-250  experimental  data,  has  been  used  to  predict  the  performance  of 
larger  scale  heaters. 
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TABLE  14  HEARC  AND  HYARC  OUTPUT  NOMENCLATURE 


Parameter 

Definition 

P 

ATM 

Arc  pressure  (atm) 

AMPS 

Arc  current  (A) 

DCONS 

IN 

Constrictor  diam  (in.) 

DTHRT 

IN 

Nozzle  throat  diam  (in.) 

LCONS 

IN 

Front  constrictor  length  (In.) 

ARCL 

IN 

Total  arc  length  (In.) 

FRNT 

IN 

Arc  length  in  front  electrode  (in.) 

REAR 

IN 

Arc  length  in  rear  e!ectiode  (in.) 

TSTAG 

DECK 

Stagnation  temperature  (K) 

QREAR 

Power  loss  to  rear  electrode  (Btu/s) 

QFRNT 

Power  loss  to  front  electrode  (Btu/s) 

ONOZ 

Power  loss  to  nozzle  (Btu/s) 

VOLTS 

Arc  voltage  (V) 

PIN 

KW 

Electrical  power  input  to  arc  (kW) 

GASP 

KW 

Power  added  to  gas  (kW) 

GRAD 

V/IN 

Arc  column  voltage  gradient  (V/in.) 

QRAD 

BTU/FT2S 

o 

Radiative  heat  flux  to  constrictor  wall  (Btu/ft  s) 

QCONV 

Convective  heat  flux  to  constrictor  wall  (Btu/ft  s) 

WDOT 

LB/S 

Mass  flow  rate  (Ib/s) 

ENTH 

B/LB 

Stagnation  enthalpy  (Btu/lb) 

EFF 

PERC 

Efficiency  (%) 

TARC 

Arc  column  temperature  (K) 

DEGK 

DARC 

IN 

Arr'  column  diam  (in.) 
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TABLE  15  PREDICTED  N-250  PERFORMANCE  BY  HEARC 


HELIUM  ARC  HEATER:  SCALE  N = 0.5 


t'tiES 

DCO^S 

DMrfT 

LCO^S 

AtiCL 

Fil9r 

rtKAR 

TSTAG 

QREAR 

QFR'JT 

Q>J0Z 

AH 

I >1 

IM 

I 

1 3 

I ^ 

DEG  < 

BTU/SEC- 

yJOL  i‘b 

GAS/^ 

GiiAD 

QrtAD 

QC0>JV 

wDOr 

£>m 

EFF 

TAKC 

DAHC 

KW 

KW 

O/I  ^ 

3TU/F*r 

•*2“S 

L3'l/S 

B/LB 

PERC 

DEG  < 

IM 

as 

240 

.375 

.162 

5.0 

4.5 

2.3 

2.2 

333  3 

13 

47 

12 

232 

147 

217 

731 

303 

.020 

7566 

63 

16512 

. 113 

as 

320 

.375 

.162 

5.0 

4.5 

2.3 

2.2 

4029 

19 

53 

21 

<$46 

270 

16S 

l:>0 

1 1 13 

925 

.019 

3395 

60 

16412 

• 143 

as 

400 

.375 

. 162 

5.0 

4.5 

2.3 

2.2 

47  17 

26 

63 

24 

766 

306 

17> 

172 

1433 

1013 

.017 

10530 

53 

16362 

. 169 

so 

240 

.375 

• 162 

5.0 

4.7 

2.0 

2.7 

2313 

34 

77 

24 

16:^S 

406 

262 

35:^ 

1565 

1259 

. 044 

6233 

64 

17012 

.033 

so 

320 

.375 

. 162 

5.0 

4.7 

2.0 

2.7 

3279 

43 

37 

23 

i4ya 

477 

2^2 

316 

2273 

1449 

.042 

7320 

61 

16925 

. 110 

so 

400 

.375 

. 162 

5.0 

4.7 

2.0 

2.7 

3746 

64 

1 15 

32 

USl 

S40 

3lb 

236 

2352 

1624 

.039 

3362 

53 

16737 

. 133 

vs 

240 

.375 

. 162 

5.0 

5. 0 

2.  1 

2.9 

2532 

5! 

103 

30 

2376 

S70 

375 

477 

2152 

1632 

. 070 

5733 

65 

17250 

.077 

75 

320 

.375 

« 162 

5.0 

5. 0 

2.  1 

2.9 

3044 

66 

126 

35 

205:^ 

6S3 

417 

413 

2733 

1333 

.065 

6735 

63 

17012 

. 100 

75 

400 

.375 

. 162 

5.0 

5.0 

2.  1 

2.9 

3417 

91 

153 

39 

1363 

747 

447 

375 

3357 

2145 

.061 

7629 

59 

17012 

. 113 

100 

240 

.375 

• 162 

5.0 

5.2 

2.2 

3. 0 

2417 

70 

130 

34 

3014 

723 

475 

574 

2332 

1923 

.035 

5395 

65 

17750 

.065 

100 

320 

.375 

. 162 

5.0 

5.2 

2.2 

3.0 

2316 

32 

161 

40 

2635 

343 

532 

502 

3745 

223  1 

* .090 

6233 

63 

17437 

.036 

100 

400 

.3/5 

« 162 

5.0 

5.2 

2.2 

3.0 

3195 

117 

192 

44 

2371 

V43 

574 

452 

4725 

2600 

.034 

7134 

60 

17262 

. 105 

See  Teblo  14  for  nomenclature 
See  Tcble  A 2 for  experimental  data 


Results  for  scale  N = 1 and  N ==  1. 73  heaters  corresponding  to  constrictor 
dlarneters  f/'.  0.75  and  1*3  in.,  are  shown  in  Tables  16  and  17.  Arc  lengths, 
front  and  rear,  were  scaled  from  those  of  the  N-250,  (N  = 0.5)  heater  as: 


" (0.5)  H1250  • 

The  front  constrictor  lengths  were  10  and  17.3  in.  for  N = 1 and  1.73 
respectively.  Nozzle  diameters  were  0.324  and  0.562  in.  which  were  scaled  as: 


3 

i 

3 

'I 

i 

I 

* 
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d 


* 


^*N250 


(47) 


Arc  currents  were  scaled  up  approximately  as; 


I = 


(48) 


nozzle  losses  were  evaluated  from  the  expression  shown  in  Figure  25- 


TABLE  16  SCALED  HELIUM  ARC  HEATER  PERFORMANCE  (N  = 1.0) 


rfELiJrt  AHC  rfEATEa*  SCALE  I 


i 


^HES 

AMPS 

DC0>JS 

DHRT 

LCO'IS 

ArfCL 

ATM 

1^ 

13 

ly 

yJOLlS 

PI  ^ 

GASP 

GRAD 

QRAD 

QCO^V 

liW 

KW 

BTU/FT»*2-S 

as 

400 

• 750 

• 324 

10.0 

8.9 

1767 

706 

493 

198 

528 

546 

25 

600 

.750 

• 324 

10.0 

3.9 

1A36 

861 

567 

161 

861 

664 

25 

800 

• 750 

.324 

10.0 

8.9 

1223 

978 

640 

137 

937 

761 

50 

400 

-750 

• 324 

10.0 

9.4 

3160 

1264 

875 

335 

1217 

8 09 

50 

600 

.750 

.324 

10.0 

9.4 

257*» 

1547 

1026 

273 

1705 

1037 

50 

800 

.750 

.324 

10.  0 

9.4 

2222 

1778 

1147 

235 

2117 

1219 

75 

400 

.750 

• 324 

lU.O 

10.0 

A4l23 

176^ 

1243 

444 

1629 

1036 

75 

600 

.750 

• 324 

10.0 

10.0 

3673 

2203 

1424 

363 

2790 

1301 

75 

800 

.750 

• 324 

10.  0 

10.0 

3212 

2569 

1587 

322 

3649 

1546 

100 

400 

.750 

• 324 

10.  D 

10.5 

5681 

2272 

1578 

541 

2308 

1201 

100 

600 

.750 

• 324 

10.0 

10.5 

A62A 

2774 

18  54 

440 

3080 

1595 

100 

800 

.750 

• 324 

10.0 

10.5 

A087 

3270 

2014 

38  9 

4646 

1320 

T*bl«  14  for  1 

nomtnclatur* 

REAR  TSTAG  QREAR  QFR>IT  0^02 


I M 

19 

DEG  K 

BT 

U/SFC-- 

^Dor 

LBM/S 

E9rrf 

B/LB 

EFF 

PERC 

TARC 
DEG  < 

DARC 

IM 

4.5 

.095 

4.  4 
5580 

2499 

69 

37 

16125 

128 
• 166 

35 

4.5 

.034 

4.4 

7043 

3154 

65 

61 

16000 

172 
• 231 

44 

4.5 
. 077 

4.  4 
8550 

3330 

65 

70 

15937 

197 

.293 

52 

4.  1 
.206 

5.4 

4696 

2103 

69 

106 

16712 

213 
. 125 

47 

4.  1 
. 186 

5.4 

5913 

2643 

66 

149 

16262 

233 
. 187 

60 

4.  1 
. 167 

5.4 

7184 

3213 

64 

136 

16000 

340 

.250 

M 

4.2 
• 321 

5. 3 
4354 

1950 

70 

154 

17012 

230 
• 103 

57 

4.2 
• 295 

5.8 

5252 

2352 

64 

265 

16837 

402 

.150 

70 

4.2 

.271 

5.3 

6221 

2737 

61 

346 

16550 

500 
• 199 

83 

4.4 
• 444 

6.  1 
4037 

1303 

69 

223 

17562 

364 
• 090 

64 

4.4 

.401 

6.  1 
5049 

2262 

66 

305 

17012 

484 
. 139 

32 

4.4 

.373 

6.  1 
5783 

2592 

61 

460 

17012 

634 
. 171 

94 

89 
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TABLE  17  SCALED  HELIUM  ARC  HEATER  PERFORMANCE  (N  = 1.73) 


rfELIUM 

[ AtiC 

HEATEK 

5 SCALE  >1-1. 

73 

AiHS 

DCO^S 

DT3KT 

LC0>JS 

ArtCL 

AT>1 

I 9 

IM 

IM 

V/OLTS 

r'I'J 

GASt^ 

GKAD 

QKAD 

QCOnIV 

Kb 

KW 

9/1  y 

BTU/FT 

**2-S 

25 

692 

1.300 

.562 

17.3 

15.4 

261A 

1309 

1336 

169 

411 

416 

25 

lOJd 

1.300 

.562 

17.3 

15.4 

2121 

2202 

1543 

137 

669 

513 

25 

133A 

1.300 

.562 

17.3 

15.4 

16S0 

2561 

1713 

119 

965 

537 

50 

692 

1.300 

.562 

17.3 

16.3 

A505 

3325 

2375 

293 

1 113 

603 

50 

1033 

1.300 

.562 

17.3 

16.3 

J756 

3399 

23  00 

229 

1 154 

3 05 

50 

133A 

1.300 

.562 

17.3 

16.3 

J279 

A539 

3113 

200 

1650 

929 

75 

692 

1.300 

.562 

17.3 

17.3 

6d3;^ 

A733 

3316 

396 

1733 

746 

75 

1033 

1.300 

.562 

17.3 

17.3 

5623 

5337 

3907 

325 

2471 

995 

75 

133A 

1.300 

.562 

17.3 

17.3 

AdVJ 

6773 

4374 

233 

3113 

1210 

100 

692 

1.300 

.562 

17.3 

13.2 

d6J2 

5973 

423  7 

474 

2025 

39  1 

100 

1033 

1.300 

» d62 

17.3 

13.2 

7225 

7500 

49  16 

397 

3379 

1143 

100 

1334 

1.300 

.562 

17.3 

13.2 

6362 

33  05 

5503 

349 

4420 

1331 

See  Table  14  for 

nomenclature 

FH^IT 

HEAR 

TSTAG 

QREAH 

QFR>JT 

Q>10Z 

I 'J 

I'J 

DEG  < 

era/SEC- 

WDOr 

E'JTH 

EFF 

TARC 

DARC 

L33/S 

B/L3 

FERC 

DEG  < 

7.3 

7.6 

2127 

33 

296 

63 

.310 

4749 

73 

15975 

.244 

7.3 

7.6 

2655 

144 

400 

79 

.273 

5927 

70 

15312 

.343 

7.3 

7.6 

3143 

207 

502 

93 

.255 

7C27 

66 

15737 

.431 

7.0 

9.3 

13  09 

293 

520 

35 

.663 

4033 

71 

16337 

.191 

7.0 

9.3 

2302 

304 

625 

1 1 1 

.594 

5139 

71 

15925 

.293 

7.0 

9*3 

2697 

435 

735 

129 

. 551 

6022 

63 

15337 

.370 

7.2 

10.  1 

1637 

510 

731 

100 

1.053 

3654 

70 

1697  5 

.152 

7.2 

10.  1 

2039 

705 

993 

130 

.954 

4553 

66 

16337 

.236 

7.2 

10.  1 

2417 

339 

1229 

154 

.377 

5397 

64 

16100 

.317 

7.7 

10.5 

157  9 

603 

373 

1 16 

1.423 

3526 

71 

17100 

. 141 

7.7 

10.5 

1336 

1036 

1297 

145 

1.316 

421  1 

65 

17000 

. 194 

7.7 

10.  5 

221  1 

1316 

1641 

172 

1.223 

4936 

62 

16625 

.261 
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5 . 4 HYARC  COMPARISONS  AND  SCALING 


The  results  of  the  algorithm  for  hydrogen  are  presented  in  Figures  5, 
12  and  13  and  in  tabular  form  in  Table  18  as  a function  of  pressure  with 
current  as  a parameter.  For  comparison  the  experimental  results  are  shown 
in  Table  A-1. 


TABLE  18  PREDICTED  N-250  PERFORMANCE  BY  HYARC 


HYDROGEN  ARC  HEATER:  SCALE  N = 0.5 


DCO^b 

D16H‘ 

Lcogs 

AHCL 

F«gr 

rtEArt 

TSFAG 

QfiEA6 

QFrtgr 

Qgo-c 

Af'! 

in 

in 

in 

IN 

IN 

IH 

DEG  K 

BT'J/SEC-* 

\/0Lib 

t'l^ 

UA6P 

GrCAD 

<3tfAD 

jcogo 

VD3  r 

Egn 

EFF 

TARC 

DARC 

rCw 

KW 

V/I  M 

31  J/r'I 

*•2-5  LBA/S 

B/LB 

PE6C 

DEG  < 

I g 

10 

240 

.37b 

• 162 

b.  0 

6.2 

3.6 

2.6 

3553 

12 

42 

21 

230 

146 

154 

600 

607 

.004 

35377 

64 

11562 

.160 

10 

320 

.375 

. 162 

b.O 

5.4 

3.4 

2.0 

3656 

15 

52 

23 

<504 

2b7 

160 

147 

725 

667 

.004 

36727 

62 

11675 

.206 

10 

400 

.37b 

• 162 

5.0 

4.6 

3.2 

1.6 

3755 

15 

61 

26 

272 

163 

Ibl 

1 155 

760 

. 004 

42244 

62 

11667 

.226 

lb 

240 

• 37b 

• 162 

b.  0 

6.7 

3.7 

2.6 

3472 

17 

53 

27 

iJb:5 

326 

21b 

201 

73  0 

625 

.007 

31646 

66 

1 1562 

. 162 

16 

320 

.376 

• 162 

b.  0 

5.6 

3.7 

2.  1 

3635 

21 

70 

23 

lObJ 

337 

207 

16  1 

1230 

62b 

.006 

36636 

62 

11625 

.175 

1 6 

400 

.376 

. 162 

5.0 

5.2 

3.5 

1.7 

3607 

23 

63 

27 

676 

3b0 

206 

166 

1716 

6367 

. 005 

4246  1 

56 

1 166? 

.224 

20 

240 

• 37b 

• 162 

5.  0 

7.3 

4.3 

3.0 

3566 

24 

71 

30 

IdJJ 

366 

23b 

211 

777 

6d7 

.007 

33633 

63 

11475 

. 165 

20 

320 

.37b 

• 162 

5.  0 

6.2 

4.  0 

2.2 

3656 

2o 

66 

32 

1262 

410 

2b3 

205 

1467 

7 67' 

.007 

36136 

61 

1 1562 

. 170 

20 

400 

• 37b 

* 162 

b.O 

5.6 

3.6 

1.7 

3737 

27 

107 

35 

1 131 

4b2 

27  0 

203 

2071 

1040 

.007 

36467 

57 

11662 

.206 

2b 

240 

.37b 

• 162 

b.O 

7.6 

4.6 

3.  1 

3473 

27 

70 

37 

2034 

466 

321 

261 

1156 

1145 

.01  1 

30545 

65 

11550 

• 143 

2b 

320 

.37b 

. 1 62 

5.0 

6.7 

4.4 

2.  3 

3626 

32 

106 

33 

ib74 

b03 

314 

236 

1732 

1147 

.007 

34525 

62 

1 1562 

.160 

2b 

400 

.37b 

• 162 

b.O 

5.7 

4.  1 

1.6 

3741 

35 

130 

40 

1 3b2 

b40 

324 

227 

2365 

1201 

.007 

37524 

57 

11625 

r?02 

bO 

240 

.37b 

• 162 

6.0 

10.4 

6.4 

4.  0 

3466 

65 

230 

57 

4046 

771 

576 

370 

2010 

17  17 

• 021 

2677  1 

61 

11600 

. 127 

bO 

320 

375 

• 162 

6.0 

6 . 6 

6.0 

2*6 

3565 

65 

271 

63 

3273 

10b3 

631 

373 

2636 

2021 

• 021 

30567 

57 

1 1612 

• 147 

bO 

400 

.3  75 

• 162 

6.0 

7.6 

5.7 

2.  1 

3577 

67 

321 

65 

261b 

1126 

647 

360 

3731 

2100 

. 021 

31313 

57 

11667 

. 166 

Soo  Table  14  for  nomanclaturo 
See  Table  A-1  for  exponmontal  data 
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These  results  utilized  the  experimentally  measured  flow  rates  as  input 
to  the  algorithm  with  nozzle  sonic  flow  calculation  of  mass  flow  rate  by-passed. 
When  sonic  flow  is  not  bypassed  the  code  over-predicts  the  mass  flow  rates. 

This  is  most  likely  caused  by  the  fact  that  the  experimental  nozzle  throat 
diameter  late  in  the  tests  was  not  effectively  0.162  in.  due  to  clogging.  As 
can  be  seen  from  Figures  6 to  8,  the  experimental  sonic  flow  enthalpy 
disagreement  (with  the  energy  balance  value)  increases  as  the  running  time 
increases,  indicating  a changing  throat  size.  This  was  also  seen  in  the 
computed  results  with  mass  flow  calculations.'  The  disagreement  was  greater 
for  those  experimental  cases  where  the  run  time  was  longer,  i.e. , when  a run 
consisted  of  the  sequence  240,  320  and  400  A.  The  predictions  for  enthalpy 
for  the  10  atm  point  are  not  in  good  agreement  due  to  the  short  arc  conditions 
described  in  Section  4. 1, 

The  calculated  efficiencies  are  in  reasonable  agreement  with  experiment 
indicating  that  for  proper  voltage  (power)  predictions,  the  power  to  the  gas 
has  been  predicted  adequately.  The  semi-empirical  scaling  code  HYARC  was 
used  to  predict  the  performance  of  scale  N = 1 and  N = 1.73  heaters. 

Constrictor  lengths  utilized  were  24  and  36  in. , respectively.  Nozzle 
diameters  and  currents  were  scaled  in  the  same  manner  as  for  helium.  The 
nozzle  loss  expression  was  that  shown  in  Table  8.  The  front  and  rear  arc 
lengths  were: 


L 


R 


6 D , 
c 


(49) 


L 


f 


14.36 


p0.391 


(50) 


This  particular  form  for  corresponds  to  the  rear  arc  terminus  at  a spin 
coil  location  placed  six  constrictor  diameters  behind  the  inlet  housing. 

The  front  arc  length  corresponds  to  a total  arc  length  correlated  with  data 
for  a preliminary  version  of  HYARC,  These  arc  lengths  scale  directly.  Sonic 
flow  was  used  in  these  predictions  to  obtain  the  mass  flow  rate  since  no 
experimental  values  were  available. 

The  scaled  results  are  shovm  in  Tables  19  and  20.  Figures  47  and  48  show 
the  power  supply  requirements  for  these  predictions  along  with  the  character- 
istics of  the  MDRL  power  supply. 
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1 


TABLE  19  SCALED  HYDROGEN  ARC  HEATER  PERFORMANCE  (N  = 1.0) 

H^DHOGE^  AiiC  "lEAlEKS  SCALE  D*  = 0.32A  I ; LCO>JS=a^  I 

t'dES  Mt^S  DCOgi>  DHrtl  LCOgS  AriCL  Krfgi  rtEArt  TSTAG  2AEAA  QF^gi  Q^OE 

Afi  ig  ig  ig  ig  ig  ig  deg  < 3tu/sec 

X/OLIS  r»ig  GASt'  GriAD  JrtAD  OCOgV  wDJf  EgTH  EFF  TAHC  DArtC 

Kw  Kw  s;/ig  31J/FI**2-S  L3g/S  B/L3  PE.-^C  DEG  •<  ig 


2d 

AOO 

173/ 

. 750 
1 lOb 

.32A 

2>2 

24.  0 
>71 

1 5.  3 
335 

2b 

Jb2J 

600 

211A 

^7b0 

k2bb 

.33A 

350 

24.  0 
14>5 

14.  1 
>56 

2b 

JU25 

300 

2A20 

.750 

1331 

.334 

227 

24.  0 
2116 

13.3 

1064 

SO 

61bb 

AOU 

3502 

.7b0 

220> 

.324 
43  7 

24.0 

1533 

20.  0 
1443 

50 

600 

4I1A3 

.750 

2A67 

.324 

374 

24.0 

2533 

13.4 

1660 

50 

5:^10 

300 
A 7 23 

.750 

2677 

.324 

33> 

24.  0 
3634 

17.4 

1323 

75 

13J16 

AOO 

b326 

.750 

3272 

• 324 
563 

24.  0 
2541 

23.4 

1>>7 

7b 

10360 

600 

6316 

. 7 50 
37A6 

• 3 24 
47> 

24.  0 
33>3 

21.6 

2331 

7b 

3331 

300 
7 06A 

. 750 
A036 

• 324 
432 

24.  0 
4355 

30.4 

2553 

100 

1752A 

AOO 

7ooy 

. 7 50 
AA03 

.324 

663 

34.  0 
23  53 

26.2 

2557 

100 

13630 

600 

3303 

. /bO 
A>A6 

.324 

565 

24.  0 
41>3 

24.2 

2>37 

100 
1 1633 

300 

>306 

. 750 
5315 

.334 

50> 

34.0 

5>43 

22.  3 
3l>6 

See  Table  14  fo*  nomenclataro 


10.3 
. 070 

4.5 

16672 

2472 

61 

71 

11675 

4>5 
. 176 

79 

>.6 

.066 

4.  5 
1>>00 

23  03 
5> 

110 

11537 

60> 

.241 

94 

5.3 
. 062 

4.  5 
22320 

3040 

57 

155 
1 1500 

722 

.294 

106 

15.5 
. 136 

4.  5 
17170 

253> 

63 

1 16 
1 1662 

>70 
. 152 

133 

13.  > 
. 123 

4.5 

20012 

2341 

5> 

1 >0 
1 1537 

1242 

.205 

160 

12.  > 
. 122 

4.  5 
22625 

3071 

56 

271 
1 1550 

1494 

.251 

178 

13.  > 
.206 

4.5 
163  72 

2511 

61 

137 
1 1837 

1572 
• 130 

137 

17.  1 
. 1 >2 

4.  5 
20243 

2334 

60 

250 
1 157  5 

1367 
. 133 

223 

15. > 
. 133 

4.5 

22720 

3111 

57 

357 
1 1537 

2266 

.230 

246 

21.7 

.36> 

4.  5 
1 /3>> 

2567 

62 

210 
1 1737 

2020 
. 125 

239 

l>.7 

.251 

4.5 

30451 

2>22 

60 

303 
1 1575 

2504 
. 176 

273 

15.3 

*233 

4.  5 
22>34 

U63 

57 

437 
1 1525 

3033 

.213 

307 

i 

i 

I 

93  t 
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TABLE  20  SCALED  HYDROGEN  ARC  HEATER  PERFORMANCE  {N  = 1.73) 
rf/DrtOCE^  ArtC  HEAlE.-iS  SCALE  D*=0.b625  M.;  LCOVS»36 


Aif^S 

DCOgS 

DI-lKl 

Lcogs 

AHCL 

Ai.i 

1 V 

ig 

I g 

VOL  IS 

/^1>J 

GPiSt' 

G.-lAD 

OHAD 

jcogv 

KW 

v/i  g 

Di  J/FI 

♦♦2-S 

2b 

6^2 

1*300 

*562 

36*0  , 

26*4 

65S2 

A555 

3067 

24d 

751 

663 

25 

IDJS 

1*300 

*562 

36*0 

24*  4 

SISO 

5376 

3431 

212 

1249 

766 

25 

13S4 

1*300 

*562 

36*0 

23*  0 

AA2i 

6122 

3754 

192 

1761 

d53 

50 

6^2 

1*300 

* 562 

36*0 

34*  7 

UU50 

9031 

6165 

376 

1232 

1 IdO 

5'' 

103d 

1*300 

*562 

36*0 

32*  0 

1026A 

10654 

6d60 

320 

2056 

135d 

dO 

I3d4 

1*300 

* 562 

36*0 

30*2 

S7JS 

120:^3 

7441 

2d  9 

29  17 

149d 

75 

692 

1*300 

*562 

36.0 

40*6 

W669 

13611 

9397 

4d3 

160d 

1634 

75 

103d 

1*300 

*562 

36*0 

37*5 

1550A 

16093 

I049d 

413 

26d2 

1927 

75 

13d4 

1*300 

*562 

36.  0 

35*4 

13176 

ld235 

11312 

372 

3d55 

2114 

IDO 

692 

1*300 

*562 

36*0 

45*5 

261S0 

ldll6 

12592 

575 

1962 

2126 

100 

103d 

1*300 

*562 

36*0 

41*9 

20550 

21330 

13932 

490 

3290 

2437 

100 

1334 

1*300 

*562 

36*0 

39*6 

17A45 

24144 

15001 

440 

4691 

266d 

Sm  Tkota  14  for  nomtncloturo 


Frigr 

tfEAii 

ISl  AG 

3jiEAA 

3F3gi 

Qg02 

ig 

ig 

DEG  K 

BT  J/SEC-- 

wD/l 

Egn 

EFF 

fAac 

DA.3C 

LBg/S 

3/LB 

FEHC 

DEG  < 

I g 

ld*6 

7*3 

2243 

166 

1074 

169 

*224 

14301 

67 

11225 

*230 

16*6 

7*3 

2544 

276 

1370 

196 

*210 

17319 

63 

11162 

*373 

15*2 

7*3 

2733 

339 

1632 

221 

• 1 99 

19725 

61 

11112 

*465 

26*9 

7*3 

2321 

272 

2145 

297 

.430 

153J6 

63 

1 1212 

-243 

24*2 

7*3 

2613 

455 

2793 

343 

*406 

17331 

64 

11150 

*329 

22*4 

7*3 

2350 

645 

3332 

332 

*335 

20104 

61 

11100 

*406 

32*d 

7*3 

2369 

355 

3213 

419 

*639 

15743 

69 

1 1 137 

*224 

29*7 

7*3 

2663 

593 

4226 

434 

*606 

13215 

65 

1 1 125 

*303 

27*6 

7*3 

2393 

353 

5175 

534 

*577 

20335 

62 

! 1037 

*373 

37*7 

7*3 

2423 

434 

4269 

533 

*332 

16145 

69 

1 1 175 

*212 

34*  1 

7*3 

2720 

723 

5675 

609 

*737 

13579 

65 

11112 

*237 

31*d 

7.3 

2956 

1033 

6956 

672 

*749 

20775 

62 

11062 

*355 

94 
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Arc  curient  (A) 


Figure  47  HYARC  predicted  Huels-type  arc  heater  characteristics  50  atm  pressure 
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6,  CONCLUDING  REMARKS 


This  work  has  shown  that  significantly  higher  gas  stagnation  enthalpies 
can  be  achieved  with  a Huels-type  arc  heater  operating  on  hydrogen  and  helium 
than  are  possible  using  air.  Specifically,  hydrogen  enthalpies  of  five  to 
seven  times  those  on  air  can  be  achieved  at  pressures  to  50  atm,  and  helium 
enthalpies  41%  to  69%  higher  than  air  can  be  achieved  at  pressures  to  100  atm. 
This  helps  verify  the  original  hypothesis  that  the  AEDC  DET  can  achieve  higher 
dust-particle  velocities  using  these  gases  as  accelerators. 

A standard  Huels-type  arc  heater  can  be  operated  on  pure  hydrogen  at 
pressures  to  50  atm  and  enthalpies  to  40,000  Btu/lb  without  major  redesign. 
Significant  instrumentation  and  stability  problems  had  to  be  solved  before 
achieving  accurate  data  over  a wide  range  of  arc  pressures  and  arc  currents, 
but  no  arc  heater  modifications  are  necessary  other  than  gas  injector  sizing, 
spin  coil  insulation,  and  proper  nozzle  design.  Hydrogen  performance  data 
are  not  highly  repeatable.  The  arc-jet  appears  more  stable  on  hydrogen  than 
the  arc  current  and  voltage  fluctuations  indicate. 

A standard  Huels-type  arc  heater  can  be  operated  on  pure  helium  at 
pressures  to  100  atm  and  enthalpies  to  9700  Btu/lb  without  major  redesign. 

The  arc  is  quite  stable  and  the  data  are  repeatable. 

Electrode  erosion  and  thus  test  stream  contamination  are  a problem  on 
hydrogen  at  low  pressures  where  the  arc  is  not  properly  stabilized  in  the 
cathode  (front  electrode).  However,  hydrogen,  helium  and  air  all  exhibit  less 
than  0. 4%  test  stream  contamination  at  pressures  above  25  atm.  Anode  erosion 
is  severe  on  air  at  pressures  above  75  atm.  Electrode  erosion  on  helium  is 
negligible  in  the  anode  and  minor  in  the  cathode  over  a wide  test  range. 

A computer  code  based  on  a rod-arc  model  with  radiation  losses  and 
convective  turbulent  heating  of  the  gas  can  be  used  to  characterize  high 
pressure-high  mass  flow  rate  arc  heaters.  The  performance  of  an  N =*  0.5 
Huels-type  arc  heater  can  be  predicted  for  hydrogen,  helium  or  air. 

The  performance  of  air  arc  heaters  to  scale  N = 4 including  the  AEDC 
DET  (N  = 1. 73)  heater  can  be  accurately  characterized  by  this  code  thus 
verifying  its  scaling  credibility  and  giving  confidence  to  the  scaled 
predictions  for  hydrogen  helium. 
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APPENDIX  A:  DATA  SUMMARY  AND  COMPUTER  CODES 


A-1  DATA  SUMMARY 

Tables  A-1,  A- 2,  and  A-3  summarize  the  hydrogen,  helium  and  air  arc  heater 
data  from*  this  test  program.  The  data  shown  were  extracted  from  the  printed 
data  sheet  output  of  the  RDAS.  Each  point  represents  a stable  condition  in 
near  equilibrium.  A complete  volume  of  the  printed  data  sheets  was  furnished 
to  USAF-AEDC  under  separate  cover. 


The  column  shown  in  Tables  A-1  through  A-3  a^^e  as  follows. 


Symbol 

Parameter 

Units 

I 

Arc  current 

A 

P 

0 

Arc  pressure 

atm 

V 

Arc  voltage 

V 

m 

Gas  flow  rate 

Ib/s 

Anode  energy  loss 

Btu/s 

^C 

Cathode-body  energy  loss 

3tu/s 

Nozzle  energy  loss 

Btu/s 

Energy  balance  enthalpy 

Btu/lb 

Sonic  flow  enthalpy 

Btu/lb 

Note  that  the  sonic  flow  enthalpy  is  not  tabulated  for  hydrogen  since  it  was 
consistently  in  error  because  of  nozzle  throat-area  redaction  (see  Section  4.1) 

In  addition  to  the  data  given  in  Tables  A-1  through  A-3,  the  printed 
data  sheets  included  the  arc  power,  thermal  efficiency,  flowmeter  pressure, 
injector  pressure,  iniector  pressure  ratio,  differential  coolant  temperatures 
and  coolant  flow  rates  on  the  anode,  cathode-body,  and  nozzle.  The  data 
reduction  program  assigned  a value  of  1*F  to  the  differential  temperatures 
prior  to  and  following  the  test.  Thus,  since  the  specific  heat  of  water  is 
unity,  the  tabulated  energy  losses,  pre-  and  post- test,  were  numerically  the 
water  flow  rates  in  units  of  Ib/s. 


100 


AEOC-TR-76-2S 


TABLc  A-1  N-250  HYDROGEN  DATA 


Test 

I 

Po 

V 

m 

°A 

°N 

37 

237 

10.4 

876 

0.0042 

11.1 

54.3 

14.3 

29306 

296 

10.0 

741 

0.0041 

10.7 

56.3 

17.3 

32175 

400 

633 

0.0043 

14.7 

67.6 

17.2 

34410 

38 

236 

10.5 

1004 

10.0 

61.1 

14.4 

31151 

298 

9.8 

796 

12.0 

59.9 

15.7 

33312 

400 

10.0 

549 

0.0040 

15.9 

57.8 

15.2 

31706 

59 

240 

15.5 

1309 

0.0068 

19.7 

55.1 

24.3 

31107 

245 

15.9 

1300 

0.0068 

18.5 

55.6 

25.1 

31646 

314 

15.2 

992 

19.3 

67.4 

24.3 

34180 

323 

15.4 

970 

0.0057 

18.8 

67.6 

24.4 

34257 

398 

15.2 

757 

0.0048 

19.6 

76.6 

22.4 

36052 

60 

241 

15.1 

1419 

0.0065 

20.1 

61.8 

22.3 

35538 

242 

15.1 

1389 

0.0065 

19.4 

62.6 

23.3 

34525 

311 

14.5 

1109 

0 0052 

21.8 

78.3 

23.0 

38656 

392 

15.0 

939 

0.0052 

24.5 

95.3 

23.7 

41574 

56 

240 

25.0 

1984 

0.0106 

32.1 

96.2 

34.3 

313 

24.1 

1616 

0.0091 

32.5 

116.3 

33.9 

34413 

395 

24.9 

1421 

0.0086 

36.0 

149.2 

37.7 

38050 

57 

233 

25.2 

2124 

BEl 

35.7 

81.0 

34.6 

32647 

310 

24.8 

1727 

Bil 

34.2 

111.7 

39.8 

393 

22.6 

1229 

32.7 

128.4 

34.0 

37561 

64 

303 

49.3 

3442 

0.0194 

74.9 

266.0 

54.4 

32363 

74 

304 

49.8 

3578 

0.0208 

75.6 

307.1 

71.3 

30008 

63 

302 

10.1 

2924 

0.0179 

51.2 

250.3 

49.0 

29014 
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TABLE  A-2  N-250  HELIUM  DATA 


Test 

I 

Po 

V 

m 

Qc 

Qn 

hb 

hsf 

95 

240 

25.3 

968 

0.0203 

14.4 

54.0 

12.8 

7508 

7805 

318 

24.8 

871 

0.0194 

20.8 

72.5 

15.2 

8593 

8243 

396 

25.7 

807 

0.0186 

27.2 

91.6 

17.6 

9669 

9663 

106 

238 

25.2 

957 

0.0196 

13.1 

64.4 

14.2 

7012 

8332 

316 

25.3 

845 

0.0137 

19.4 

82.3 

16.5 

7891 

9234 

400 

25.3 

767 

0,0178 

25.8 

99.4 

17.9 

8938 

10204 

98 

237 

49.8 

1573 

0.0442 

34.8 

79.4 

19.7 

5615 

6392 

313 

50.0 

1473 

0.0409 

49.2 

107.4 

24.7 

7106 

7563 

397 

50.3 

1393 

0.0390 

62.0 

133.8 

29.4 

8356 

8395 

99 

237 

50.3 

1613 

0.0449 

32.9 

81.3 

19.4 

5762 

6351 

316 

49.9 

1518 

0,0421 

47.5 

110.8 

25.2 

7118 

7087 

399 

50.2 

1415 

0.0397 

61.1 

137.8 

29.5 

8392 

3088 

100 

242 

75.3 

2144 

0.0721 

52.4 

98.8 

21.1 

5092 

5513 

317 

75.1 

1997 

0.0668 

71.4 

138.8 

27.2 

6105 

6375 

397 

74.5 

1857 

0.0618 

92.1 

167.8 

34.6 

7204 

7335 

104 

236 

74.9 

2214 

0.0696 

49.8 

110.5 

24.3 

5115 

5839 

315 

75.0 

2059 

0.0642 

69.8 

148.3 

31.0 

6359 

6888 

398 

75.2 

1910 

0.0580 

90.3 

182.1 

36.0 

7747 

8479 

103 

238 

96.6 

2600 

0.0908 

61.8 

120.7 

25.0 

4834 

5706 

317 

94.0 

2415 

0.0804 

85.6 

163.1 

33.1 

6195 

6898 

399 

94.6 

2233 

0.0737 

109.9 

191.0 

39.5 

7514 

8304 

105 

238 

100.8 

2677 

0.0951 

66.2 

129.6 

27.3 

4663 

5665 

316 

100.9 

2556 

0.0870 

91.7 

171.8 

36.3 

6005 

6780 

396 

100.8 

2414 

0.0786 

119.3 

215.1 

45.8 

7361 

8308 
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TABLE  A-3  N-250  AIR  DATA 


T«st 

I 

Po 

V 

m 

Qc 

Qn 

•-b 

77 

238 

10.1 

0.0135 

7.0 

34.4 

5150 

5300 

7 ' 

315 

10.3 

0.0124 

9.5 

45.4 

10.3 

5760 

6730 

.7 

401 

10.1 

420 

0.0145 

11.4 

55.5 

9.8 

7230 

7700 

239 

9.9 

570 

0.0133 

7.3 

33.1 

8.0 

6140 

5130 

3C 

320 

10.1 

'480 

0.0132 

8.5 

41.5 

7.5 

6030 

5510 

BO 

9.9 

400 

0.0128 

9.7 

47.8 

7.4 

6800 

5660 

7-:> 

240 

15.1 

760 

0.0201 

15.C 

44.9 

10.0 

5?b0 

5290 

/5 

320 

14.8 

660 

0.0189 

19.6 

59.1 

11.8 

5910 

5910 

75 

393 

15.6 

550 

0.0178 

17.5 

67.7 

13.3 

76 

237 

15.2 

660 

0.0199 

10.9 

43.4 

11.5 

4240 

5510 

76 

317 

14.8 

540 

0.0171 

13.0 

53.4 

11.2 

5090 

7550 

76 

394 

15.1 

480 

0.0158 

15.5 

63.4 

10.9 

5850 

9680 

84 

239 

15.0 

730 

Bl 

34.9 

7.3 

5370 

4310 

81 

320 

15.0 

660 

waM 

55.1 

9.9 

6030 

5160 

82 

400 

14.8 

600 

19.3 

65.6 

1 1.2 

6800 

5360 

78 

242 

25.2 

1000 

0.0338 

15.2 

54  8 

18.0 

4340 

5160 

78 

322 

24.9 

830 

0.0299 

19.3 

70.7 

16.6 

5030 

6830 

78 

405 

25.0 

720 

0.0279 

22.1 

83.0 

16.2 

5650 

8260 

84 

238 

24.5 

0.0357 

15.8 

55.2 

13.5 

4280 

4220 

83 

320 

24.7 

0.0336 

22.9 

70.1 

15.2 

5100 

84 

850 

0.0326 

33.3 

86.7 

16.6 

5830 

5580 

92 

238 

49.5 

1662 

0.0691 

27.4 

77.8 

18.2 

3770 

4680 

92 

317 

50.4 

1474 

0.0600 

37.8 

21.0 

4440 

5520 

92 

400 

50.9 

1343 

0.0630 

485 

124  9 

23.9 

5070 

6340 

107 

236 

50.5 

1663 

0.0721 

28.5 

76.7 

19  7 

3560 

107 

320 

51.2 

1453 

0.0683 

36.9 

22.4 

4250 

5260 

107 

400 

50.9 

1307 

0.0652 

42.1 

123.6 

25.0 

4800 

5820 

93 

239 

75.8 

2263 

0.1059 

38.0 

92.8 

25.0 

3500 

4700 

93 

320 

75.2 

1976 

50.9 

122.2 

28.7 

4090 

5210 

93 

393 

75.1 

1804 

0.0966 

62.6 

149.9 

32.5 

4550 

5780 

94 

317 

73.8 

1973 

0.1017 

49.6 

118.8 

26.7 

4040 

4850 

108 

234 

3008 

0.1444 

47.7 

111.2 

32.2 

3430 

4330 

108 

319 

2462 

0.1355 

60.1 

142.2 

33.7 

3880 

5160 

108 

398 

100.2 

2260 

0.1313 

72.8 

179.2 

40.6 

4400 

5510 

109 

314 

100.1 

2450 

0.1311 

55.4 

140.7 

33.1 

3950 

5510 
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A-2  HYARC  AND  HEARC  COMPUTER  CODES 

The  computer  codes  HYARC  ani  HEARC  are  essentially  identical  with  the 
exception  of  those  parameters  ?aat  characterize  each  individual  gas.  Thus, 
to  avoid  unnecessary  duplication,  only  a complete  description  of  the  HEARC 
is  given;  those  aspects  which  differ  for  HYARC  are  discussed  in  the  following 
subsection. 

HEARC  Program 

A description  of  the  function  of  the  main  program  and  each  subroutine 
is  given  in  Table  A~4.  The  computer  variable  descriptions  are  listed  in 
Tables  A- 5 and  A- 6. 

The  FORTRAN  HEARC  program  listing  is  presented  in  Table  A-7  and  the 
corresponding  flow  chart  is  presented  in  Figure  A-1. 

HYARC  Program 

A description  of  the  function  of  the  main  program  and  each  subroutine 
is  given  in  Table  A-4.  The  differences  between  HEARC  and  HYARC  are  those 

TABLE  A4  HYARC-HEARC  PROGRAM  FUNCTIONS 


Main  Program 


Subroutines: 

INPUT 

GASPOW 

ARCTEM 

ARCDIA 

OUTPUT 

GAMMAX 

TEMP 

VISCX 

ENTHX 

RHO 

SIGMAX 


RADX 


TCONDX 


Sets  up  pressure  and  current  calculational  loops  and  initial  values  for  TA*RC,  TGAS,  and  DARC. 
Iterates  gas  temperature  until  gas  power  from  enthalpy-mass  flow  relation  agrees  with  energy 
balance  value.  Iterates  mass  flow  rate  until  it  agrees  with  value  calculat  ?d  for  sonic  flow  at 
nozzle  throat. 


Redds  in  required  data  in  engineering  units  and  converts  to  MKS  system 

Calculates  gas  power  uy:  (1)  energy  balance  on  the  qas  and  (2)  enthalpy-mass  flow  of  arc 
column  and  bulk  gas  annulus 

Calculates  temperature  of  the  arc  column  for  minimum  power  input 
Calculates  arc  column  diameter  bv'  solving  the  arc  column  power  balance  equation 

Computes  voltage,  gradient,  efficiency,  and  electrode  heat  loads  and  presents  results  of 
solution  in  engineering  units 

Finds  ratio  of  specific  heats  for  given  temp'jrature  and  pressure 
Finds  temperature  that  corresponds  to  a given  enthalpy  and  pressure 
Finds  viscosity  from  table  interpolation  for  a given  temperature  z nd  pressure 
Finds  enthalpy  from  table  interpolation  for  a given  temperature  and  pressure 
Finds  gas  density  from  table  interpolation  for  a given  temperature  and  pressure 

Helium.  Finds  electrical  conductivity  from  table  interpolation  for  a given  temperature  and 
pressure 

Hydrogen:  Finds  electrical  conductivity  by  a curve  fit  for  a given  temperature  and  pressure 

Helium:  Finds  volumetric  radiation  intensity  from  electron  mole  fraction  table  interpolation 
and  electron  density  calculation 

Hydrogen:  Obtains  volumetric  radiation  intensity  from  table  interpolation  for  a given 
temperature  and  pressure 

Finds  thermal  conductivity  from  table  interpolation  for  a given  temperature  and  pressure 
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' parameters  that  are  characteristic  for  each  gas.  The  thermodynamic  and 

^ transport  properties  are  from  the  sources  shown  in  Table  4.  The  descriptions 

I of  the  HYARC  main  program  and  suBro^*  tines  tnPUT,  OUTPUT,  GASPOW,  ARCTEM, 

' ARCDIA,  and  TEMP  are  shown  in  Table  A-4.  The  thermophysical  property 

determinations  follow  the  saiae  approach  as  HEARC.  The  front  and  rear  arc 
lengths  are  shown  in  Table  7 for  the  N-250.  The  cold  gas  inlet  enthalpy 
and  Prandtl  number  are  listed  in  Table  5.  The  correlation  results  utilized 
are  discussed  in  Section  5.4. 


TABLE  A-5  HYARC-HEARC  VARIABLE  DESCRIPTIONS 


'/triable  Name 

Description 

Variables  common  to  more  than  one  subroutine 


AMPS 

Aic  current,  A 

AMPINC 

Current  increment,  A 

ARCLEN 

Total  length  of  arc,  m 

I 

CNVF01 

39.3701  in./m 

t 

CNVF02 

2.2046  Ibm/kg 

1 

CNVF03 

0.0254  m/in. 

\ 

CNVF04 

4.2292  X 10““^  (Btu/lbm|/(J/kg) 

CNVF05 

2.6839  X 10~®  (Btu/ft^)/(J/m^) 

i 

CNVF06 

8.8056  X 10~®  (Btu/ft^s)/(W/m^) 

CNVF07 

0.2048  (Ibm/ft^sj/lkg/m^s) 

1 

CNYF08 

3.2808  (ft/s)/(m/s) 

i 

CONSLN 

Length  of  front  constrictor,  m 

> 

1 

DARC 

Dtameter  of  arc,  m 

* 

DCONS 

Diameter  of  constrictor,  m 

DTHRT 

Diameter  of  nozzle  throat,  m 

I 

ENTH 

Net  gas  exit  enthalpy,  J/kg 

i 

1 

FARCLN 

Arc  length  in  front  electrode,  m 

i 

1 

FNU 

Nusselt  number 

i 

lARCD 

Arc  diameter  iteration  counter 

lARCT 

Arc  temperature  iteration  counter 

lERR 

flag  indicating  that  a preset  iteration  count  has  been  exceeded 

ISONIC 

Mass  flow  rate  iteration  counter 

ITGAS 

Bulk  gas  iteration  counter 

.NAMPS 

Number  of  current  increments 

NCASE 

Index  indicating  first  call  to  output  subroutine  for  each  input  case 

NPRES 

Number  o^  pressure  increments 

PCONV 

Power  convected  from  gas  to  constrictor  wall,  VV 

PGAS 

Power  in  gas  from  energy  balance,  W 

PGASEN 

Power  in  gas  from  enthalpy-mass  flow,  W 

PI 

3.14159 
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TABLE  A-5  (CONTINUED) 


Variable  Name 

Description 

PIN 

Electrical  power  into  arc,  W 

HINC 

Pressure  increment,  atm 

PRES 

Pressure,  atm 

QCONV 

Constrictor  wall  convective  heat  flux,  W/m 

QN 

Nozzle  energy  loss,  W 

QRAD 

Constrictor  wall  radiative  heat  flux,  W/m^ 

RARCLN 

Arc  length  in  rear  electrode,  m 

REV 

Reynolds  number 

SIGMA 

Arc  electrical  conductivity,  mhos/m 

TARC 

Arc  temperature,  K 

TGAS 

Bulk  gas  temperature,  K 

TSTAG 

Stagnation  temperature  of  exit  gas,  K 

WDOT 

Mass  flow  rate,  kg/s 

WDOTAR 

Mass  flow  rate  in  arc  column,  kg/s 

WDOTGA 

Mass  flow  rate  in  oulk  gas  annulus,  kg/s 

Y 

o 

Volumetric  radiation  of  arc  x (7T/4),  W/m 

Variables  used  in  MAIN 

ASTAR 

9 

Area  of  nozzle  throat,  m 

DELN 

DELPGA  value  from  previous  iteration 

DELPGA 

Normalized  gas  power  difference 

DELWDT 

Normalized  r iss  ficv/  rate  difference 

GAMMA 

Ratio  of  gas  specific  heat,  Cq/Cw 
, 2 , RATIO 

GR 

NEXT 

Index  for  new  case  or  end 

RATIO 

(7+1)/(2(7'-D) 

RHOSTG 

Gas  density  corresponding  to  stagnation  temperature,  kg/m 

SCALE 

Scale  N of  heater,  ratio  of  constriction  diameter  to  0.^5  in. 

TABS 

Absolute  value  of  difference  between  present  and  previous  bulk  gas  temperatures,  K 

TGASN 

Bulk  ges  temperature,  from  previous  iteration,  K 

WDOTSF 

Mass  flow  rate  calculated  from  nozzle  sonic  flow  relationship,  kg/s 

XP 

New  bulk  gas  temperature  by  method  of  false  position,  K 

Variables  used  in  subroutine  GASPOW 

AARC 

o 

Arc  column  area,  m 

AGAS 

Bulk  gas  annulus  area,  m^ 

ENTHAR 

Arc  gas  enthalpy,  J/kg 

ENTHGA 

Bulk  gas  enthalpy,  J/kg 

1C6 


AEDC-TR-76-25 


TABLE  A-5  {CONTINUED) 


Variable  Name 

Description 

H 

Turbulent  heal  transfer  coefficient.  W/m^K 

HIN 

Cold  gas  inlet  enthalpy.  J/Kg 

PR 

Prandtl  number 

REY 

Reynolds  number 

RHOARC 

'3 

Arc  gas  density,  kg/m 

R HOG AS 

Bulk  gas  density,  kg/m 

TA 

Dummy  temperature  variable.  K 

TCOND 

Thermal  conductivity.  W/m  • K 

TWALL 

Constrictor  wall  temperature.  K 

VEL 

Gas  velocity  in  constrictor,  m/s 

Variables  used  in  ARCTEM 

DIRCH 

Equals  minus  one  if  PiN  greater  than  PINP.  otherwise  equals  plus  one 

DT 

Arc  temperature  increment.  K 

PINP 

Previous  calculated  value  of  power  input.  W 

SIGNF 

Sign  of  arc  temperature  incremen 

Variables  used  in  ARCDIA 

ALPHA 

2/Y,  used  in  solving  energy  balance 

BETA 

X/Y,  used  in  solving  energy  balance 

CORR 

Term  that  takes  into  account  fractional  variation  of  DARC 

DELDIA 

Normalized  arc  diameter  difference 

DELD 

Arc  diameter  iteration  increment,  m 

D1 

Previous  iteration  value  of  DARC,  m 

DP 

Used  in  CORR,  D^,  + Dg 

DM 

Used  in  CORR,  - Dg^ 

F 

Arc  column  energy  balance  function,  P,  -(Pr+pQ)  = 0 

FNU 

Nusselt  number 

FP 

First  derivitive  of  F.  used  for  arc  diameter  iteration 

H 

Turbulent  film  heat  transfer  coefficient.  W/m  K 

INN 

Subroutine  entry  indicator  flag 

PR 

Prandtl  number 

PXX 

Dummy  variable  for  power  input,  W 

RAD 

Volumetric  radiation,  W/m^ 

TCOND 

Thermal  conductivity,  Vy//m  • K 

TXX 

Dummy  variable  for  arc  temperature.  K 

TX 

Dummy  temperature  for  evaluating  TCONDX.  K 

VSC 

Viscosity  of  bulk  gas,  N s/m^ 

X 

Power  input  per  unit  length  - D^,  , W/m 

2 

Convective  heat  flux  loss  from  arc  per  unit  length  t Dg,  W/m^ 
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TABLE  A-5  (CONCLUDED) 


Variable  Name 

Description 

Variables  used  in  subroutine  OUTPUT 

DCONSX 

Constrictor  column  diameter,  in. 

DIARCX 

Arc  column  diameter,  in. 

DTHX 

Nozzle  throat  diameter,  in. 

EFF 

Efficiency 

ENTHX 

Exit  gas  bulk  enthalpy,  Btu/lbm 

FRNT 

Front  arc  length,  in. 

GRAD 

Arc  column  voltage  gradient,  V/in. 

HFRNTX 

Front  electrode  (constrictor)  power  loss,  Btu/s 

HREARX 

Rear  electrode  power  loss,  Blu/s 

INN 

Subroutine  entry  flag 

lOUT 

Logical  device  for  output 

ICONS 

Constrictor  length,  in. 

PGASX 

Net  power  added  to  gas,  kW 

PINX 

Power  input,  kW 

QCONVX 

Constrictor  wall  convective  heat  flux,  Btu/ft  s 

QNX 

Nozzle  heat  load,  Btu/s 

QRADX 

2 

Constrictor  wall  radiative  heat  flux,  Btu/tt  s 

REAR 

Rear  arc  length,  in. 

VOLTS 

Arc  column  voltage  drop,  V 

WDOTX 

Tjtai  mass  flow  rate,  Ibm/s 

Variables  used  in  subroutine  TEMP 

CONV 

Normalized  temperature  convergence  parameter 

DT 

Temperature  increment,  K 

ERR 

Convergence  error 

H 

Enthalpy,  J/kg 

HTN 

Enthalpy  corresponding  to  temperature  TN,  J/kg 

HTNP 

Enthalpy  corresponding  to  temperature  TN  + DT,  J/kg 

K 

Iteration  index 

P 

Pressure,  atm 

T 

Temperature,  K 

Tl 

Initial  starting  temperature,  K 

TN 

Previous  iteration  temperature  value,  K 

TNP1 

Temperature  calculated  for  this  iteration,  K 
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TABLE  A-6  HEARC  GAS  PROPERTY  VARIABLE  DESCRIPTIONS 


Variable  Name  Description 

Variables  used  in  subroutine  VISCX 


ATM!  Data  temperature  array  for  1 atmosphere 

ATM  10  Data  temperature  array  for  10  array  for  10  atmospheres 

ATM  100  Data  temperature  array  for  100  atmospheres 

ATMIOOO  Data  temperature  array  for  1000  atmospheres 

DT  Temperature  difference  between  table  value  J and  T 

I Pressure  index  for  table  values 

J Temperature  index  for  last  table  value  before  T 

P Pressure,  atm 

PR  Pressure  values  in  table  array,  atm 

PP  Logarithmic  pressure  increment 

T Temperature,  K 

2 

TC  Viscosity  table  array,  dyne  — s/cm 

TCI  Viscosity  interpolated  to  temperature  T at  pressure  PR  (I) 

TC2  Viscosity  interpolated  to  temperature  T at  pressure  PR  (I  + 1) 

TVISC  Viscosity,  Nt  s/m^ 

Variables  used  in  subroutine  ENTH 


ATM! 

ATM  10 

ATMIOO 

ATMIOOO 

DT 

I 

J 

P 

PR 

PP 

RHE 

TENTH 

T 

TC 

TCI 

TC2 

TO 


Data  temperature  array  for  1 atmosphere 
Data  temperature  array  for  10  atmospheres 
Data  temperature  array  for  100  atmospheres 
Data  temperature  array  for  1000  atmospheres 
Temperature  difference  between  table  value  J and  T 
Pressure  index  for  table  values 
Temperature  index  for  last  table  value  before  T 
Pressure,  atm 

Pressure  values  in  table  array  TC,  atm 
Logarithmic  pressure  increment 
Helium  gas  constant,  J/kg  K 
Enthalpy  corresponding  to  T,  P,  J/kg 
Temperature,  K 

Enthalpy  table  array,  normalized  to  Rue^o 
Enthalpy  interpolated  to  temperature  T at  pressure  PR  (I) 
Enthalpy  interpolated  to  temperature  T at  pressure  PR  (I  +•  1) 
Standard  temperature,  273.16  K 
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TABLE  A-6  (CONTINUED) 


Variable  Name 

Description 

Variables  used  in  subroutine  BHO 

ATMl  Data  temperature  array  for  1 atmosphere 

ATM  10  Data  temperature  array  for  10  atmospheres 

ATM  100  Data  temperature  array  for  100  atmospheres 

ATM  1000  Data  temperature  array  for  1000  atmospheres 

DT  Temperature  difference  between  table  value  J and  T 

I Pressure  index  for  table  values 

J Temperature  index  for  last  table  value  before  T 

P Pressure,  atm 

PP  Logarithmic  pressure  increment 

PR  Pressure  values  in  table  array  TD.  atm 

RHOO  0.17859,  density  normalization,  kg/m^ 

T Temperature,  K 

TO  Density  table  array  normalized  to  RHCO 

TOENS  Density  corresponding  to  T and  P,  kg/m 

TD1  Density  interpolated  to  temperature  T at  pressure  PR  (I) 

T02  Density  interpolated  to  temperature  T at  pressure  PR  (1+  1) 

TO  Standard  temperature,  273.16  K 


Variables  used  in  subroutine  SIGMAX 


ATMl 

ATMIO 

ATMIOO 

ATMIOOO 

DT 

I 

J 

JSHIFT 

P 

PR 

PP 

T 

TS 

TSIG 

TS1 

TS2 


Data  temperature  array  for  1 atmosphere 
Data  temperature  array  for  10  atmospheres 
Data  temperature  array  for  100  atmospheres 
Data  temperature  array  for  1000  atmospheres 
Temperature  difference  between  table  value  J and  T 
Pressure  index  for  table  values 
Temperature  index  for  last  table  value  before  T 
Temperature  index  shift  parameter 
Pressure,  atm 

Pressure  values  in  table  array  TS,  atm 
Logarithmic  pressure  increment 
Temperature,  K 

Electrical  conductivity  array,  Oog^Q  to  13,000  K)  mho/rm 
Electrical  conductivity  corresponding  to  T and  P,  mho/m 
Electrical  conductivity  interpolated  to  temperature  T at  pressure  PR  (I) 
Electrical  conductivity  interpolated  to  temperature  T at  pressure  PR  (I  + 1) 
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TABLE  A-6  (CONCLUDED) 


Variable  Name  Description 

used  in  subroutine  RADX 


ATM!  Data  temperature  array  for  1 atmosphere 

ATMIO  Data  temperature  array  for  10  atmospheres 

ATMIOO  Data  temperature  array  for  100  atmospheres 

ATMIOOO  Data  temperature  array  for  1000  atmospheres 

DT  Temperature  difference  between  table  value  J and  T 

1 Pressure  index  for  table  values 

J Temperature  index  for  last  table  value  before  T 

JSHIFT  Temperature  index  shift  parameter 

KP  Helium  equilibrium  constant 

KPP  Square  root  of  KP  -r  Pressure 

o 

NE  Electron  number  density,  -/cm 

NT  Total  particle  number  density,  -/cm 

P Pressure,  atm 

PP  Logarithmic  pressure  increment 

PR  Pressure  values  in  table  array  TR,  atm 

T Temperature,  K 

TR  Electron  mole  fraction  array  (logic  lo  16,000  K) 

TRl  Mole  fraction  interpolated  to  temperature  T at  pressure  PR 

TR2  Mole  fraction  • nerpolated  to  temperature  T at  pressure  PR 

TRAD  Volumetric  radiation  intensity  at  T and  P,  W/m^ 

XE  Electron  mole  fraction 

Variables  used  in  subroutine  TCONDX 


(I) 

(1+1) 


ATM1 

ATMIO 

ATMIOO 

ATMIOOO 

DT 

I 

J 

P 

PP 

PR 

T 

TC 

TCI 

TC2 

TCOND 


Data  temperature  array  for  1 atmosphere 
Data  temperature  array  for  10  atmospheres 
Data  temperature  array  for  100  atmospheres 
Data  temperature  array  for  1000  atmospheres 
Temperature  difference  between  table  value  J and  T 
Pressure  index  foi  table  values 
Temperature  index  for  last  table  value  before  T 
Pressure,  atm 

Logarithmic  pressure  increment 
Pressure  values  in  the  table  array  TC,  atm 
Temperature,  K 

Thermal  conductivity  array,  erg/cm  s K 

Thermal  conductivity  interpolated  to  temperature  T at  pressure  PR  (I) 

V 

Thermal  conductivity  mterpolated  to  temperature  T at  pressure  PR  (I  + 1) 
Thermal  conductivity  for  T and  P,  W/m  K 
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TABLE  A-7  HEARC  LISTING 


C MINiNUP  ENCRCy  ADDITION  APC  HEaT6».  CKAQ  ACT  E P I / AT  I ON  ALCOBITHH 
C FOR  HELIUM 

I MCDONNELL  DOUGLAS  RESEARCH  LABORATOoIE S,  ST. LOUIS,  mq. 
i PRUCRAMtD  BY  J.F.SHAEFPcP 
C 

BLOCK  data 
COMMON/  Cf.NVF  / 

♦ CNVFOl,  CNVF02.  CNVF03,  CNVPO<., 

♦CNVF05,  CNVF06,  CNVF07,  CNVF08,  PI 

r 

DATA  CNVFOl, CNVF 02. CNYF03.CNVF0A .CNVF05 .CNVF06 , C V/F07, CNV FOB , P I 
♦/39.3701,  2-20<.fe,  .0254,  a. 29926-4,  2.68396-5,  8.80566-5, 

♦0.2043,  3.2308,  3.14159  / 

END 

C 

: CNVFOl*  M • 39.3701  • INCHES 

: CNVF02=  KG/ScC  • 2.204©  « lBM/SEC 

: CNVFOi*  VOLTS/M  ♦ 2.546-2  » VOLTS/lNCM 

: CNVF04*  J0UL6S/KG  • 4.29‘>?6-4  - BTU/L8M 

C CNVFOS*  JOULES/M  3 • 2-68396-5  flTU/FT  3 

: CNVF06*  HATTS/M  2 * 8.8056E-5  ■ BTU/FT  2-SEC 

: CNVF07-  KG/M  2-SEC  ♦ 0.2046  » LHM/FT  2-SEC 

C CNVF08*  M/SfcC  ♦ 3.2808  » FT/ScC 
COMMON  / CONVF  / 

♦CNVFOl,  CNVF02,  CNVF03,  CNVFQ*., 

♦ CNVF05,  CNVF06,  CNVF07,  CNVF08,  P! 

COMMON/  BLKl  / AHPS.ARCLEN, OARC ,0C3NS, OTHPT , ENTH, PC.^iV , PCA< , 
♦PGA  SEN,  PIN,PRA0,RRE  S,QC0NV,URA0, SIGMA, T ARC, T GAS, woof, Y 
COMMON/  BLK2  / ON.tSTAC 

COMMON/  ERROR  / I ARCO, UftCT , I ERR, I SONIC , I TCAS 

COMMON/  LOOPS  / AMPINC ,NAHPS,PI nC.NPRES 

COMMON/  LENGTH  / C ONSLN,FAPC LN , RARCl N 

COMMON/  NCASfc  / NCAS6 

COMMON/  REY  / R E Y.WOOTAR .WOOTGA ,FNU 

IN«5 

TAKC*18000. 

ICAS»2900.0 
100  CALL  INPUT 

hdotga*mdgt 

SCAL  6*1  DC0NS/(0.  75*0. 025*.  n 
NCASE-O 

ampst*ahps 

PREST.PRES 

0ARC*0.9*0CUNS 

DU  160  1P*1,NPRES 

OU  160  IA»1,NAMPS 

AHPS-AMPST  ♦ AMPINC*(IA-1I 

PR6S*PR6ST  ♦ PINC*tIP-l) 

r.  ARC  LENGTH  CORRELATION  FOR  HELIUM  N-250  HEATER 

RARCLN«IPRES/(0.288*PR6S  ♦ 4.22 ) ) ♦CNVF03 
FARCLN*(4.45  ♦ 0.0106* IPPES-26. 1-RARCLN/CNVF03I *CNVF03 
FARCLN*( SCALE/0.5I*FARCLN 
RAPCL.N*!  SCALE /O.  S)*RARCLN 
ARCLEN-FARCLN  ♦ RARCLN 
IERR*0 
I SON  1C *0 
no  1TGAS*0 
120  CONTINUE 

: CALCULATE  ARC  TEMPERATURE  AND  DIAMETER 

CALL  ARCTEM 
IF(  lERR.EO.UGO  TO  150 
call  GASPOM 
ITCAS*! TGAS+l 
IF!  ITGAS.GT.lOO)  16RR*I 
IF!  IcRR.EO.l)  GO  TO  150 
delpga  *|PGAS-PGASEN»/PGAS 
IF(  ITGAS.NE.U  GO  TO  130 

oeln*delpga 
TGASN*TGAS 
TCAS-rCAS  - 100.0 
GO  TO  120 

130  XP*nGASN*DELPGA  - TGA  S*OEL N)  / ( DEI PGA-OELN) 

TA8S»AbS(XP  - .'GAS) 

IFl  TABS. GT, 300.  ) XP«TGAS  ♦ ( XP  - TC  A S)  / TABS  *100  , 

TCASNaTCAS 

CELN*06LPGA 

IFIABSIOELPC  ).LE..On  GO  Tr  \uq 
TGAS-XP 
GO  TO  120 
140  CONTINUE 

C CALCULATE  HASS  FLOM  T()  SAflS'Y  SONIC  FLOW  RELATION 
ENTH*PGAS/WD0T 
CALL  TcHP(ENTH.PP[ S.TSTAG) 

CALL  GAHMAXI TSTAG.PRES.GAMMAT 
CALL  RHOl TSTaG.PRE S.RMOSTG) 
kATI0«(GAMMA^1.0)/i 2.0* (GAMMA-1.0) ) 

GP «( 2, 0/ (GAMMAS  1.0) )**RATIO 
ASTAR«P!*OTHRT*DTHRT/A.o 

WOOTSF  -ASTAP*SQRT(CAHMA*RH(JSTG*PRES*1.013E5|*GR 

IS0NIC*IS0NIC*1 

IF!  ISONIC.GT.lOO)  I£RR*l 

1F(  lERR.EO. 1)  GU  TO  150 

DfcLHOT*(MDOTSF  - WOOT)  / wOOT$F 

IF(AbSU)ElWOT).LT.O.Ol)  GO  TO  1^0 

WDDT*wOOTSF 

GO  TO  no 

150  CONTINUE 

CALL  OUTPUT 

160  continue 

PEAOnh.175)  NFXT 
175  F0RMAT(9X, I 1) 

CO  TO  ( 100.  1^0)  ,NE  XT 
170  continue 
END 


subroutine  input 
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TABLE  A-7  (CONTINUED) 


CARD  NO 


C'INTENTS 


109 

no 

ni 

112 

113 

!U 

H? 

lie 

m 

123 

12^ 

125 

126 
127 
120 
120 

130 

131 

132 

133 
136 

135 

136 

137 

138 

no 

160 

161 

162 

163 

166 

165 

166 
167 
160 
169 
150 

pi 

152 

153 
156 

155 

156 

157 

158 

159 

160 
161 
162 
163 
166 

165 

166 
167 
166 

169 

170 

Hi 

176 

175 

176 

177 

178 

179 

180 
181 
182 
103 
106 
105 
186 
187 
108 
169 

190 

191 

192 

193 
196 

195 

196 

197 

198 

199 

200 
201 
202 
203 
206 

205 

206 

207 

208 

209 

210 


i 


C READS  REQUIRES  DATA  IN  ENGINEERING  UNITS  AND  CONVERTS  TO 
C KKS  SYSTEK  FOR  IMF  CChPUTAIIONS 
CO^'J'ON/  CCNVF  / 

♦ CNVh01,CNVF02»CNVF03tCNVF06  ,CNVCQ5  ,CWC06,CNVF07  ,CNVF03t  PI 
CCHSLN/  BLKl  / AMP S , ARC DARC  t DC0NS» OTHAT  , ENTH.  PCCl4Vt  ^CAS i 

♦ PCASEN,  PIN.PRAD,PRBS,OCONV,CRAOfSIGMA,TARC|TGAS.WOOT  »Y 
COMMON/  imps  / AMPINC.NAHPSiPINCiNPRES 

COMMON/  TITLE  / 7ITLEii6).OAT(6| 

COMMON/  LENGTH  /CONSLN,  FARCLN,  RARCLN 
IN-5 

REA0(lN,100MTITLE(J)tJ-ltld) 

100  50RHAT( l8A6) 

KLAOI IN. llOlAHPSiAMPlNC.NAMPS 
no  FORMAT!  2E10.1.6X,  12) 

READ!  IN,nO)PRES|PlNC,NPRES 
READ! IN. l20)CUNSLNtFARCLN.RARCLN 
120  F0RHATI3E10.2) 

REAOl IN. 120)WDOT,DCONSiOTMRT 
C0NSLN*C0NSLN*CNVF03 
FARCLN -FARCLN»CNVf 03 
RARCLN»RARCLN*CNVF03 
ARCLEN-FAPCLN  * RARCLN 
0R6AR«0REAR*CNVF03 
wr0T-.<00T/CHVF02 
0C0NS-DC0NS6CNVF 03 
OTHRT-OTHR  T*CNVFQ3 
CALL  DATIHE(OAT) 

RETURN 
END 

SUBROUTINE  GASPOR 

COMPUTES  POWER  IN  CAS  BY  CALCULATING  CONVECTIVE  LOSS  TO  CONSTRICTOR 

CUHHUN/  BLKl  / AHPS.ARCLEN.OARC.OCONS.OTHRT.ENTH.PCONV.PCAS. 
♦PGASEN,  PIN.PRAOfPREStQCONVtQRAD.SlGMA.TARC«TGAS.MOOT,Y 
COMMON/  BLK2  /ON.TSTAG 

COMMON/  ERROR  / lARCO, lARCi .If RP.ISCNIC. ITCAS 
COMMON/  lENCTM  /CONSLN . FARCLm,  RAKCLN 
COMMtyj/  REY  / REY,WOOTAR,W>:dTCA,FNU 
PI-3.16159 
TWALL-590. 

: CALCULATE  CONVECTION  LOSS  TO  CONSTRICTOR  HALL 

TA-l.OO*! TGAS  - TWALU  ♦ THALL 
CALL  TCONDU  T A ,PRE  S.  TCCNO  I 
M-FNU*TCONO/  U)C6nS-0ARC) 

QCONV-H*! TGAS-TWALL) 

L RADIATED  HEAT  FLUX  OUT  FRCM  ARC 
ORAD-Y*(OARC*OARC) /(PI  *OCONS) 

: POWER  IN  CAS  FROM  ARC  (BY  CONVECTION)  IS  Z *DARC* ARCL EN 

: HIN-COLO  CAS  ENTHALPY  (JOULES/KG) 

MIN-670. 0/^.2992E“6 

: LOSSES  RADIATION  OVER  ARCLFNGTM.  CCNV.  OVtR  FRONT  CONSTRICTOR 

P&AS-PIN  - QRAO*PI*OCONS*ARCLlN  - OCONV*P!*DCONS*CONSLN 
44  MOOT*H1N 

: fWZZLE  LOSS  COPRElATICN  FOR  HELIUM, CN  IS  LOSS  IN  WATTS 

6NTM-(PGAS/WDnT)*6.3E-6 

ON-0.06*(INOOT*2.2066)**0.7)*(ENTM-1325.) *105^. 

PCAS-PGAS  - QN 
C 

C CALCULATE  POWER  IN  CAS, ARC  ^ ROM  ENTHALPY-HASS  FcOH  RELATION 
C ASSUME  UNIFORM  VELOCITY 

CALL  RHO(YGAS,PRES,RHOGAS) 

CALL  RH0(TARC,PRES»RH0ARC) 

AGAS-PI*CDCONS**2-0ARC**2)/6.0 
AARC-P I *OAPC*OARC /6.0 
VEL-WDOT/(RHOGA$*AGAS  4 RNOARC^AARC) 

WDOTAR -RhOARC*VEL*AARC 
W00TCA«RH0GA$*V£L*ACAS 
CALL  ENTMX(TCAS,PRES,ENTHGA) 

CALL  ENThXITARC.PRES.ENTHAR) 

PCAS£N»ENTMGA*W00TCA  4 ENTHAR*WOOTAR 
RETURN 
END 
C 

c 

subroutine  arc  TEH 

C 

C CALCULATES  ARC  TEMPERATURE  F CR  POWER  INPUT  TO  BE  A MINIMUM 
C 

CUMMUN/  BLKl  / AMP S . APCl 6 N , OARC .OCONS r OTHRT , ENT H, PCONV , PGAS , 
♦PGASFN.PIU.PPAD.PRE  S.OCaNV.ORAO.SIOMA.TARC.fCAS ,HOOT ,Y 
COMMON/  ERROR  / URCO,  I ARC  T .1 1:  RR,  I SONi  C > 1 TCAS 
SICNF-l.O 
lARCT-O 
OT-200.0 

100  IARCT-IARCT41 

CALL  ARCOfAlTARC  .PIN  I 

)F(IERP,f0.l)RETURN 

IF( IARCT.NE.11G0  TO  50 

TAHC-TARC  ♦ OT 

PINP-P IN 

CO  TO  100 

50  DIR  CM- SI  GM  1.0,(PINP-PINM 
SIGNF-UIRCH-SIGNF 
DT«nT/( !,♦( l.-riRCM)/2.I 
1F(AHS(DT  I.LT.IO.O)  return 

TARC-TARC  ♦ SICNF*0T 
PINP-P  IN 

IF(  lARCT.GT.lOO)  GO  TO  UO 
CU  TO  100 


jJP  JERR’l 


0 return 

END 


SUBROUTINE  ARCOIAITXX.PXXi 
CALCULATE  A^C  DIAMETER  FOR  FNERCY  BALANCE 
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TABLE  A-7  (CONTINUED) 


COMM(>J/  BlKI  / AMPS  iAPCLEN,DA«C,OCCNS,UrMRT, ENT H.PCCNV.P&AS. 
♦PGASEN,  PIN-, P^AOfPRESfUCONVtORAO, SIGMA, TAPCiTGAS.WOOT  ,Y 
CUPMON/  ERROR  / lARCn, lARCT,tERR,lSCNIC,ITCAS 
COMMON/  REY  / REY,WOOIAR,MOOTGA,FNU 
DATA  IKN/0/ 

IF(  INN.NE.O)  GO  TO  100 
PI-3.l<.159 
01«0ARC 
100  INN«1 
IARCO-0 

POWER  IN  » X(TARC)/I)ARC**2  per  UNIT  LENGTH 
CALL  SIGHAXITXX  ,PRES,SIGHAI 
X«A.O»AHPS*AMPS/(PI*SIGMA| 

CALL  RAOXITXX  ,PRES,RaO) 

POWER  RAOUTFO  ■ Y(T  ARC ) •0ARC**2  PER  UNIT  LENGTH 
Y»P I*RAO/A.O 

TX«0.70*(TAPf  - TCASI  ♦ TCAS 
CALL  VISCXiTCASfPRES.VISO 
CALL  TCONDXITX.PRES.TCGNO) 
llO  CONTINUE 

lARCO-  lARCO-*^! 

PXX»(X/0APC**2I*ARCLEN 

POWER  OUT  Of  ARC  COLUMN  CONVECTEO  « 2(TARC1*CAVC  PER  UNIT  LENGTH 
REY«IA.O*WOOTCA) /|P| ♦VISC*(OARC*DCCNSI I 
PRANOTL  NO  FOR  HE  IS  0.672 
PR-0.672 

FNU-0.0265*(RE Y**0.6)*PR**0.3 
H-FNU*  TCOND/IOCCNS-OARC) 

2*PI*H*(TXX  -TCASI 

ALPMA-2/V 

BETA— l.*X/Y 

F«(0I**4)  ♦ ALPHA*|0I**3I  ♦ BETA 
OP-CCONS  ♦ 01 
OM«OCONS*OCGNS  - 01*01 

CORR-1.0  ♦ 0. 2*01/3. 0/OP  ♦ ?.0*f)I  *DI /3  .O/OM 
FP-4.0*I0I**3)  ♦ 3.0*ALPHA*OI*OI*CORR 
OELO— l.*F/FP 
CARC-OI  ♦ OELO 
OELOIA»(UARC-DII/OARC 
1F(  1ARCD.GT.100)  GO  TO  120 
|fUBSI0ELt)lA),LE.l.E-5)  CO  TO  130 
Di -OARC 
GO  TO  110 
120  lERR-l 
RETURN 
130  continue 


<^UBRCUTINE  OUTPUT 

: lALCULATE  REMAINING  ARC  PARA^JETERS 

: CONVERT  MkS  to  ENGINEERING  UNITS  A-.O  PRINT  CUT  RESULTS 

I EAL  LENX, ICONS 

COMMON/  BLKI  / AMP  S • ARCLt  N,  OARC,l)CONSf  OTHRT  ,£NT  H,  PCCNV  , PGAS, 

♦ PGASEN,  PIN,PHAO,PRES,  (iCCNViORAOtSIGMA,  T ARC,  TCAS, WDOT.Y 
COMMON/  BLK2  / ON,TSTAG 

COMMON/  fRROR  / I ARCO, I ARC T , I ERR , I SCNI C , H G AS 
COMMON/  TITLE  / TI TLE ( 18)  ,OAT I4 ) 

COMMON/  CONVF  / 

♦CNVFOl , CNVFO?  CNV^03 , CNVC04, 

♦CNVF05,  CNVF06,  CNVf 07,  CMVF08,  PI 
COMMON/  OAO  / KNUP,KNU 
COMMON/  LENGTH  /CONSLN ,F ARC LN , R ARCLN 
COMMON/  NCASE  / NCASE 
COMMON/  RET  / RE  Y,  WOOTAR  ,WOOTCA  ,<=NU 
iour-6 

GPAO-I  AMPS/C  SI  CHA*OAP'’  *OARC  *P I /4 . 0 ) ) *CNVf03 

LCCNS«C0NSLN*CNVF01 

FPNT*FARCLN*CNVF01 

REAP«RARCLN*CNVF01 

LENX«APCLEN*CNVFOl 

VOLTS-GRAO*LEMX 

: HELIUM  ROOM  TEwp  FMHALPY  * 670  RTU/Lb 

HIN-I6  70./CNVF  04|*WOOT 
EFF»((PCAS  - HIN)/PIN)*100. 

WD0TX«WD0T*CNVF02 

0RA0X-QRAt>*CNVF06 

UC0NVX-0C0NV*CNVF06 

CNX*QN*0.948E-3 

ENTMX-fNTH*CNVFO‘t 

DIARCX«DARC*CNVF01 

OCONSX*OCQNS*CNVFOl 

OTHX-OTHP T*CNVF0I 

PINX-P IN*0.001 

PGASX-IPGAS  - H!r(I*0.om 

HREARX»QRA0*PI*I0CDNS)  *kARCLN*0-R-*8F-j 

HFRNTX*(ORAO*PI*OCONS*FARCLN  ♦ QCONV *P I *DCONS ♦CONS LN I • .948E-3 
IF(NCASE.NE.C)  go  to  100 
WRITE! inUT,l05) 
lOS  FORMAT! IHI,//|IX 

♦61HHCARC  HELIUM  MINIMUM  ENERGY  ARC  MfATFR  CHARACTERIZATION  CnnC) 
WRITE!  inuT,110IDAT, TITLE 

100  NCASf*l  , 

no  FORMAT!  IX,  4AA,//U18A-,  //,^M  POrS  ,6H  AMPS,6H  OCONS , 

♦ 6M  0THRT,6H  ICONS, 6H  APCL,6M  RPAR,bH  TSTAG, 

f6H  OKEAP.oH  UFRNT, t>H  JNJ2, 

♦ /,6M  ATK,6X,6!4X2MIN),6H  DEG  K,IRH BTU/SEC ,//, 

♦ 6M  VOLTS, 6H  PIN,6H  GASP,6H  CRAt),6M  0KAC,6H  OCCNV, 

♦ 6H  WDI«T,6M  INTH,6H  £rr,^H  TAPC,6H  0ARC,/,tX,6M  KW, 

♦6M  KW,6H  V/IN,12H  0TU/TT**2-S  , fcH  L6M/S, 

♦6H  b/LB,6H  PCRC.6M  OEG  K,6H  IN,/) 

IF! 1 ERR ,E0. l)hP I TH IDUT, 126) I SUNIC, f TGAS,  I ARCO,  I AECT 
170  FORMAT! //,  »ERPOR  1S)N1C*M4*  ITCaS*M<.»  |ARC0*M4 

♦ • IARCT-M4  J 

WRITE!  1UUT,130)PRLS,AMPS»I/CCNSX,OTHi(,  LCONS  , L C'JX,  f RNT  , REAR , 

♦ tstaG,hrearx,htrntx,cnx 
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CARO  NO 

TABLE  A-7  (CONTINUED) 

• •••  CCNf£S't$ 

130  FURHATI/,2l6t2F6.3*6F6.1.6j  6) 

33^ 

WRITE!  10UTtK0)V0LtS,PINX.PCASXfCRA0,QRA0X>0C0NVX, 

>wpotx,enthx,eff,tarc,durcx 
160  FURMAT(6I6t  F6.3.3 (6 ,F6. 31 

333 

33A 

335 

RETURN 

336 

END 

337 

c 

338 

c 

339 

SUBROUTINE  GAHHAXI TEMP iPRhStCAHHA  ) 

3A0 

GAMHA-1.6667 

361 

RETl^N 

362 

END 

363 

c 

366 

c 

365 

c 

366 

c 

367 

c 

368 

SUBROUTINE  TEHPIH.P.TJ 
GIVEN  H.P  FIND  T USING  ITERATION  SCHEME 
T IN  JEG  K 

369 

c 

350 

c 

351 

c 

P IN  ATH 

352 

c 

M In  joule S/KC 

353 

c 

TI  INITIAL  TEMPERATURE  GUESS 

356 

TI-5000. 

355 

tn»t; 

356 

100  CONTINUE 

357 

EkR-.OOOOl 

358 

K-0 

359 

110  CONTINUE 

360 

K»K  + l 

361 

CALL  ENTHXlTNrPfHTNI 

362 

ot*,oi»tn 

363 

CAIL  ENTHXI TN«OT,PtKTNP) 

TNP1«TN^(H-HTHI /! IHlNP-HTNI /OT  I 

366 

365 

IF  IK.GT.1000.1  GO  TO  120 

366 

CONV«! tnpi-tni/tn 
CONV»A8S(CONVI 

367 

368 

TN»TNF 1 

369 

IF  ICONV.GT.ERR)  go  to  110 

370 

T«TNP1 

371 

RETURN 

372 

120  CONTINUE 

373 

ENO 

376 

c 

375 

SUBROUTINE  VISCX(T,P,TVISC) 

376 

c 

HELIUM  VISCOSITY  I-IOOO  ATM.  TO  30.0000EG  K 

DIMENSION  TC(  6,311  ,PR 461  ,ATMH30*  #ATM1 0(301  ,*TM100C 

data  pr/1.0, 10.0,100.0,1600.0/ 

377 

378 

379 

CATA  lN/0/ 

380 

CATA  ATHl 

381 

+/.6^C9E-03t  .6957E-03,.96l*.E-03,.l2l2E-02..1676E-02 

382 

♦.1729E-02,*1980F-02,.2226F-02,.2668£-02,.2705F-02, 

383 

♦•2933E-02,.3l37E-02,.3290F-02f.3362E-02,.3230E-02, 

386 

♦ .292  3fc-02,.2652E-02,.l9l6F-02,.  1609E-02  ».99‘t3E-03, 

385 

♦ .687-.E-0  3,.67l  7E-03,.3268E-03, ,2 3886-03,.  1869C-03, 

386 

♦.15^66-03,. 1306E-O3,. 131 8E-03,.l307t-03,.1330E-03  - 

337 

DATA  ATM10/.6-.09E-03, 

388 

♦.6957£-03,.9tl6£-03*.1212E-02,.l67AE-o2, 

389 

♦.1729E-02».l9B0E-02,.2226E-02,.2668E-O2t.2707E-02» 

390 

♦.2960E-02».3166E-02,.3369E-02,.3537E-02,.3636E-02, 

391 

♦. 36286-02,. 3677E-02,. 3 1066-02,. 2 7686 -02 f.2306E-O2, 

392 

♦ .H66E-02,.  163  2E-02,.1099E-02,.8356E-03,.  63806-03, 

393 

♦.6966fc-03,.3933E-03,.325lE-03,.2802c-03,.2522E-03  , 

CATA  ATHl00/.‘»609E-03, 

396 

39  5 

♦•6957|-03,.9616E-03,. 121 2E-02,. 1676 E-02, 

396 

♦.l729E-02i .19806-02,, 222 66-02,. ?66BE-02,.2707t-02. 

397 

♦. 296 2E-02,.3173F-02,. 33956-02,. 36066-02 #. 3786 E-02» 

398 

♦ . 39266-02..  39')7E-02t.3982F-02,.  3 856E-02,.  36 196-02, 

399 

♦.3288E-02* . 2390E-02,.26B56- 02,. 2089E-02 ,.17316-02, 

600 

♦ . 16  236-02,. 1168E-02,.9630E-03,.  802 16-03,. 67786-03  , 

CATA  ATl000/.6409E-03f 

601 

602 

♦.6957E-03,.96l6E-03,. 121 26-02.. 16766-02, 

603 

♦ . 17296-02,.  l9aOE-32,.?22t)6-O2,.2660F-O2,.  27076-02, 

606 

♦.29636-02, .3176E-02,.3606E-02,.3626fc-02,.3a38E-02, 

605 

♦ . 603 26-02,. 62016-02,. 6 330E-02,. 6403E-02  ,.6..0t>E-02, 

606 

♦ . 6:  276-02,.  6162E-02,.?  91  7E-02..  3609'^ -02,.  3 2656-02, 

607 

♦ . 2)066-02,.  25516-02..2213E-02,.  l92*.E-02  ,.i665E-02 

608 

1F( IN.NE.O)  GO  TO  <20 

609 

00  100  J«l,30 
TCi  l,J^ll«ArHlU) 

610 

611 

TC(2,J^1)*ATM10(JI 

612 

TCI  3. J^n-*ATHIOO(J» 

613 

TC(6,J  + n-AT1000(J) 

100  continue 

616 

615 

00  110  1-1.6 

616 

TC( I|l)-0.0 
110  CONTINUE 

617 

618 

lN-1 

619 

120  CONTINUE 

6 20 

IF(P.GT.PR(  3M  1-3 

621 

|F(P.L6.P«(3n  1-2 

622 

lF(P.Lt.PP(?n  1-1 

623 

J>A1NT(T/1000.U1 

626 

OT»T-(J-l)*lOOO, 

625 

i)T-OT/1000. 

626 

PP-ALOOlO(P/PPm» 

627 

TCl-TC(I  ,J)*(TC(I,J^n-TC(I  ,JM*OT 

628 

TC2-Tcn*i,j:^(TC(ui,  j>u-TC(i^i,  jn*0T 

629 

TVISC-TCI^(TC2-TC1)*PP 

630 

TVISC-TVISC60.10 

631 

RETURN 

632 

END 

633 

SUBPCiUTINl  LNTHxn,P,TFNTMJ 

636 

c 

HELIUM  enthalpy  i-lOCO  ATM  T)  30.000  OLG  K 

635 

c 

636 

DIMENSION  ATHK2U  .ATMlO(2l»  . A T M100(  21 ) . AT  lOOOl  2 1 ) 

637 

DIMENSION  TC(6^1)  ,PR(6) 

UATA  PR/1. 0.10^, 100.0*1000.0/ 

638 

6 39 

CATA  lN/0/ 

660 

DATA  RME/2077.02/ 

661 

DATA  TO/273.16/ 

• ••• 


ICOOI ?OI 
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TABLE  A-7  (CONTINUED) 


CARO  NO 


CONTENTS 


AA2 


DATA  ATKl  / 


443 

♦.9160E 

02, .10106 

03,.1110£ 

03,.  12236 

03,-13616 

03 

444 

♦.15476 

03, .18116 

03,.2196C 

03, .27556 

03,-35476 

04 

445 

♦ .<-6286 

03, .6002E 

03,.7o22E 

03,. 93316 

03, .10936 

04 

446 

♦.1226E 

04,-13306 

04,.1407E 

04,. 14646 

04,-15086 

04 

447 

♦ .154  3E 

C4  / 

448 

DATA  ATMIO  / 

449 

♦ 

.91556 

02,. 10086 

03,-ll02E 
03,. 17656 

03, .12006 

03 

450 

♦.13076 

03,-l4?9E 

03, .15776 

03,.2011E 

03 

451 

♦.23406 

03,.2773E 

03, .33446 

03, .40806 

02  ,.4988F 

03 

452 

♦.6080E 

03,. 73226 

03, .36526 

03,. 99896 

03.. 11236 

04 

453 

♦.1234E 

04..1327E 

04  / 

454 

OATA  ATMIOO  / 

455 

♦ 

.91536 

02,-10076 

03,. 10996 

03: 

456 

♦.1193c 

03,. 12906 

03,. 13916 

03,. 15016 

03,-16256 

03 

457 

♦.1768E 

03,. 1939E 

03, .21466 

03,. 24056 

03, .27256 

03 

458 

♦.31136 

03,. 35916 

03,.4165f 

03,. 48436 

03,.5&266 

03, 

59 

♦.65C6E 

03,.7465E 

03,.64%2E 

03  / 

460 

CATA  ATIOOO  / 

461 

♦.9152E 

02,.1007E 

03, 

462 

♦.1099E 

03,.1191E 

03,.1284E 

03,-13796 

03,. 14776 

03 

463 

♦.1579C 

03,. 1669E 

03,. 18086 

03,. 19416 

03, .20916 

03, 

464 

♦-2265E 

03,.2466E 

03, .26976 

03,.2970E 

03, .32876 

03 

465 

♦.36556 

03,.4063E 

03. .45566 

03, .50936 

03  / 

466 

467 

468 

469 
4TO 

471 

472 

473 

474 

475 

476 

477 

478 

4 79 
C80 
461 

482 

483 

484 

485 

486 

487 

488 

489 

490 

491 
*»92 

493 

494 

495 

496 
4Q7 

498 

499 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 

510 

511 

512 

513 

514 

515 
5*6 

517 

518 

519 

520 

521 

522 

523 

524 

525 

526 

5 27 

528 

529 

530 

531 

532 

533 

534 

535 

536 

537 
536 

539 

540 

541 

542 

543 

544 

545 

546 

547 

548 

549 

550 

551 

552 


1E(  IN.NE.O)  CO  TO  llO 
00  100  K-1,21 
TC(  1,K)-ATH1CK) 

TC12,K)»AT«10(K) 

TCI3,KJ*ATM100(K) 

100  TCUtK  )*ATIOOO(K) 

IN*1 

110  CONtlNUE 

IFIT.LT. 10000.)  CO  TO  120 
IFtP.GT.PROM  1^3 
IF|P.LE.Pk(3)J  1*2 
IF(P.LE.PR12»)  I«1 
J«AINT(T/1000.)-9 
0T«T-( J^9)* 1000. 

OT»OT/lOOO. 

PP*AL0G10(P/PP (I ) I 
TCl»TCn,J)*(TCn,J»l)-TCM  ,J»  »*0T 
TC2-TU  Ul»J)^acll*l.J»U-tC(I^l»JI  » 
TENIH-7C !♦( TC2-TC1)*PP 
TENTH«TcNTH*RHE*TO 
RETURN 

120  TENTM«2.5*T*RME 
RETURN 
END 

subroutine  RH0(T,P,T0ENS) 

MELIU?'  DENSITY  1-1000  ATH. 


TO 

DIMENSION  ATMK21)  .ATM10I21I  tATMlOOl 
DIMENSION  T0(4,21l .PR14! 

DATA  PR/1.0, 10. OtlOO.OtlOOO.O/ 

DATA  IN/0/ 

DATA  TO, RHOO/273. 16, 0.17359/ 

CATA  ATHl  / 

♦ .27  2 IE-01,. 2483E-0l,.2274E-0l,.  20956 

♦.1795F-01,.lo59E-0l ,.1525E-01 ,.1389E 
♦.1110E-0l,.S758E-02,.8536E-02,. 750O£ 
♦.6097e-02,.5647E-02,.5304£-02,.5033E 
♦.4619E-02  / 

DATA  ATHIO  / 

♦ .2M2E  OO,.2403E  00,.2276F 

♦ .I'M.'e  00, .18126  00,.1691E  00,.1579E 

♦.1369E  00,.1267F  00,.1165E  00,.1063E 
♦.E671E-OI,.77^0E-Ol,.70llE-Ol,.6348fc 
♦. 53696-01, .5018E-01  / 

CATA  ATHIOO  / . 

♦ .2732E  01,.2483E 

♦ .2101E  01,.1950E  01, .13186  01  , .17026 

♦ •IS02E  0U.1414E  01, .13316  01,.1251E 

♦•1097E  01,.1022E  01, .9^866  00,.b770E 
♦.7A29E  00,.6825E  00,.6276E  00  / 

CATA  ATIOOO  / 

♦.2732E  02,.2463E  02, 


30,000  CEG  K 
21)  ,ATIOOO(2P 


-01,,19‘'8£“01, 
-01, .12616-01, 
-02,.6703f“02, 
-02,.4810E-02, 


00,.2099E  00, 
00, .14736  00, 
00, .96206-01, 
-01,.«806E-01, 


01  , .22766  01, 
01  ,.15976  01, 
01  , .11736  01, 
00,.8080E  00, 


.2276E  02, .21016 
'2, .15126 


02.. 1951E  02,.1320E 

02. . 14296  02,.1353E 


♦ .1604E  02 _ _ _ 

♦ .1216E  02..H53E  02, .10936  02,.1035E 
♦.9226E  01,.8682E  01,-81636  01,.7656E 

1F{ IN.NE.O)  GO  TO  110 
DC  100  K«l,21 
TD(  l,K)«ATHl(Kt 
T0(2,K)«ATH10IK) 

T0( 3,K)-ATM100(K) 

100  T0(4,K)«AT1000|K) 

IN«1 

110  CONTINUE 

IFir.LT. 10000.)  GO  TO  120 
IF|P.GT.PR13)I  I»3 
IF(P,L6.PR(3))  1-2 

IF(P.L6.PR{2n  1*1 
J*AlNm/1000.)-9 
0T«T-IJ*9l*1000. 

DT.OT/1000. 

PP-AIDG10<P/PR  MI) 

T01»TO(  I ,J)^(T0II  .J^U-TOn  ,J)  I *0T 
T02«T0IU1  • toll  ♦!,  J*  l»-ton  ♦!  , J)  i*or 
TDENS*TOl*l  TO?“TOU*PP 
TDENS^  m^lNS*P*i00 

return 

120  TO£NS'P*(TCM)*»<HaO 
RETURN 
END 

SUBPOUTINC  SICMAXI T,P,TSIGJ 
HELIUM  electrical  CCNOJCTIVITY 


02  ,.17056  02, 

02.. 12B3E  02, 

02.. 07t!0£  01, 
01  / 


l- 

DIMENSIUN  ATMH2  5)  ,ATM10(25)  ,ATM100I 
DIMENSION  TS(4,26) ,PR(4) 

DATA  PR/1.0,10.0,100.0,1000.0/ 

CATA  IN/0/ 


lOOOATM.  TO  30,000  060 
25)  ,AT1000(2‘‘» 
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TABLE  A-7  (CONTINUED) 


CARO  KO 


CCNIENTS 


553 

554 

555 

556 

557 

558 

559 

560 

561 

562 

563 

564 

565 

566 

567 

568 

569 

570 

571 

572 

573 

574 

in 

577 

578 

579 
500 

III 

583 

584 

585 

586 

587 

588 

589 

590 

III 

593 

594 

595 

596 

597 

598 

599 

600 
601 
602 

603 

604 

605 

606 

607 

608 

609 

610 

612 

613 

614 

615 

616 
6p 
616 

619 

620 
621 
622 

623 

624 

625 

626 

627 

628 

629 

630 
6>l 

632 

633 

634 

635 

636 

637 
636 

639 

640 

642 

643 

644 

645 

646 

647 

648 

649 

650 

651 

652 

653 

654 

655 

656 

657 

658 

659 

660 
661 
662 
663 


DATA  ATMl  / 

♦•22 34E-04,*7520£-03i. 10706-01 ,.a643£-0l . 
♦•20716  01i«71026  01fl77  56  02  f31**ee  02 » 

♦ •55556  02», 65746  02f7544E  02f8482E  C2. 

♦.1027E  03,.1109E  03fll86E  03fl256E  03» 

♦.13896  03, .14526  03,. 15136  03,.1573t  03, 

DATA  ATM 10  / 

♦.7065g-05,. 23806-03, .33856-02, .27176-01, 
♦•60726  00, .20656  01,. 59866  01,. 14^56  02, 

♦•44846  02, .61496  02,.76796  C2,.9tl4E  02, 

♦.11766  03,.1303E  03..1424E  03f.l540E  0>, 

♦.1753E  03..1650E  03..1941E  03..2026E  05, 

DATA  ATM 100  / 

♦. 22356-05,. 7534E-04,.1073E-02,. 86336-02, 
♦•1882E  00,.6162E  00, .17476  01, .43806  01, 
♦.19636  02,. 34856  02, .54096  02,. 75836  02, 
♦•11976  03,. 14076  03, .16056  03,. 17986  03, 
♦•21616  03,. 23326  03r*2495E  03, .26486  03, 
DATA  ATIOOO  / 

♦.75006-06,. 23876-04, .34076-03, .2 7526-02, 
♦•6011E-01,.1960E  00,.5397E  00,.1320E  01, 

♦•6064E  01, .11496  02, .20556  02,.342SE  02, 

♦•77556  02,. 10616  03,. 13766  03,. 17056  03, 

♦.23736  03,. 27016  03,.30276  03..3329E  03, 

IFUN.NE.OI  GO  TO  150 
DO  100  N-1,25 
TS( l,S)«ArHl(NI 
TS(2,N)«ATM)0(N) 

T$( 3,N1-ATM100(NI 
100  T$(4,N}*AT1000(N) 

110  00  140  K«l,4 
00  120  L«l,18 
M-26-L 

120  TS(K#M^l)-:S{K,M) 

00  130  L-1,8 

130  TS(K«l.;>AL0C10(TS(K,Ln 
140  CONTINUE 
tN-1 

150  CONTINUE 

I6(T.LT.6000.) 


4.793E  00, 
4434E  02, 
9395E  07, 
13256  03, 
1633E  03 

1474E  00, 
2321E  02, 
10466  03, 
16506  03, 
2108E  03 

46506-01, 
98536  01, 
9799E  02, 
1983E  03, 
2792E  03 

14066-01, 
2959E  01, 
53246  02, 
20396  03, 
3624E  03 


GO  TO  180 
1»3 
I«2 

I-l 

JSHI6T-5 

IFIT.Ce.13000.1  jSHirT*4 
J«AINT(T/i000.)-JSHlF7 
OT-T-(J^JSMIFTM1000. 
07-07/1000. 

160  PP-ALOGIOIP/PRU  ») 


IF«P.C7.PRt3I) 
IFIP.LE.PftOU 
IFCP.L6.PR(2)I 
FT«5 


TSl-7SCI,JI^(TSn,J^l)-TSn  ,JM*OT 
TS2-TSI  I^I.JI^ITSIUI,  j^i)-tsn^i,jn*OT 
TSIG-TSI^ITS2-TSU*PP 
IFiT.LT,13000.»  GO  TC  170 
TSlC*TSlG*lOO.O 


RETURN 
170  TSIG-10**TS1G 
TSIG-TSIG*100.0 
return 

180  TSlG-0.0 
RETURN 

IuDROJTINF  RAUXlT,P,TRAOI 


,ATM10(2l) ,ATM100(21) ,Ari000(21) 
PRI4I 


PEAL  NE.NT 
DIMENSION  ATHK2 
PIHENSION  TRU,2^, 

CATA  PR/1.0, 10.0, 100.0, 1000.0/ 

DATA  IN/0/ 

CATA  ATMl  / 

♦ . 73 10E-0t,.3073E-03,. 1 01 66-02,. 781 7C-02,. 67786-07, 

♦ .14576-01 , .28426-01 , .50926-01 ,.  8«'^36-01  ,.13016  00, 

♦•18706  00,.2499E  00,.3126t  00, .36786  00,. 41126  00, 
♦.4420E  00,.A626E  00,.4758F  00, .48416  00, .48946  00, 
♦.49206  00  / 

DATA  atm  10  / 

♦ . 23666-04, .9809F-04,. 32566-03, .90746-03, 
<•2 19 *6-02, ,4 7626-02,. 9419F-02,. 1 7296-01,. 29436-01, 

♦ . 43516-01,. 72136-01  ,.10-*4E  00,. 14426  00,. 18936  00, 
♦.23626  00,.2871E  00, .33276  00, .37266  00, .^0496  00, 
♦.4301E  00..44Q0E  00  t 

DATA  ATHIOO  / 

♦ . 75566-05, .31A6E-04, .10456-03, 

♦.29426-03,.7l87E-03,.l5’2E-02,.3l20E-02,.5804£-02, 
♦.10016-01,. 164 16-01 ,.256 16-01 r. 38276-01,. 55 19E-01, 
♦•76C76-Oi. .10176  CO, .13176  00, .16526  00,.7014E  00, 
♦•2386E  00*.2755E  00, .31086  00  / 

DATA  ATIOOO  / 

♦.24296-05, .10166-06, 

♦.34136-04,, 95816-04,. 2360F-03,. 52526-03,. 10596-02, 
♦. 19656-02,, 3A51E-02, .57386-02,. 91016-02,. 1303 6-01, 
♦. 20296-01, .28606-01,. 39656-01, .5302F-01 ,.69156-01, 
♦•880^6-01. ,:i006  00,.1334£  00,.1592E  00  / 

IFt  »N,N£.0)  00  TO  130 

00  100  N*l,21 
TKl L,N)«ATM1(N) 

TRI 2,N)«ATM10(N1 
7R(  3,N)*ATH100IN) 

100  TK(4.N}*AT1000(N1 
00  l50  R'*l,4 
no  L»i,r» 

M-22-L 

no  TPlK.M^n«TklK,M,' 

00  lio  L-1.7 
TR(K,L)*AL0G10( TP(K,U ) 

120  continue 

1N»1 

130  CUNTINUE 

1FIT.LT.  10000.)  GO  TU  160 
IF(P.Gr.PR*3)l  1*3 
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TABLE  A-7  (CONCLUDED) 


CARO  NO 


CCNT6NTS 


664 

665 

666 

667 

668 

669 

670 

671 

672 

673 

674 

675 

676 

677 
6T8 

679 

680 
681 
682 

683 

684 

685 

686 

687 

688 
669 

690 

691 

692 

693 
69^ 

695 

696 

697 

698 

699 

700 

701 

702 

703 

704 

705 

706 

707 

708 

709 

710 

711 

712 

713 

714 

715 

716 

717 

718 

719 

720 

721 

722 

723 

724 

725 

726 
7 27 
728 
7 29 
7 30 

731 

732 

733 
^'34 

735 

736 

737 

738 
7 39 

740 

741 

742 

743 

744 


IF(P.LE.PRO))  1-2 
IFJP.LE.PP(2))  I-l 
JSm!FT»9 

IF( T.CF. 16000, ) JSHIFT.8 
j»a;n7it/iooo,)-jshift 
Or«T-( J+JSHIFT»*1000, 

CTi>DT/1000. 

140  PP»ALOClO(P/PR(n) 

TRl-TRU  f J)^(  TRU  f J*l)-TRJI  .JU^OT 
TR2-TR(Ul,J)»(TRiUl,>l)-TR(I^l,  Jl  )*0T 
XE  ■TRlt(7R2'TRl)*PP 
IFCT.GF. 16000,)  GU  TO  150 
X6»10**X6 

150  NT»(P*1. 01326  5/U,  386-23*T|)  •UE-6 
NE-X6*NT 

TRA0-(4.34E-28)*NE**2/SQKTm 

TRA0»TRAD*1,0E^6 

RETURN 

160  KP»0,1 352* T** I 5/2)^EXP (-285500, /TI 
KPP*S0KT(KP/(  1.0132E  5*PM 
XE»XPP/( 1.0+KPP) 

GO  TO  150 
END 

SUBROUTINE  TCONDXC TrP.TCONOI  ^ 

HELIUM  thermal  CONOUCtiVITV  I-IOOO  ATM.  TO  3O,O0ODEC  K 
9IMSNSI0N  TC(<-f  31)  ,PR(4) 

DIMENSION  ATHlOO)  tATM10(30)  tATMlOOOO)  iAT1000(30) 

DATA  PR/l.0tlC*0tl06.0, 1000.0/ 

DATA  IN/0/ 

DATA  ATM  I / 

4,3<.C5E  05».5378E  05. .72776  05..9366E  05,,1139e  06, 

♦.13376  06..1530E  06..1721E  06..1910E  06,.2104£  06, 

♦.2321E  06, .26006  06, .29856  06,.3459E  06, .39906  06, 

♦ .tbllE  06, .54286  06, .65566  06,. 80116  06, .96186  06, 

♦.1098E  07..1165E  07,,ll36£  07,.1030E  07,.8982C  06, 

♦ .7858E  06, .71056  0t>,. 67106  06, ,65936  06..6670E  06  / 

DATA  ATHIO  / 

♦ .3408E  05, .53786  05, .72776  05,. 93666  05,.1139F  09, 

♦ .1337E  06,.153C:  06, .17216  06,.1909E  06, .20966  96, 

♦ .2289E  06,-25036  06,.2764C  06,.3109E  06, .35576  06, 

♦.4079E  06, .^6456  06,-52636  06, .59886  06,-68776  06, 

♦ .7957E  06, .91866  06,. 10466  07,.U61E  07, 

♦ .120*6  07,-12786  07, .12416  07,.U90E  07, 

DATA  ATMIOO  / 

♦.34C8E  05,.5378E  05,.7277£  05,-03666  05, .11596  06, 
f.l337E  06,.1530r  06,. 17216  06,. 19086  06, .209* 6 06. 

♦.2280E  06, .24716  06..2676E  06..2913E  06..321CE  06, 

♦.3593E  06,-4084E  06, .46676  06,-53206  06, .60356  06, 

♦-6320E  06, .70716  06,.8645F  06..9739E  06, .1093-  07, 

♦ .1217E  07..1341E  07, .14576  07..1558E  07,. 1639k  07 


-U61E  07,-12456  07, 

- 07 


DATA  ATIOOO  / 
♦*34086  05, .53786  05, 


.11396  l6, 


..72776  C5,. 93666  05, _ 

♦.13376  06,.1530E  06, .17216  06,.1908E  06, .20936  0^  . 
♦ .22776  06, .24006  06,.2648E  06, .28456  06, .30626  06 
f.  33146  06, .36226  06*. 40146  06, .45136  06, ,.51416  06, 
♦.5902E  06, .67996  06, .782^6  06, .88826  06, .10116  07, 
♦.lltSe  07f. 12916  07,. 14436  C7,. 16086  07, .17746  07 
1F(  IN.N6.0)  GO  TO  l20 
OU  100  N»l,30 
TCJl,Nfl)=ATMl(N) 

TC(2,N^n»ATHlO(NI 
TC(  3,NfU*ATMlOO(N) 

100  TC(4,N4l)«AT1000(N) 

DO  no  1*1-4 
TCI [*1)-0.0 
llO  CONTINUE 
IN-l 

120  CONTINUE 

IFIP.GT.PROn 
IF(P.LE.PP(3n 
IF(P,LE.PR(2n 
J-AINT(T/1000.)^l 
OT*T-( J-1)*1000. 

OTaOT/IOOO. 

PP*AL0G10(P/PR  (in 

TC1«TC(  I iJI^ITCn  ,JtU-TC(I  ,J)  )*0T 
TC2*TC(ni,J)  + (TC(l  + l,  J*U-TC(!  MfJ)  l•OT 
TC0ND-TCl*(TC2-TCl)*PP 
TCnNn-TC0N0*l,0E-9 
RETURN 
END 


1*3 

1-2 

I»l 


118 


AEOC-TR.76-25 


iuTCft’W  CMlt*  StI  - 


CHlKI  TIUl  - rOCfJUkES 


I Nort  O'" 

'i  NOU  0» 

*•••••••••• 

• HfCIN  00  ♦ 

• 160  16  • I*  • 


I .♦  ‘ 

I - n 1 


MtLiOk  N-?*D  Mf6l€k 

i 


I |P*CS/lO./(6*  ! 

I f*fs  ♦ ! 

I t,2?M*CHVfO)  j 

1 P*«CVN  • <••.6^  ♦ I 

I 0.0t06*l“RCS  - j 

I KA6ft6?CNVfO)l  I 

I •CNVF05  1 


Figure  A-1  HEArC  flow  chart 
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Figure  A-1  H E ARC  flow  chart  (Continued ) 
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Figure  A-1  HE/  kC  flow  chart  (Continued) 
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Figure  A-1  HEARC  flow  chart  (Continued) 
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Figure  A-1  HEARC  flow  chart  (Continued) 
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Figure  A-1  HEARC  flow  chart  (Continued) 
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Figure  A-l  HEARC  flow  chart  (Continued) 
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Figure  A*1  HEARC  flow  chart  (Continued) 
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Figure  A-1  HEARC  flow  chart  (Continued) 
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Figure  A-1  HEARC  flow  chart  (Continued) 
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Figure  A-1  HEARC  flow  chart  (Continued) 
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Figure  A-1  HEARC  flow  chart  (Continued) 
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Figure  A<1  HEARC  flow  chart  (Continued) 


I 


131 


AEDC-TR.76-25 


07/U/7S 
CMIM  TlUt 


su^KOuTiNF  siCHAxi r,^,rsici 


AUtO^LOM  CHAAT  s»1  - 


JA 


oz 

"l 


J • 

AINT  V|C00.I  “ 


OT  • (JT/IOOO. 

■ i 

pp'I"“  I 

AL);«10(P/P(«f  m ( 


T$n,j)  * I 
j ♦ h - I 
J)l*OT 


I TS2 

• T$(|  f 1 

liJl 

* (fSfl  » I 

1 ItJ  » 

n - ISM  * 1 

j ' 

iJ>M0T  j 

■'ll 

• TSl  ♦ i 

I 

- TSIMPP  1 

i 

• 1 -tT.  IJOOO.  •" 


I > TsiG*too.u  i 

_ * 

I 



•*"exFt  • 


>— — — »• 

^ I’O  ! If* 

I ”t  sT~“  10**  TSi  r,’“  j 

I T$|f.  . T^if,*ioo.o  I 

• 

i 

n 

• <X|T  < 


Figure  A-1  HEARC  flow  chart  (Continued) 
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Figure  A-1  HEARC  flow  chart  (Continued) 
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Figure  A-1  HEARC  flow  chart  (Continued) 
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NOMENCCATURE 

Cross  sectional  area 

Radiation  constant,  mass  flow  constant 

Diameter 

Gravity  constant 

Electrical  voltage  gradient 

Enthalpy  Ccomputer  codes) 

Turbulent  heat  transfer  coefficient 

Enthalpy 

Arc  current 

Arc  length 

Mass  flow  rate 

Arc  heater  scale 

Electron  density  in  arc  column 

Nusselt  number 

Pressure 

Power 

Prandtl  number 
Gas  constant 
Reynolds  number 
Heat  loss 
Heat  flux 
Temperature 

Arc  column  volumetric  radiation 
Arc  voltage 

Gas  velocity  in  constrictor 
Water  flow  rate 
Input  lower  parameter 
Radiation  parameter 
Convection  paramei:er 
Ratio  of  specific  heats 
Thermal  conductivity 
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n 

P 
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a 

A 

Subscripts 

a 

a 

A 

b 

c 

C 

F 

g 

H 

He 

i 

N 

N-250 

o 

r 

R 

sf 

w 

X 

Superscript 

* 


Efficiency 

Viscosity 

Density 

Electrical  conductivity 
Differential 


Arc  column 

Air  (mass  flow  rate) 

Anode 

Bulk  gas  at  nozzle 
Constric'  r or  convection 
Cathode 

Front  electrode 

Bulk  gas  around  arc  column 

Hydrogen 

Helium 

Input  injector  or  value 
Nozzle 

N-250  arc  heater 
Stagnation  condition 
Radiation 
Rear  electrode 
Sonic  flow  value 
Wall  value 

Weighted  conductivity  - temperature 


Nozzle  throat  value 
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